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TRANSLATION EDITOR’S PREFACE 


The first, Russian, edition of Linear Differential Operators by N. A. 
Naimark appeared in 1952. The English edition of this important book has had 
the benefit of advice from, and additional material supplied by, Professor 
Naimark, who has also written a preface for this new edition. 

Some explanation may be called for in presenting at this time an English 
edition of this book. Since the original Russian edition, several presentations 
of the theory of linear ordinary differential operators have appeared, notably 
those of E. C. Titchmarsh (second edition, 1962), Coddington and Levinson 
(1955), K. S. Miller (1964), G. Hellwig (1964), and Dunford and Schwartz 
(1964). Some of these authors, notably Dunford and Schwartz, use methods of 
functional analysis; others, in particular, Titchmarsh, restrict themselves to 
classical analytical methods; none combines these two approaches so effectively 
as Naimark does in this book. There is no better introduction to the wide range 
of ideas and techniques required for the study of linear differential operators 
than that given here. . 

Part I presents a complete account of the so-called elementary theory of 
linear differential operators, in which the methods of functional analysis are 
kept in the background whilst use is made of the properties of functions of a 
complex variable. There is a detailed investigation of the general, regular, 
eigenvalue problem for ordinary differential equations with complex co- 
efficients and regular boundary conditions. No comparable account exists else- 
where in the literature. 

Part II of the book, which is to be published separately, begins with a com- 
pact account of the theory of linear operators in separable Hilbert spaces; 
these results are then used to present the propcrties of symmetric operators 
derived from ordinary differential expressions with real coefficients. This, in 
turn leads to an account of the spectral theory of such operators and the corre- 
sponding eigenfunction expansions in square-integrable function spaces. Next 
follows a long chapter on the determination of the deficiency indices and the 
spectrum of these differential operators. These results, whilst presented in the 
terminology of Hilbert space, are obtained by the use of methods of classical 
analysis, in particular, the development of certain asymptotic expansions for 
systems of differential equations. There is no comparable account of these 
expansions, which are largely due to Naimark himself, anywhere else in the 
literature. This second part of the book ends with an account of the work of 
Russian mathematicians on the inversion of the Sturm-Liouville problem. 
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v1 TRANSLATION EDITOR’S PREFACE 


The book has been expertly translated by Mr. E. R. Dawson, who has also 
brought the Bibliography up to date. 

Mr. Dawson and I wish to express our thanks to several of our colleagues in 
Queen’s College, Dundee, for their help in reading the manuscript and in 
checking the proofs. 

W. N. EVERITT. 
Queen’s College, Dundee. 
July 1966. 


AUTHOR’S PREFACE TO THE ENGLISH 
EDITION 


This edition incorporates the editorial amendments and additions which were 
made in the German edition of the book at the suggestion of the German 
editors, Professors H. O. Cordes and F. Riihs; the author takes this occasion to 
express his cordial acknowledgment to them. 

In addition, a number of new changes and additions have been made in the 
English edition; these take into account some of the important advances in the 
theory of linear ordinary differential operators which have been made in the 
past few years. In this work the author has been assisted by Professor V. E. 
Lyantse, who has also written a Supplement to Part II. This Supplement deals 
with the author’s work in the theory of non-selfadjoint singular differential 
operators, work which has later been extended by Professor Lyantse’s own 
researches. It is the author’s pleasant duty to express his deep gratitude to 
Professor Lyantse. The author also expresses his gratitude to Professor. W. N. 
Everitt for undertaking the task of editing the English translation of the book. 


M. A. NAIMARK. 
Moscow, May 1966. 
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PREFACE TO THE RUSSIAN EDITION 


Many topics in mathematical physics lead to the problem of determining the 
eigenvalues and eigenfunctions of differential operators and of expanding an 
arbitrary function as a series (or an integral) of eigenfunctions. A problem of 
this sort is encountered, for example, whenever the Fourier method is used to 
find a solution of a partial differential equation under prescribed initial and 
boundary conditions. Consequently a lot of attention has been paid to differ- 
ential operators and they are the subject of much current research. 

Most of the published work deals withthe main problem just mentioned—the 
spectral theory of differential operators, i.e., the investigation of the spectra of 
differential operators and the expansion of given functions in terms of the 
eigenfunctions of these operators. Interest in such problems has been particu- 
larly stimulated by the development of quantum mechanics. The spectral 
theory of differential operators appears as the basic mathematical method for 
investigating many topics in quantum mechanics. In particular, quantum 
mechanics demands a detailed study of “‘singular’’ differential operators, of 
operators, for example, which are defined on some unbounded interval. Such 
operators may, in general, have a continuous spectrum as well as a discrete 
spectrum; and then an expansion in terms of eigenfunctions takes the form of a 
Stieltjes integral. 

It should be mentioned that, in recent years, highly important results in this 
field have been obtained by Soviet mathematicians. Results such as the proof 
of the existence of an expansion in eigenfunctions of a singular, self-adjoint 
operator of arbitrary even order, the solution for the inversion of the Sturm- 
Liouville problem (see Chapter VIII), the completeness of the system of eigen- 
functions and associated functions for a wide class of differential operators 
which are not self-adjoint, and many other results, are due, either entirely, or 
almost entirely, to the efforts of Soviet mathematicians. 

Notwithstanding these fundamental achievements, the problems of spectral 
representation of differential operators cannot yet be regarded as exhausted. 
Particular attention should be drawn to the problem of determining the multi- 
plicity of the spectrum and the deficiency indices of a differential operator in 
relation to the behaviour of its coefficients, and to the problem of selecting and 
normalizing a minimal system of eigenfunctions when the spectrum is con- 
tinuous. In particular, very little work has been done for the case of an operator 
generated by a system of differential expressions, that is, for a differential 
operator in the space of vector-functions (cf. Chapters V and VII). 
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An equally important and interesting problem is that of expanding a given 
function in terms of eigenfunctions of a differential operator which is not self- 
adjoint. A noteworthy advance with this problem was made in a paper by 
M. V. Keldysh [47] for the case of a regular differential expression with, in 
general, non-regular boundary-conditions, but so far no one has investigated 
the case of a singular differential operator which is not self-adjoint. 

Among other outstanding problems in the theory of differential operators 
we may also mention that of putting on a proper mathematical basis the so- 
called “perturbation” theory, which has many applications in quantum 
mechanics (see, in this connection, [39e], [56], [61], [96a], [112b], and [27b)). 

Despite the importance of the theory of linear differential operators, due 
to its applications, and despite the abundance of fundamental results already 
achieved, the literature, both Soviet and foreign, is extremely poor in books 
devoted to this theory. Thus, in the foreign literature there is just one book in 
which a complete and exact account of the fundamental topics in the theory of 
second-order differential operators is given; this is [112a] by Titchmarsh, 
published in 1946. In the Soviet literature there is one monograph by B. M. 
Levitan [63a], published in 1950, in which another, and in fact a very ele- 
mentary, presentation of the theory of differential operators is given. No book 
is yet available dealing with differential operators of higher order. Moreover, 
a whole range of questions arising in the theory of second-order differential 
operators are not discussed at all in the books by Titchmarsh and Levitan. 

In the present book, the theory of linear, ordinary, differential operators of 
arbitrary order is developed. Extensive use 1s made of the ideas and theorems 
of functional analysis, especially those of the theory of linear operators in 
Hilbert space. For the reader’s convenience all the necessary theorems in 
functional analysis are developed in the text itself, to make the book, as the 
author hopes, self-contained and accessible to a wide range of readers. Many 
topics in the theory of linear differential operators can also be discussed without 
using functional analysis, as Titchmarsh’s treatment [112a] shows. However, 
the author does not consider such a treatment to be appropriate, since it 
is only by using the ideas and methods of functional analysis that a deeper 
understanding of the theory, and its most general results, can be achieved. 

The book is in two parts. In the first part, which may be described as an 
elementary theory of differential operators, the use of the methods of functional 
analysis 1s kept to a minimum. This part presents the theory of differential 
operators defined on a finite interval, including the case of differential opera- 
tors which are not self-adjoint, the theory being developed on the assumption 
that the coefficients of the operators are sufficiently smooth, i.e., sufficiently 
differentiable. To understand Part I the reader needs only an elementary know- 
ledge of the theory of ordinary differential and integral equations and of the 
theory of functions. 

Part Il develops the theory of differential operators, using Hilbert space 


PREFACE TO THE RUSSIAN EDITION X1 


methods. Here, in addition to the subjects already mentioned, the reader will 
be expected to have some knowledge of the basic results in the theory of 
the Lebesgue integral. 

Many new results are given in this book, most of them due to the work of 
Soviet mathematicians; in particular, the results relating to the deficiency 
indices and the spectrum of a differential operator. 

The author wishes to express his thanks to M. I. Vishik for reading the 
book in manuscript and making a number of valuable suggestions for improv- 
ing the text. 

M. A. NAIMARK. 
July, 1952. 
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CHAPTER I 


FUNDAMENTAL CONCEPTS AND 
THEOREMS 


$1. Definition and Principal Properties of the Linear Differential Operator 
1. General Definition of a Linear Vector-space and of a Linear Operator 


A set R of elements x,),... is called a linear vector-space if it is constituted as 
follows: 


|". For any two elements x,y €R there is a sum x + y having the following 
properties: 
(a,) if x,y eR, thenx + yeR; 
(bx +yH=y+x; 
(Cy) (x+y) +2=x+(4+2); 
(d,) there is a “‘null element” O in R such that x +0 =x forall xe R. 


2°. For each element x € R and for each real or complex number X there is a 
product Ax having the following properties: 
(a.) if x € R, then Axe R; 
(De) Aux) = (Apn)x; 
(Co) lL.x =x; 
(d,) 0.x = 0 (wherve, on the L.H.S., 0 means the number 0, and on the 
R.H.S., 0 means the null-element); 
(C2) Ax + y) = Ax + Ay; 
(fe) (A + wx = Ax + px. 
The element (—1)x is denoted by —x, and by properties (c,), (f,) and (d.) 


x +(—x) =U +(-—1))x =0.x = 0. 


The elements x,y,... of the space R are called vectors in R. 

In this definition the nature of the elements x,y,... is entirely arbitrary, and 
itis equally a matter of indifference how the concepts ofthe sum x + yand the 
product of an element x with a number A are interpreted in particular cases. It 
is required only that these concepts shall satisfy the conditions enumerated 
above. Any two operations whatever which satisfy these conditions may be 
regarded respectively as addition of vectors, and as multiplication of a vector 
by a number; and the set of elements for which these operations are defined 
may then be called a linear vector space. 

If in the space R multiplication only by real numbers is allowed, then R is 
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called a real space; but if multiplication by any complex number 1s allowed, 
then R is called a complex space. We shall take R to be a complex space unless 
the contrary is specifically stated. An expression of the form A,x, + Agx, +... 
+ X,x, is called a linear combination of the vectors x,,X9,...,X,; a linear com- 
bination is said to be trivial if all the numbers 4,,,,...,A, are zero, otherwise it 
is non-trivial. The vectors x1,X»9,...,X, are Said to be /inearly dependent if at least 
one non-trivial combination of them vanishes, and to be /inearly independent 
if no non-trivial combination of them vanishes. 

A space R is finite-dimensional, or more precisely, n-dimensional if in it there 
are n but not more than av linearly independent vectors. Such a system of n 
vectors is called a basis in R. If, on the other hand, there are arbitrarily many 
linearly independent vectors in R, then R ts infinite-dimensional. 

A subset R’ of the finite-dimensional or infinite-dimensional space R is 
called a subspace’ of Rif every linear combination of elements of A’ is itself an 
element of R’. 

Let Z be a subset in the linear space R. A function A which assigns to each 
element x of ZY some element, x’ = A(x), of Ris called an operator in the space 
Rand @ is the domain of definition of A. Often we write Ax instead of A(x) if 
this is not likely to be misunderstood; and we write GY, instead of Y when it 
is necessary to emphasize that @ is precisely the domain of definition of the 
operator A. 

The set of all vectors Ax,x € ,, is called the range of values of the operator 
A and is denoted by 2, or AD. 

If, generally, & is an arbitrary set, then A& denotes the set of all vectors 
Ax,x € &, 1.e. the set of all vectors which are obtained by applying the operator 
A to all the vectors of the set &. 

An operator A is said to be linear if Y, is a subspace and if, for arbitrary 
vectors x,y € ZY, and for any number J, the relations 


A(Ax) = AA(x), 


A(x + y) = A(x) + AQ) 
hold. 

Two operators A and B in the space R are regarded as equal if and only if 
they are defined on the same domain@ and if Ax = Bx forall x €Q. 

An operator A 1s called an extension of an operator B and we write A > B 
or Bc A if DB, > Dy and if the operators are equal in Y,, i.e. if Ax = Bx 
for all x € Dg. In such a case we also call the operator B a restriction of the 
operator A to Dy. 

In this book we shall consider only operators which are linear, and therefore, 
to avoid constant repetition of the word “linear”, we shall always use the one 
word “‘operator’’, but mean by it “linear operator’. 


1 In the English literature this is usually called a linear manifold. 
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2. Linear Differential Expressions 
A linear differential expression means an expression of the form 


Ky) = po(x)y + pix)y—P +... + p,lx)y. (1) 


The functions po(x),p,(x),...,p,(x) are called the coefficients and the number n 
is called the order of the differential expression. 


In this first chapter we shall assume that the functions Sy PaCsnPalA), 
0 


D,{x) are continuous on a fixed, closed, finite interval [a,b]; in some cases we 
shall impose further conditions. 

We denote by C the set of all functions y(x) which have continuous 
derivatives up to the nth order inclusive on the interval (a,)]. 

For every function y¢ C™ the differential expression ((y) is well defined, 
and represents a function continuous on the interval [a,b]. 


3. Boundary Conditions 


We denote the values of the function y and its first (7 — 1) successive 
derivatives at the boundary points a and b of the interval [a,b] by 


VarVaree VEY; VorVors Vor? (2) 
Let U(y) be a linear form in the variables (2): 


U(y) = % Va + %GVa te. of, yr + Bo¥p + Biyy +... + Bn». (3) 


If several such forms U,(y) have been specified, v = 1,...,m, and if the con- 
ditions 
U,(y) = 0, vy = l,...,™m, (4) 


are imposed on the functions y(x) € C™, we call these the boundary conditions 
which the functions y must satisfy. 

We denote by J the set of all functions ye C™ which satisfy a specified 
system of boundary conditions of the form (4). Clearly is a linear subspace 
in C™ which coincides with C™ only if the conditions (4) are entirely lacking 
or if all their coefficients vanish. 

Suppose a certain differential expression ((y) and a particular subspace Q 
defined by conditions of the form (4) are given. To each function y € QZ we let 
the function u = I(y) correspond. This relation is a linear operator with Das 
its domain of definition; we denote it by L. Using the notation previously 
introduced, we write 

u = Ly. 


The operator L is called the differential operator generated by the differential 
expression I(y) and the boundary conditions (4). 


4 §1 LINEAR DIFFERENTIAL OPERATOR 


In this way—and this fact will be very important in the sequel—one and the 
same differential expression can generate various differential operators depend- 
ing on the way in which the boundary conditions (4) are chosen. 

if, in particular, the conditions (4) are absent, we obtain a differential 
operator, denoted by L,, with the domain of definition Z=C™. L, is 
obviously an extension of all other operators L which can be generated by the 
same differential expression /()). 

For many questions it turns out to be profitable to investigate not merely 
the operator L, which has the widest domain of definition, but also the 
operators considered above, which are restrictions of Ly. 

Certain of the forms U,(y) may be expressible as linear combinations of the 
remaining forms. The corresponding conditions U,(y) =0 are then im- 
mediate consequences of the remaining conditions, and therefore, being 
superfluous, can be discarded. Hence we may, ab initio, regard the forms 
U,(y) as being linearly independent. This implies, of course, that the rank of 
the matrix formed from the coefficients of these forms is equal to m. 

For m = 2n, the equations (4) are equivalent to 


y=yoen. SH yt May = ys. = yO =O. 


The particular differential operator generated by /(y) and these conditions is 
denoted by Lp. 

We can also give the following geometrical interpretation of the boundary 
conditions. Let 7 stand for the whole of system (2): 7 = (VgsV%45-0.5 "7 PLY p, 
Vhyeeey?~), Clearly we may regard 7 as a vector in the 2n-dimensional space 
R,,. The conditions (4) separate out from the space R,, a definite subspace 
M of dimensions (2n — m). These conditions imply that the vector 7 which 
corresponds to a function y belongs to the subspace -@. 


4. The Homogeneous Boundary-Value Problem 


The problem of determining a function y ¢ C” which satisfies the conditions 


(y) = 0, (5) 
U,(y) = 9, y = 1,2,...,71. (6) 


is called the homogeneous boundary-value problem. lf L is the operator which is 
generated by the differential expression /()’) and the boundary conditions (6), 
then the homogeneous boundary-value problem amounts to finding, in the 
domain of definition 2 of the operator L, a function y for which L vanishes. 

Clearly, any homogeneous boundary-value problem always has at least one 
solution, viz. the solution y = 0. This is called the trivial solution. A homo- 
geneous boundary-value problem can, however, have other ‘‘non-trivial’’ 
solutions, i.e. solutions which are not identically zero. 
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We now want to ascertain under what conditions a homogeneous boundary- 
value problem has non-trivial solutions. 

Let }'1,Vo,---,¥, be linearly independent solutions of the differential equation 
iy) = 0. Then, as is known from the theory of linear differential equations 
(see [108] or [17], for instance), any solution of the equation /(y) = 0, and 
therefore the solution of the homogeneous boundary-value problem, too, 
must be expressible in the form 


Y= Cy + Cove + ee BO 


where ¢},C,....C, are certain constants. By substitution in the conditions (6) 
we obtain a system of linear, homogeneous equations 


CU.) + ¢2Ui0'2) + + OUI) = 9 
Cy U,()1) + Co U2(y2) + coe + Ch U(Vn) — 0 (7) 
Cy Un(Vi) ++ C2U (J 2) ++ vee + Cr n(n) — 0 


for the determination of the constants c¢,,Co,...,c,. Let r be the rank of the 
matrix 


Ui.) U2) Ui(yn) | 
U(y1) U2(y2) U(r) 
U= (8) 
U (V1) UY 2) Un Yn) j 


Then there are precisely (n — r) independent solutions of the equations (7) for 
the ¢1,C2,...,C,; and to these correspond (n — r) independent solutions y of the 
boundary-value problem. 

Hence the following result: 
I. If the rank of the matrix U is equal to r, the homogeneous, boundary-value 
problem has exactly (n — r) independent solutions. 

In particular, 

II (a). The homogeneous, boundary-value problem has a non-trivial solution if 
and only if the rank r of the matrix U is less than the order n of the 
differential expression l. 

(b). For m <n, the homogeneous, boundary-value problem always has a 
non-trivial solution. 

(c). Form =n, the homogeneous, boundary-value problem has a non-trivial 
solution if and only if the determinant of the matrix U (in this case, a 
square matrix) vanishes. 


The rank of the matrix U does not depend on the choice of the basic system 


6 §1 LINEAR DIFFERENTIAL OPERATOR 


Vi,-+->¥,- For, the transition from one such system y,,...,y, to another fy,...,9, 
is effected by means of a linear transformation 
XN 


n 
Sp; = \ Aizy jr i= 1,2,...,7, 
j=l 


with a non-zero determinant. Since in this transition the matrix U is multiplied 
by the matrix [a;,] (ij = 1,...,7), the rank of U remains unchanged. 

The rank of the matrix U is called the rank of the boundary-value problem. 

The definition which we have given of the homogeneous boundary-value 
problem may be generalized as follows. 

Let U,(y),U2(y),...,U,(y) be linearly independent functionals in C™. (As 


the norm in C™, we can take |y| = )) max |y“(x)|.) Let /(y) be a differen- 
k=0 aSxSb 


tial expression of the mth order. 
Then by a generalized homogeneous boundary-value problem we mean the 
problem of determining a function y ¢ C™ which satisfies the conditions 


Ky) = 0, U,(y) = 9, vy = 1,2,...,n. 


The results of this section will clearly go over completely to the case of a 
generalized boundary-value problem. 

Certain special cases of the generalized boundary-value problem are treated 
in Tamarkin’s articles [110] and by other writers (see [106] and [45a], Part IJ, 
§5 in which a detailed list of references is given); but the question has not so 
far been discussed in the full generality given here. 


5. Lagrange’s Formula. The Adjoint Differential Expression 


We now assume that the coefficients p,(x),k = 0,1,2,...,1, of the differential 
expression 


Ky) = pax) &% + ry 


dxn-t 


+... + py(x)y 


have continuous derivatives up to the order (n — k) inclusive on the interval 
[a,b]. Further let y and z be two arbitrary functions in C”. By & partial 
integrations we get 


b - th - _ -\/ _ 
|, Pn—n2y Pax — | Pane" vp (Pn—KZ) yp 2) + cee 
_ _ _ x=b _ fb _ 
+ (DEM BPE + (DEY? rp, Dax. 9) 


(Here z denotes the complex number conjugate to z, and Z = z(x) denotes the 
function whose values are conjugate to those of z(x).) 
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If we put k =n, n —1,...,0 in (9) and add the resulting equations, we 
obtain the formula 
b _ b ——— 
|, Won2dx = PO.) +f" yF@az, (10) 
where 
(2) = (—1)(Boe) + (= IBZ) + (= 1a) 
+... + p,z, (11) 


and P(7,¢) is a certain bilinear form in the variables 


1 = (Vardar IF? Vbs-- 8?) 
and C= (Za:Zas ve zt DZ psZbs one zi v), 
The differential expression /*(z) defined by the formula (11) is called the adjoint 
differential expression of /(y), and (10) is called Lagrange’s formula. 
If we apply the preceding considerations to the integral 


[, 1*(z). ydx, 


we obtain an expression of the form 


b b 
I, I*(z)ydx = Q(y,6) + I, z.Uy)dx. 
Hence the differential expression /(y) is adjoint to /*(z): 


I**(y) = Ky). 
In other words, the differential expressions Iy) and [*(y) are mutually adjoint. 
From the definition (11) of adjoint expressions we see at once that 


(7, + 2)* = if + Fy, (12) 
and, if Ais any number, 
(Al)* = XI*. (13) 
A differential expression /(y) is said to be self-adjoint! if /* = 1. It follows 
from (12) and (13) that: 
III (a). A sum of self-adjoint differential expressions is also a self-adjoint 
differential expression. 
(b). The product of a self-adjoint differential expression with a real number 
is also a self-adjoint differential expression. 
We now want to determine the general form of all self-adjoint differential 


expressions. 


IV. Any self-adjoint differential expression is a sum of differential expressions 


of the form 
LY) = (py) 5 lay aly) = HG py? PY + py), 


where p in each case is a real-valued function on [a,0}. 


1 In the English literature such an expression is often said to be formally self-adjoint. 
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Proof. We consider the integrals 
[. l,,(y)zdx and [ loy—a(y)Zax, with y,ze C™. 


By partial integration it is at once clear that /,,(y) and /,,_,(y) are self-adjoint 
differential expressions. Consequently a sum of expressions of the form 
l,,(y) and /,,_4()’) is likewise a self-adjoint differential expression. 
Conversely, let 
ly) = poy? + ry"? +e + Da 


be a self-adjoint differential expression; we shall prove that /(y) can be 
represented as the sum of expressions of the form /,,(y) and /,,_4(y). 
By definition, /(y) must agree with the adjoint expression 


*(y) = (—1)"(Boyy + (— 1) ay) + 4 Bay 
= (—1) "oy + [(—1) impo + (— I) Ay"? +... 
Hence, in particular, 
Po = (—1)"Do. (14) 


If nis even and = 2p, (14) implies pp = po, and so py takes only real values. 
If we subtract from /(y) the self-adjoint expression 


lou(Y) = (Do¥?Y” = poy” + (u — Ipoy*-P +... 


we obtain the self-adjoint expression (y) — /,,(y) of order n — 1. 

Next, if ” is odd and = 2u — 1, (14) implies pp = —p,, and therefore 
Po = ip, where p takes only real values. By subtracting from /(y) the self- 
adjoint expression 


lou 1) = ALCipy~ PY + (ip)%- Y 
= (ip)yP#"? + Spip'yP") +4... 
= poy” + tupoy"-? + ..., 
we obtain the self-adjoint expression /(y) — /,,_4(y) of order (n — 1). 

It follows from this argument that if self-adjoint expressions of the form 
In(y) and /5,-4(y) are subtracted successively from /(y), then finally a self- 
adjoint expression of zero order will be obtained. This must necessarily agree 
with /)(y) = py, where p denotes a real function. This proves assertion IV. 


If, in particular, /(y) is a self-adjoint expression with real coefficients, then it 
can have no addends of the form /,_ ,(y). Hence 


V. Any self-adjoint differential expression with real coefficients is necessarily 
of even order and has the form 


KY) = (Poy?) + (piy*- PHP 4 + (py + py 


where Po,Pi,-+-.P, are real-valued functions. 
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6. Adjoint Boundary Conditions. The Adjoint Operator 


Let Uj,...,U,, be linearly independent forms in the variables y,,y/,...,y"7 », 
VosVooeees Ves if m <2n, we shall supplement them with other forms 
Um+1s++-»U>, to obtain a linearly independent system of 2” forms U,,U,,...,Us,- 
Since these forms are linearly independent, the variables y,,y,,...,.y™~, 
VooVoo-¥"—-P can be expressed as linear combinations of the forms 
U1, U2,..., Von: 

We substitute these expressions in the bilinear form P(7,f) which occurred in 
Lagrange’s formula (see formula (10)). Then P(y,f) becomes a linear, homo- 
geneous form in the variables U,,U,,...,U.,, and its 2n coefficients are them- 
selves linear, homogeneous forms, which we denote by V,,,,V.,—1,---,V, 1n the 
variables z,,2/,...,20°~ ),z,,Z,,...,26"" ». Lagrange’s formula now takes the 
form 


[7 ).2dx = UyVan + UsVenna + 1 + Une +? yP@dx (15) 


The forms V,,V,...,Vo, are linearly independent. 

To prove this assertion, We first recall that the form P(7,¢) arose from the 
substitution of the limits x = a and x = bin the integrated terms obtained in 
the integration by parts, and so it has the form 


P(,6) = Pi(y,$) — Paln, 4), 


and P,(y,f) and P,(y,¢) involve only the values of the functions of y,y’, 
1. YP Y,z,2',...,2%-» at the points x = a and x = b respectively. 
Consequently the matrix of P(7,2) is of the form 


Po a 
; (16) 
0 A, 


where A, and A, denote the matrices of P,(y,¢) and P,(7,¢) respectively, and 0 


the null matrix. 
But from formula (9) we see at once that each of the matrices A, and A, has 


the form 


(—1)"~ "Po 
.(—1)"~ "po 0 
—py .0 0 
Po 0 ..0 0 


where p, is to be evaluated at x = aand x = b for A, and A, respectively. (The 
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elements above the diagonal in this matrix have not been written down, since 
they play no part in the argument.) 

The determinants of both these matrices .are non-zero, and from (16) it 
follows that the determinant of P(n,¢) is also non-zero. 

The transition from the variables y,,y/,...,.V"—-Y,y5,Vps-- Ye” P to the 
variables U,,U,,...,U., is a linear transformation with non-zero determinant. 
Hence the determinant of the matrix of P in the variables 


U4, Uo5...,U on ANG ZqpZjy- 0g Z YZ, Zp yee egZer P 
is also non-zero. But this matrix is also the coefficient matrix of the forms 
Von» Von—-++9 V4: 


These forms are therefore linearly independent. 
The boundary conditions 


V, — OV, = 0,...,Von—m = 0 (17) 


(and all boundary conditions equivalent to them) are said to be adjoint to the 
original boundary conditions. 


U, =0,U, = 0,...,U,, = 0. (18) 


Boundary conditions are self-adjoint if they are equivalent to their adjoint 
boundary conditions. 

Let L be the operator generated by the expression /(y) and the boundary 
conditions (18). The operator generated by /*(y) and the boundary conditions 
(17) will be denoted by L* and called the adjoint operator to L. 

It follows from formula (15) and the boundary conditions (17) and (18) that 
the equation 


[. Ly.2dx = |" y. D¥zdx (19) 


holds for the operators L and L*, for all y in the domain of definition of Z and 
all z in the domain of definition of L*. 
If we introduce the contracted notation 


b __ 
(9,2) = | v(x). 2@)dx, 
then (19) takes the form 
(Ly,z) = (y,L*2z). (20) 
It follows from the definition of adjoint operators that the operator L is 
adjoint to the operator L*: 
L** = LJ, 


An operator is self-adjoint if L* = L. 
It also follows from the definition of L* that: 
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An operator L is self-adjoint if and only if it is generated by a self-adjoint 
differential expression and self-adjoint boundary conditions. 
For a self-adjoint operator L the formula (20) becomes 


(Ly,z) = (y,Lz). (21) 


7. The Adjoint Boundary-Value Problem 


If L* is the operator adjoint to L, the homogeneous boundary-value problem 


L*z =0 (22) 
is said to be adjoint to the boundary-value problem 
Ly = 0. (23) 
In more detail, the boundary-value problem (22) can be written as: 
/*(z) = 0, (24) 
V (z) = 0, v = 1,2,...,2n — m. (25) 


where /*(z) is the differential expression adjoint to /(y), and (25) are the 
boundary conditions adjoint to those of L. 

We intend now to obtain a relation between the ranks of adjoint boundary- 
value problems. 

Let 2;,Z5,...,Z, be linearly independent solutions of the equation /*(z) = 0, 
and let r’ be the rank of the matrix 


V,(z,) Vi(z,) 
V2(2;) we Vo(2n) (26) 
Von - m(Z i) cee Von - m(Zn) 


Then the adjoint boundary-value problem (24),(25) has exactly (7 — r’) 
independent solutions. 

If, on the other hand, y is an arbitrary solution of the boundary-value 
problem (23), then the left-hand side of Lagrange’s formula (15) vanishes for 
z=z,. Also, U,(y) =... =U,(y) =0. Hence in this case Lagrange’s 
formula simplifies to 


Umn+1(¥) Ven —m(Zy) + oe + U2,Ay)Vi(2,) — 0. 


Putting v = 1,2,...,7, we obtain for the forms U,,4107),--->U2,() the system of 


equations 
Un -1V)Von- m(Z1) ee U2, y¥)Vi(z1) = 9, 


Un +1) Von~m(Ze) +» + Unaly) Vie) = 0, (27) 


Un+iV)Von—m(Zn) Fee oF U2) Vi(Z,) = 9. 
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This system has at least (n — r) independent solutions, viz. 


Un ti poses Von(V,)s J + 1,2,...,7 — fF, 


ifr is the rank of the boundary-value problem (23), and )4,...,y,-,18 a system of 
linearly independent solutions of (23). Therefore the rank of its matrix is not 
greater than 2n —m —(n—r) =n—m-r. Since the matrix of the 
system (27) corresponds—up to the arrangement of rows and columns—to the 
matrix (26), we have 


/ 


ro <cn—m+r. (28) 


We now recall that the two boundary-value problems are mutually adjoint, 
and therefore we can interchange their roles. In such an interchange, r and r’ 
are interchanged, and m goes over into (2n — m); hence r <n — (2n — m) 
+r’, 1.e. 

r<m—n-+r', (29) 


Comparing (28) and (29), we have 
r’=n—m-+r. (30) 


VI. The rank r of a boundary-value problem is related to the rank r° of the 
adjoint boundary-value problem by the relation 


r’=n—m+r. 
For the particular case m = n, formula (30) gives r’ = r. Hence: 


VII. If the number of independent boundary conditions is equal to the order of 
the differential expression, then the rank of the boundary-value problem is 
equal to the rank of the adjoint boundary-value problemi. 

In particular it follows from this that: 


VIII. If the number of boundary conditions is equal to the order of the differential 
expression, and if the associated homogeneous boundary-value problem has only 
the trivial solution, then the adjoint boundary-value problem has only the 
trivial solution. 


Additional information on the topics discussed in this section can be 
found in chapter 11 of Coddington and Levinson’s book [17]. 

We have been considering here problems where the boundary conditions 
are prescribed at the end-points of an interval. In mathematical physics and 
in applications, however, problems are encountered where the boundary con- 
ditions are prescribed at 7 arbitrary points of an interval [a,b]. Chapter 4 of 
Sansone’s book [97] may be consulted on these so-called n-point problems. 
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§2. Eigenvalues and Eigenfunctions of a Differential Operator 


1. Definition of Eigenvalues and Eigenfunctions 


A number Ais called an eigenvalue of an operator L if there exists in the domain 
of definition of the operator L a function » # 0 such that 


Ly = Aj. (1) 


The function y is called the eigenfunction, of the operator L, for the eigehvalue X. 
The operator L may be generated by the differential expression (y) and the 

boundary-conditions 
Ui’) = 0,...,U,0) = 0. (2) 


Since an eigenfunction y must belong to the domain of definition of the 
operator L, it must satisfy the conditions (2). Moreover, Ly = I(y), and 
therefore (1) is equivalent to 

Ky) = dy. (3) 
Hence: 
The eigenvalues of an operator L are those values of the parameter for which 
the homogeneous boundary-value problem 


l(y) = Ay, U,(y) = 0, vy = 1,2,...,7, (4) 


has non-trivial solutions; each of these non-trivial solutions is an eigenfunction 
belonging to A. 

A linear combination of eigenfunctions which belong to one and the same 
eigenvalue A is itself an eigenfunction belonging to 4. For, if Ly, = Ay, and 
Ly, = Ayo, then L(cyy, + Coe) = A(C1)1 + Coy.) for any constants c, and Cs. 

Since a homogeneous equation (3) can have, for a given A, not more than n 
linearly independent solutions, it follows that the aggregate of all eigenfunctions 
which belong to one and the same eigenvalue form a finite-dimensional space of 
dimension <n. The dimension of this space is simply the number of linearly 
independent solutions of the boundary-value problem (4) for the given 
eigenvalue A; this number is called the multiplicity of the eigenvalue. 

We shall specify conditions for the determination of eigenvalues. 

Let 

Yy(X,A),Vo(X,A), ++ Vn(%A) (5) 


denote a fundamental system of solutions of the differential equation (3) 
which satisfy the initial conditions 


Oforj # », 


| for j =» Jv = 1,2,...,n. (6) 


yy Paar) = 


From general theorems on the solutions of linear differential equations it 
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follows that, for each fixed value of x in [a,b], the functions (5) are integral, 
analytic functions of the parameter A. 

From the results of §1.4, the boundary-valte problem (4) has a non-trivial 
solution if and only if the rank r of the matrix 


Ui). UO) 
Vi 


UV) +++ Un (Yn) 

is less than n: 
and ifm <n,thenr <n; and in this case the boundary-value problem (4) has 
a non-trivial solution for any value of A. Hence, if m <n, any value of A is an 
eigenvalue. 

If i > n, the rank of the matrix Y will be less than n if and only if all its 
minors of order n vanish. But each of these minors is an integral, analytic 
function of A. Hence the only possible alternatives are: 


1°. All the nth-order minors of the matrix Y% vanish identically, and so, by the 
previous result, any value of X is an eigenvalue. 


2°. At least one nth-order minor of YW does not vanish identically. In this case, 
only the zeros of these minors can be eigenvalues, and further, a zero of one 
particular minor can be an eigenvalue only if it makes all the other, not 
identically zero, minors of order n of Y also vanish. 

Now, a non-vanishing integral function has at most denumerably many 
zeros (it need not have any at all), and moreover, these zeros have no finite 
point of accumulation. Hence: Jn case 2°, the operator L has at most denumer- 
ably many eigenvalues (it may have none at all), and these eigenvalues have no 
finite limit-point. 

Summarizing, then, we have the following alternatives. 

I. For any differential operator L, only the following two possibilities can occur: 

1°. Every number A is an eigenvalue for L; 


2°. The operator has at most denumerably many eigenvalues (in particular 
cases, none at all), and these eigenvalues can have no finite limit-point. 
The case when m = n is of particular interest, and in the sequel we shall 
consider this particular case, if nothing to the contrary is stated. 
We put, then, 
U,(1) Ui(y,) 


A(A) = , _ . (7) 
U,V) UY nd 


By the preceding discussion, A(A) is an integral, analytic function of 4, the 
so-called characteristic determinant of the operator L (or of the boundary- 
value problem Ly = 0); and the following theorems hold: 
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II. The eigenvalues of the operator L are the zeros of the function A(A). If A(A) 
vanishes identically, then any number 4 is an eigenvalue of the operator L. 

If, however, A(A) is not identically zero, the operator L has at most de- 
numerably many eigenvalues, and these eigenvalues can have no finite limit- 
point. 

If, in particular, the function A(A) has no zeros at all, then the operator L has 
no eigenvalues. 

An eigenvalue A may be a multiple zero of A(A). 


III. If Ay is a zero of the characteristic determinant A(A) with multiplicity v, then 
the multiplicity of the eigenvalue 4) cannot be greater than v. 
Proof. Let r be the rank of the matrix belonging to the determinant A(A,); then 
the multiplicity of the eigenvalue A, is equal to (7 — r). On the other hand, from 
a well-known rule for the differentiation of determinants, all the derivatives of 
the determinant A(A) up to the (n — r — 1)th order inclusive vanish for A = Ap. 
Since A, is a zero of multiplicity v, it follows that n —r—1 <v —1, and 
hence n —r <yv, 

If, in particular, vy = 1, then n —r <1; on the other hand, n —r 231, 
since A(A,) = 0. Hence n — r = 1, and we have proved the following: 


IV. If A) is a simple zero of the characteristic determinant A(A), then the 
multiplicity of the eigenvalue A, of the operator L is also unity. An eigenvalue 
is called simple if it is a simple zero of the characteristic determinant A(). 


2. Various Generalizations of Eigenvalue Problems 


In applications, more general eigenvalue problems occur from time to time. 
The most general formulation of the problem is as follows: 

The coefficients of a linear differential expression I(y), and perhaps also those 
of the boundary conditions U,(y) = 0,v = 1,2,...,n, may depend on a certain 
complex parameter X. The values of are to be determined for which the 
boundary-value problem 


Ky) =0, Uy)=0, v= 1,2,...,1 (8) 


has a non-trivial solution. 

These values of the parameter A are called the eigenvalues of the boundary- 
value problem (8), and the corresponding non-trivial solutions y are called its 
eigenfunctions. 

Using the same considerations as in the preceding subsection, we construct 
the characteristic determinant A(A) according to formula (7); the zeros of A(A) 
will again be eigenvalues. 

If the coefficients of the expression /(y) and of the boundary-conditions are 
analytic functions of the parameter A which are regular in a certain domain G 
of the A-plane, then the determinant A(A) will also be a regular, analytic function 
of Ain the domain Y. There are then the two following possibilities: 
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1°, A(A) = 0 in the domain G; then every number 4 € & is an eigenvalue. 


2°. A(A) ¥ 0; there then exist in the domain Y at most denumerably many 
eigenvalues and these have no point of accumulation within Y 


If, in particular, the coefficients of the expression /(y) and the forms U,(y) 
are integral, analytic functions of A, then A(A) 1s also an integral function; and 
hence in this case the theorems enunciated at the end of the last subsection 
remain valid. The most commonly occurring case is where ((y) = /,(y) + AL(y), 
and both /,(y) and /,{y) and the forms U,(y) are independent of A. The 
coefficients of the expression /(y) are then linear (and therefore integral) 
functions of the parameter A. So in this case, too, the theorems states in sub- 
section 2 hold. 

We may mention, further, the following important generalization of the 
eigenvalue problem defined above: 


L(y) + Al,(y) — 0, Uy) —= 0, y= 1,2,...,7, 


where /,,/, are self-adjoint differential expressions and U,(y) = 0(v = 1,2,...,n) 
are self-adjoint boundary-conditions relative to the differential expressions 
,(y) and /,(y). This problem has been analyzed in detail by E. Kamke (see 


[45]). 

3. Associated Functions 

We now turn our attention to the generalized eigenvalue problem 
Ky) =0, Uy) =0, v= 1,2,...,n 


for which the coefficients of the differential expression /(y) and of the forms 

U,(y) are integral analytic functions of the parameter A. Let Ay be an eigenvalue 

and (x) an eigenfunction belonging to A, for this problem. The functions 
P1(X),Po(X),.- «54 (%) 


are said to be associated with the eigenfunction g(x) if all the functions 
P(X), Po(X),..-59,(x) Satisfy the boundary conditions 


Yr 1 so 
pl OnP U.(Py- p) — Q, Po - Ff = 0,1,...,4; v= 1,2,...,72, (9) 


for A = Ao, and if for A = A, the following relations hold: 


Kp) =0 
Kp) +i2 it aA ° Kg) = 
(10) 
Kp.) + fF 9) ° Kee p+... $— sa oI “ (9) = =0 
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ak 

Here aya 

derivatives with respect to A of the corresponding coefficients of the differential 

expression /(%4). The number k is called the /ength of the system of associated 

functions. An eigenfunction ¢(x) is said to have multiplicity m if there is a 

system of functions associated with ¢(x) of length (#7 — 1) but no system of 
length m. 


i(ys) denotes the differential expression whose coefficients are the gth 


V. If Ay is a v-fold zero of the function A(A), then the multiplicity of any eigen- 
function is not greater than v. 

Proof. We assume the contrary; there is a system of functions associated with 
y(x) of length v: 


P1(X), PoC), vee 1Py(X). 


Let y,(x,A) be a solution of the equation /())) = 0, which, for x = a, satisfies the 
same initial conditions as the function 


Y(x,A) = p(x) + (A — Apdpr(x) + (A — Ao) Ppo(x) + «+. + CA — Ao)” py). 


Then the first (v + 1) terms in the expansion of the function y,(x,A) in powers 
of (A — Ay) agree with the corresponding addends in the representation of the 
function #(x,A). For, if we put A = A, in the relations 


Ky) = 0 
yar) =¥%ar, g =0,1,...(n — 1] 
it follows that the function y,(x,A,) satisfies the same equation and the same 
initial conditions for x = a as the function ¢(x) does; hence y,(x,A) = ¢(x). If 


we differentiate the relations (11) with respect to A, and then putA = Ag, we get, 
for A = Ao, 


(11) 


oy, o _. 
i() + 2 he) =o, 


(q) 
Wa _ gO@, — q =0,1,...(n — 1), 


OA 
ne) 
OA 


and comparison with (10) shows that | = @, (Xx). 
A 


Repeating the same argument, we arrive at 
1 67, (x | 
— = ¢,(x), = 1,2,...,v. 


Now we construct solutions yo(x,A),V3(X,A),-+-:Yn(%5A) of the equation (y) = 0 
such that the functions 


YX, A), Yo(,A), -° »Vn(XA) 
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are linearly independent for 4 = Ay, and with the help of these functions we 
form the determinant A()). 

Since, by hypothesis, the functions 9,¢,,...,9, satisfy the boundary con- 
ditions (9), the power series for U,(y,) in powers of (A — Ap) must begin with 
(A — A,)’*1. The first column of the determinant A(A) would therefore contain 
the common factor (A — A,)’t!, 1.e. the number Ay would be a zero of A(A) with 
at least (v + 1)-fold multiplicity, and this contradicts our hypothesis. Hence 
the theorem. 

Now let ¢, be an eigenfunction of maximum multiplicity, p, an eigenfunction 
of maximum multiplicity among the eigenfunctions linearly independent of 
P,, and, in general, y,; an eigenfunction of maximum multiplicity among the 
eigenfunctions linearly independent of ¢,92,...,p;-1- Clearly, the total number 
of such eigenfunctions is equal to the multiplicity p of the eigenvalue Ay. We 
denote the multiplicity of the eigenfunction y,; by m,, and the system of 
associated functions belonging to it by 


P j12Pj29°++9Pj,mj—-1° 

It follows immediately from the definition of the numbers m, that 

mh 2AM, 2... 2M, 
The system of functions so obtained 

PjPjrrPjarr--Pjmy-v J = 1,2,...5D, 
is called a canonical system of eigenfunctions and their associated functions. 
VI. The sum of the multiplicities my + m, + ... + m, is equal to the multiplicity 
of \) as a zero of A(A): 
1.e. m, +m, +... +m, =». 
Proof. Let 
by 41(X) Py 4 o(X),..-Bn(x) 


be any solutions of the equation /(y) = 0 for X = Ao, which together with the 
eigenfunctions 9,92,...,9, form a fundamental system of solutions of the 
equation ((y) =0 for A =A,. As above, we can construct the solutions 
YisYoo++sYn Of the equation ((y) = 0 whose expansions in powers of (A — Ao) 
have the form 


VA) = BAX) + PAA — Xo) Foe HE Pime1Q — Ag"? 
+ bx)A —Agy™ + ..., for i = 1,2,...,p, (12) 


yAx,r) = ${x) + ... for i = (p + 1),...,4, 


and these functions represent a fundamental system of solutions of the 
equation (y) = 0 for A = Ay. We form with them the determinant A(A). The 
first p columns of A(A) contain respectively the factors (A — Ag)™,(A — Ag)™,..., 
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(A — A,)"», and so the determinant A(A) contains the factor (A — A,/™tmt tmp 
Therefore 
A(A) — (A _ Xo) t mat Ft MpA (2), 

and A,(Ao) = det[U;,], i,7 = 1,2,...,.n, with U,;; = U,y;) for i = 1,...,n 
and j = (p + 1),(p + 2),...,”, 


1 o 1 @ 
but U;; — U(;) -+- iT ar U(P jm )-1) + 1 Ar? ULPj,m,-2) + _— 
4M uy ) fori = 1 dj=1 
my! evn Pj i= 1,...,4 and J = 1,...,p. 


To prove that v =m, + m,-+... + m,, we have only to show that 
A,(Ao) # 0. Let us assume the contrary, that A,(A,) = 0. Then one column of 
this determinant, the column with suffix j for example, is a linear combination 
of the columns with suffixes greater than /; let «,;4,,...,0, be the coefficients of 
this linear combination. 

First, let 7 > p. Then this implies that 


us — ab ~ 4 Pi4d re On, 
satisfies the boundary conditions U,() = 0, v = 1,...,n, for A = Ao. Since in 
this case all the ¢,,5;43,....%, are solutions of /(y) = 0, 4 would be an eigen- 
function and would have to be representable as a linear combination of the 
functions 9,,92,...,p,- But this contradicts the linear independence of the 


SYSTEM $1,Po5--- Ppp +i. oPn: 
If, on the other hand, j < p, then we put 


Pi(x,A) = (A — A)” ™y(x,A) fori <p, 


Pilx,d) — (A a Ao)™Iy (2,0) for i > DP; 
and Wx,A) = J (x,A) — O51; 41(%,A) — 1. = OP y(X,A). 
Then Ly) = 0. 


By differentiating this equation m, times with respect to A and then putting 
A = dj. it follows that the functions 


1 dy(x,A) 1 2?y(x,A) 1 amsy(x,A) __ 
— 9 goers — u(x) 
1! ocrdA 2! ada? mi oaA™ 


satisfy the equations (10) with k = m,. On the other hand, it also follows that 
the boundary conditions (9) are satisfied by the same functions, again with 
k = m,, and in particular the last equation 


m1 PU, (| a) =0 
op! ar? \(m; — p)! ade 


holds, again because of the vanishing of A,(A,) and the special construction of 
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y(x,A). These functions would therefore form a system of functions (of length 
m,) associated with the eigenfunction y(x,A9). But this is impossible, because 
y(x,A9) is a linear combination of the eigenfunctions 9;,9;+41,-..,pp, and there- 
fore, by the definition of these functions, it has a multiplicity of at most m,. 

Theorem VI means that a canonical system of eigenfunctions and their 
conjugate functions generate a subspace, the dimension of which 1s equal to the 
multiplicity of the zero of A(A) (for it can readily be seen that the eigenfunctions 
and associated functions of the canonical system are linearly independent), so 
that the associated functions complete the set of eigenfunctions to form a space 
of the necessary dimension. 

The theory of associated functions is analogous to the theory of elementary 
divisors of linear operators in a finite-dimensional space (cf., e.g. [28a]). The 
presentation of the theory of associated functions given here was originated 
by Keldysh [47]. 


4. Relation between the Eigenvalues and the Eigenfunctions of Adjoint Operators? 


THEOREM |] 
If X is an eigenvalue of multiplicity p of the operator L, then d is an eigenvalue 
of the adjoint operator L* and has the same multiplicity p. 
Proof. The operator L is generated by the differential expression ((y) and the 
boundary conditions 

U,(y) = 0, v = 1,2,...,n, 


while the operator L* is generated by the differential expression /*(y) and the 


conditions 
Vi(y) = 0, v = 1,2,...,n. 


We put h(y) = Uy) — Ay; 
then [¥(y) = T*(y) — Ay. 


If A is an eigenvalue of multiplicity p of the operator L, then the boundary- 
value problem 
hy =0, Uy) =0, v=1,2,...n 


has p linearly independent solutions. 
By the results of §1.7 (Theorem VII) it follows that the adjoint boundary- 
value problem 
i*(y) = 0, Vi(y) = 0, y = 1,2,...,n, 


also has p linearly independent solutions. This means that 4 is an eigenvalue of 
the operator L* and has multiplicity p. 
Two functions )(x) and z(x) are said to be orthogonal if (y,z) = 0.2 


' We remind the reader that only the case 2 = n is now being treated. 
2 We remind the reader of the notation introduced in §1: (y,z) = f; y(x).z(x) dx. 
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Eigenfunctions of the operators L and L* corresponding to the eigenvalues 2 
and w respectively are orthogonal if X # ji. 

Proof. Let y be an eigenfunction of the operator L for the eigenvalue A, and z 
an eigenfunction of the operator L* for the eigenvalue pu, so that 


Ly = Ay 
and L*z = pz. 
Hence (Ly,z) = (Ay,z) = A(y,z), 
and (y,L*z) = (y,nz) = 20,2). 


But from §1.6 (equation (20)), (Ly,z) = (y,L*z); so on subtracting the last 
two equations we obtain 


If A # ji this implies (y,z) = 0, i.e. the functions y and z are orthogonal. 


5. Eigenvalues and Eigenfunctions of a Self-adjoint Operator 


THEOREM 3 


All eigenvalues of a self-adjoint operator are real. 
Proof. For a self-adjoint operator L we have 


(Ly,z) — (y, Lz) ) 


(Ly,y) = (y,Ly). 


But (y,Ly) = (Ly,y); hence (Ly,y) is real. 

Now let A be an eigenvalue of the operator L, and y an eigenfunction belong- 
ing to A. The equation Ly = Ay implies that (Ly,y) = A(y,y). 
(Ly,y) 
(V9) 


in particular, 


must also be real. 


Since (y,y) > 0 and (Ly,y) is real, A = 
Theorems 2 and 3 yield the 


COROLLARY 


Eigenfunctions of a self-adjoint operator which belong to different eigenvalues 
are orthogonal. 

Eigenfunctions which belong to one and the same eigenvalue can also be 
chosen so as to be mutually orthogonal. To do this, we simply choose an 
orthogonal basis in the subspace of eigenfunctions which belong to one and the 
same eigenvalue. 

Altogether, we obtain in this way a system of pairwise orthogonal eigen- 
functions of the self-adjoint operator L. 
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6. Examples of Eigenvalue Problems 
L. —y"=dy; pO) =y(1)5> yO) = y'(). (13) 
The general solution of the differential equation is 

y =Acosx./A + Bsin x./d. 


By substitution into the boundary conditions, we see at once that the eigen- 


values are given by 
A, = (2n7)?, n = 0,1,2,3,... 


For n #0, there are, corresponding to the eigenvalue Aj, two linearly 
independent eigenfunctions 


cos (2n7x) and sin (27x); 


A, is therefore a double eigenvalue for n # 0. 
For n = 0, there is, corresponding to the eigenvalue Ay, only one eigen- 
function y = 1 (up to a constant factor); Ay is therefore a simple eigenvalue. 


2. Flexure of a Rod with one end fixed and one end free. 


We consider arod with one end, x = /, firmly fixed (Fig. 1) and its other end 
free, x = 0. At the free end a compressive force P acts along the rod. 

It is a familiar fact that for sufficiently small values of P the straight line 
shape of the rod is stable, but there is a critical value Py of the force P such that, 
for P > Po, the rod loses its straightness and bends. 

We shall investigate this bending by considering an equilibrium position 
differing but little from the straight-line position (Fig. 2). 

The equation for the bent axis of the rod y = y(x) is then such that 


Py = —Ely", (14) 


where 7 is the moment of inertia of the cross-section of the rod and E is the 
modulus of elasticity. The two sides of this equation represent the bending 
moment due to the applied force and the bending moment in the rod at 
distance x from the origin (see [111], for instance). 

We restrict ourselves to the simple case of a homogeneous rod of constant 
cross-section; then the stiffness E7 is constant. Using the substitution A = - 
we obtain the equation 


From Fig. 2 it is clear that the function (x) satisfies the conditions 


y0)=0, y()=0. (16) 
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Fig. 1. Fig. 2. 


So, by (15) and (16) we have formulated the eigenvalue problem. By sub- 
stituting the general solution 


y =AcosxJA+ Bsinx Ja 


for equation (15) into the boundary conditions (16) we find A =O, 
cos ,/A/ = 0, and so 


__ 2 
\= ("5 ‘| -  n=123... 


These eigenvalues are simple, and to them correspond the eigenfunctions 
2n — 1 
2/ 
The critical loads corresponding to the eigenvalues are 
2n — 7 2 
2/ 


Apart from a numerical factor, the eigenfunctions describe the corres- 


ponding equilibrium forms of the rod (Fig. 3). 
For arod with varying cross-section or for a non-homogeneous rod, E/ is a 


TX. 


y(x) = B sin 


P, = EI( 


function of x. Writing Fi = p(x), we obtain the following eigenvalue problem: 


—y"=Pp(x)y yO=90, y()=0. 


In the general case for an arbitrary function p(x), only approximate methods 
are possible for calculating the eigenvalues and eigenfunctions. 
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P, P, * P, 
v(x) yo(x) Ya) 


Fig. 3. 


3. Flexure of a Rod with one end fixed and one end supported. 


In this case there is a horizontal reaction H at the end x = 0 (Fig. 4). So the 
bending moment is M = Py — Hx = —Ely". 
By differentiating this equation twice, we get 
(Ely" 4 —_— —Py"; (17) 
and the function )(x) must satisfy the boundary conditions 


yO =0, yO=0, yp)=0, yO=9. (18) 
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The first and third of these boundary conditions are obvious; the second means 
that the bending moment vanishes at the point x = O(since the rod is supported 
there); the fourth holds because the end x = /is fixed and therefore the tangent 
there is parallel to the x-axis. 

So we have the eigenvalue problem (17), (18). Suppose now that EV is 
constant; we write P/E] = k?, and also, for simplicity, take / = 1. Then 
equation (17) can be written in the form 


yo — —k*y", (19) 
Its general integral, fork # 0, 1s 
y=A+Bx+Ccoskx + Dsinkx. (20) 


The boundary conditions (18) give directly that 4 = C =0 and that the 
eigenvalues are the roots of the equation 


tank =k. 
These are the abscissae of the points of intersection of the graphs (Fig. 5) 
y =tanxandy =x. 


The eigenvalues so obtained are simple, and the corresponding eigenfunctions 
are 
y =sinkx — xsink. 


— tan k 


paul 


Fig. 5. 
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For k = 0, the general integral of the differential equation Is 
y=A+t Bx + Cx? + Dx; (21) 


N 
but the boundary conditions (18) can be satisfied onlyif A = B= C= D=0. 
Consequently & = Ois not an eigenvalue of the problem. 


4. y= ky yO =V" =H) =y'() =0. 


For k = 0, the general integral (21) of the differential equation satisfies the 
given boundary conditions only if A=B=C=D=0. If k #0, the 
general integral (20) of the differential equations satisfies these boundary con- 
ditions again only if A = B = C = D = 0. Consequently this problem has 
no eigenvalues. 


5. yp” = —k*y’; y(0) = y'(0) = y"(0) = p(X) = O. 


As above, we find that k = 0 1s not an eigenvalue. If, for k # 0, we sub- 
stitute the general solution (20) of the differential equation into the boundary 


conditions, we obtain 
snk =k. 


This equation has no real roots other than zero. So the problem considered has 
only non-real eigenvalues. 


6. Natural Vibrations of a Stretched String. 


The equation of free, transverse vibrations of a string has the form 
(22) 
where u = u(x,t) is the displacement of the point of the string with abscissa x 


at time ¢, and a is a constant. 
Solutions of equation (22) which are of the form 


y = yx) sin (awt + ¢) (23) 


Fig. 6. 
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represent the natural vibrations of the string. Substitution in (22) shows that 
the function y(x) has to satisfy the differential equation 


yi + wy =0. (24) 


If we take the length of the string as unity, then the condition that the string is 
fixed at its ends (Fig. 6) can be written in the form 


y(0) = (1) = 0. (25) 


So we have an eigenvalue problem (24), (25). 
We substitute the general solution 


y=Acoswx + Bsin wx 


of equation (24) into the boundary conditions (25) and obtain A =O, 
sinw = Q, 
Hence 
w = nn, n = 1,2,3,.... 
So 
w? = (n7)?, n = 1,2,3.,... 


are the eigenvalues. They are simple, and to them belong the eigenfunctions 
y, (x) = B, sin nx. 


The eigenvalues and eigenfunctions have a simple interpretation in 
mechanics. Up to a constant factor, they are the squares of the natural fre- 
quencies of the string and the eigenfunctions yield the modal shapes for the 
natural vibrations. 


§3. Green’s Function for a Linear Differential Operator 


1. General Definition of the Inverse Operator 


An operator B is said to be inverse to an operator A if the domain of definition 
Q, of the operator B coincides with the range of values 2, of the operator A 
and if, for all x € D,, the relation 


B(Ax) = x 
holds. The inverse operator to A will be denoted by A~?. Hence, for all x ED,, 
A7~*(Ax) =x (1) 


I. If the operator A has an inverse operator A~, then the operator A is inverse 


to the operator A~?:1.e. 
(A-t)-1 = A. (2) 


For equation (1) shows that the range of values 2,-1 of the operator A~* 
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coincides with the domain of definition Z, of the operator A. Further, if we 
apply the operator A to both sides of (1) and put y = Ax, we get 


AA“ly = y’ 
for all ye &, = Y,-1. And this implies precisely that A is inverse to A7?. 


II. The operator A has an inverse if and only if the equation 
Ax =0 (3) 
has only the “trivial” solution x = 0. 
For, suppose that the operator A has an inverse operator A~}. If we apply 
A~' to both sides of (3), we obtain x = 0. Hence the condition (3) is necessary. 
If, on the other hand, we assume that the condition (3) holds, then the 
equation 
Ax =y, (4) 
for all ye &,, has a unique solution x €Y,; because, for each such y, at least 
one solution exists by virtue of the definition of the range #,. And the solution 
is unique; for if Ax’ = y were also true, then we should have A(x — x’) = 0. 
By hypothesis, this is possible only if x — x’ =0, 1e. if x = x’. To each 
element yEé&, we assign the solution x of equation (4). But we see im- 
mediately that this correspondence 1s a linear operator, and moreover it is the 
inverse to A. 


2. The Problem of Inverting a Differential Operator 


Let L be a differential operator, generated from the differential expression 


d®—ly 


dq” 
Wy) = polx) 55, + Pa) Go + + Pad 


and the boundary conditions 
U,(y) = 0, y = 1,2,...,n. 
We assume that the homogeneous boundary-value problem Ly = 0 has 


only the trivial solution y = 0; i.e. for an arbitrary fundamental system of 
solutions y,,¥e,.--,¥, of the equation y) = 0, the rank of the matrix 


Uy) Uy(v2) «Ui n) 
U(V1) Unlye) ... Usly,) 
UY ) U,()’2) cee U.Yn) 


is equal to n, and so the determinant 


det |Ui(y,)| # 0, i,j = 1,2,...,n. (5) 


§3.3 CONSTRUCTION OF GREEN’S FUNCTION 29 


By I, then has an inverse L~1, whose domain of definition coincides with the 
range of values of the operator L. We set ourselves the task of finding an 
explicit expression for L~?. As we shall see, L~+ turns out to be an integral 
operator with a continuous kernel; the kernel is called Green’s function for the 
operator L. We shall give below an explicit expression for the Green’s function. 


3. Construction of Green’s Function 
Green’s function for an operator L is to be understood to be a function G(x,€) 
satisfying the following conditions: 


1°. G(x,€) is continuous and has continuous derivatives with respect to x up to 
order (n — 2) inclusive for all values of x and € in the interval {a,)). 


2°. For any fixed value of € in the interval (a,b) the function G(x,€) has con- 
tinuous derivatives of orders (n — 1) and n with respect to x in each of the 
intervals [a,é) and (€,b]; the (n — 1)th derivative is discontinuous at x = € with 


a jump of a 
Pol) 
gn-1 ] 


3°. In each of the intervals [a,£) and (&,b], G(x,é), considered as a function of 
x, satisfies the equation (G) =0 and the boundary conditions U,(G) = 9, 
y = 1,2,...,n. 


THEOREM 1 


If the boundary-value problem Ly =0 has only the trivial solution, then the 
operator L has one and only one Green’s function. 

Proof. Let y1,Vo5---,¥, be linearly independent solutions of the equation /(y) = 0. 
Since the function G in the interval [a,é) also satisfies this equation, we have 


G(x,£) = ayy,(x) +... F AyVn(X) fora <x < &, 
where a,,d,...,4, are certain functions of €. Similarly, 
G(x,£) = byy,(x) +... + byp(X) forE<x <b. 


The requirement for continuity of the function G(x,é) and its first (n — 2) 
derivatives at x = & gives the equations 


lay) +. + aylO) — [ayn +. + OnvnlE)) = 9; 
lav +... + ay) — (iid) +. + bay] = 9, 


e 


[aye (€) +... +a, yO) — (brit PE) + + baw PO) = 0. 
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Moreover, the condition 


gn 


ona 1 


is equivalent to 


l 
[ayP-M(O) +... tay! MO] — (ay P® + + Say MO) = ~ pe 
0 
Putting 
c, = b, —a,, vy = 1,2,...,n, 
we obtain for the c, the system of equations 
aig) +... +eynlf) =9, 
cil) +... tent) =9, 
(6) 


cy WE) +... +o, y"- (6 = 0, 


ey) +... + ey Q) = + 
PoA§) 
the determinant of which is the Wronskian determinant for the fundamental 
system jj,...,y, for x = €, and therefore does not vanish. 
Hence the system always has a solution and the functions c, are uniquely 
determined by the system (6). To determine the functions a, and b, we bring in 
the boundary conditions; we write the form U,(y) as 


U,(y) —= UA ) + Uy ly )s (7) 


where U,,(y) is the sum of all the addends containing y,,y,,...,y"~, and 
U,,(y) is the sum of all the addends containing y,,y;,...,yi"~. 
Hence 


U,(G) = aU, Ai) + + anUya Vn) + AU Oy) +... + 6,00.) = 0. 
Replacing the coefficients a, by their expressions as a, = b, — Cy, we get 
by Uy (Va) + oe + On yoVn) + (Or — Cr) yar) +--+ On = Cn) Uva 'n) = 9. 
Hence, by (7), 
DU) we Hp UVVn) = C1U Va) +--+ nT y ald). (8) 


Puttingy = 1,2,...,n, we obtain from (8) a system of equations in the unknowns 
b,,b9,...5D,, for which, by (5), the determinant does not vanish. Consequently 
the system of equations has a unique solution 5,,5,,...,5,. But then the functions 
a, are also determined uniquely from a, = b, — c,. So we have proved the 
existence and the uniqueness of the Green’s function. 
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4. Inversion of a Differential Operator by means of the Green’s Function 


Suppose that the equation Ly = 0 has only the trivial solution y = 0, so that 
the inverse L~! (see subsection 1) and its Green’s function (defined as in 
subsection 3) exist. If L-4f = y, then 


Ly =f, (9) 


1.e. yis a solution of the equation 
(y) =f, (10) 
and also satisfies the boundary conditions 
U,(y) = 0, v = 1,2,...,n. (11) 


We shall show that this solution exists for any function f(x) which is con- 
tinuous in the interval [a,b], and that it can be determined by means of the 
Green’s function. More precisely, the following theorem is valid. 


THEOREM 2 


If the equation Ly = 0 has only the trivial solution, then, for any function f(x) 
which is continuous in the interval [a,b], there exists a solution of the equation 
Ly =f; this solution is expressed by the formula 


yx) = |, Gx, esde, (12) 


where G(x,é) denotes the Green’s function for the operator L. 
Proof. Let y(x) be the function defined by formula (12). We show that y(x) 
satisfies the conditions (10) and (11), so that it is a solution of the equation 
Ly =f. 

The function G(x,€) has continuous derivatives up to the (n — 2)-th order 
inclusive; so we may differentiate with respect to x under the integral sign in 
(12) (n — 2) times. Hence 


b ov 
yay = [SEO ede, vy = 12,040 — 2) (13) 


Hence the function y(x) and its derivatives y(x) up to the (” — 2)-th order 
inclusive are continuous in the interval [a,)]. 

o"-1G 
Ox" 1’ 


The function on the other hand, has a discontinuity at x = &. So in 


calculating y*-» and y™ we may not differentiate any more under the integral 
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sign without preliminary manipulation. For this reason we write formula (13) 
for vy = n — 21n the following form: 


yorm(a) = [2 MOOD) page + ip BOE) page 


Ox n—- 2 Ox" 2 


In each of the intervals (a,x) and (x,b) the integrand and its derivative with 
respect to x are continuous; we therefore differentiate with respect to x under 
the integral sign and with respect to the upper (or lower) limit x, and obtain: 


y"- V(x) as [ a *GO6E) Heyae 4+. eT ¥ f(x) 


Ox n-1 Ox" 2 
—“1G(x,€ On 2G(x, 
+ f POCO peag — [POO] pay a8 
Ox Ox é=x+0 
n— 2 
Since ie is continuous at x = &, the two integrated terms cancel out, and 


there remains 


yir-0(x) = [SOO oa + [? 


a OCS) Fede, (15) 


ox n-J] 


je yoy = [ 2 OGOD jeae (16) 


By differentiating the formula (15) again we find, as above: 


vx) = [PED poae 4 [A OOEO) fen 


> aG(x,€) o"-1G(x,§) 
OPE ANTE ae 
+ f PSS pode — [7S] peo. a7) 
From the condition 2° in the definition of Green’s function, 
ease) Z Lad a l 
£f=x-0 €=x+0 


ox ax Pox) 


Hence (17) can be written in the form 


yx) = [ PCC fg ae + fla) (18) 


The forms U,(y) contain merely the function values of }(x) and its derivatives 
up to the order (n — 1) inclusive at the points x = a and x = b; hence, by 
(12), (13), and (16), 


Uy) = |, U(@IOdé = 0, 


since, by definition of Green’s function, U,(G) = 0. So we have now shown 
that the function y(x) satisfies the boundary conditions (11). 
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We shall now show that it also satisfies the equation /(y) =f. 
By substituting the expressions (12), (13), (16) and (18) for the function f(x) 
and its derivatives into /(y), we find: 


Ky) = |? MG)fOdé + foo. 


But the integral in the last formula vanishes, because, by hypothesis, the 

function G(x,€), regarded as a function of x, satisfies the equation /(y) = O1in 

each of the intervals {a,£) and (€,5]. This completes the proof of Theorem 2. 
An operator A defined for all continuous functions f(x) by the equation 


Af(x) = |? kf Oak 


is called an integral operator with the kernel k(x, 6). 
Theorem 2 means that the operator L~' ts an integral operator with the 
kernel G(x,é). 


5. Green’s Function for the Adjoint Operator 


If L has an inverse, then, by virtue of Theorem VIII, §1.7, the adjoint operator 
L* also has an inverse. Consequently L* likewise has a Green’s function, which 
we shall denote by H(x,€). 

Let us find the relation between the Green’s functions G(x,€) and A(x,€) of 
Land L*. 

Put L-'¥f = g, and L*-1g = 4%, then f = Le, g = L*, and the equation 
(o,L*p) = (Lo,x) [see (20), $1.6] can be written in the form 

(L- 4,2) — (f,L*~*g). 

Hence 


b —. b f(b ———$____ 
[, { OOF Oeraxds = | ; { SHE x)eedxdé 
for any continuous functions f(x) and g(x). This gives us: 
G(x,€) = H(E,x). (19) 
Kernels G and H which are related to one another by the formula (19) are 


said to be adjoint to one another. 
So we have: 


Ill. The Green’s functions of two mutually adjoint differential operators are 
mutually adjoint kernels. 

In particular, if the operator L is self-adjoint, so that L* =L, then 
H(x,£) = G(x,€). The relation (19) then gives the condition 

G(x,£) = G(E,x). (20) 

A kernel which satisfies condition (20) is called Hermitian. Hence: 
{V. The Green’s function of a self-adjoint differential operator is a Hermitian 
kernel. 
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6. Boundary-value Problems involving a Parameter and their Reduction to an 
Integral Equation . 


In many questions we have to investigate boundary-value problems of the form 
Ky) = Ay +f), (21) 
U(y) = 0, v = J,2,...,n, (22) 


where A is a certain parameter, and f(x) denotes a prescribed, continuous 
function of x. Using the operator L generated by the differential expression 
I(y) and the boundary conditions (22), we can write this boundary-value 
problem in the form 


Ly =hy +f. (23) 
For the special case f = 0, this equation takes the form 
Ly =\y, (24) 


and is therefore equivalent to the homogeneous boundary-value problem 
Ky) =ay; UO) =90, v= 1,2,...,0. (25) 


The values of A for which equation (24) has a non-trivial solution are obviously 
the eigenvalues of the operator L (or of the boundary-value problem (25)), and 
each non-trivial solution itself is a related eigenfunction. 

Suppose that L has an inverse L~’; let G(x,&) be the associated Green’s 
function. By applying the operator L~! to both sides of equation (23), we 
obtain 


y=AL"y +8 (26) 
with g = L~}f. Equation (26) implies 
(x) =A] Gx, yOdé + g(a. (27) 


Hence: 


V. If the differential expression I(y), subject to the boundary conditions 
U,(y) = 0,v = 1,2,...,2, has the Green’s function G(x,€), then the boundary- 
value problem 


Ky) = Ay + f(x), UO) =0 
is equivalent to the integral equation 


v(x) =A] GOO yOdE + g(x), 


where 


(x) = |? G2 fOdé. 
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In particular, the homogeneous boundary-value problem 


Ky) — Ay, Uy) — 0, y— 1,2,...,7 


is equivalent to the homogeneous integral equation 


Wx) =A J” G(x, yOdé. (28) 


Since G(x,£) is a continuous kernel, the Fredholm theory (see, for example, 
[86a]) can be applied to the related integral equation. Accordingly, the 
homogeneous integral equation (28) has only denumerably many eigenvalues 
Ay,Ag,Ag,-.., Which have no finite point of accumulation. For all values of A 
different from the eigenvalues, the inhomogeneous equation (27) with an 
arbitrary continuous function g(x) on the right-hand side has a solution 
expressed by the formula 


yx) =) [ PoENgOdé, (29) 


where I'(x,€,A) denotes the resolvent of the kernel G(x,€) and is a meromorphic 
function of A, for arbitrary fixed values of x and € in the interval [a,b], which can 
have poles only at the eigenvalues of the homogeneous integral equation. From 
the integral equations we go over to the corresponding boundary-value prob- 
lems to obtain the following results: 

VI. If the differential equation Ky) =0, under the boundary conditions 
U,(y) = 0, v = 1,2,...,n, has only the trivial solution, then: 

(a) the homogeneous boundary-value problem 


iy) = Ay: Uy) = 9, vy = 1,2,...,n, 


has at most denumerably many eigenvalues ),,Ao,A3,... with no finite point of 
accumulation. 
(b) the inhomogeneous boundary-value problem 


Ky) =Ay +f: UO) =0, v =1,2,...0, 
has a solution for all values of 2 different from the eigenvalues for the corres- 


ponding homogeneous problem and for any continuous function f(x). 


We had already obtained proposition VJa otherwise, without using integral 
equations, in §2.1; and we could have obtained VIb by considerations similar 
to those there used. 


7. Green’s Function for the Operator L — Al 


As before, let L be an operator generated by the expression /(y) and the con- 
ditions U,(y) = 0, v = 1,2,...,7. We want to find an expression for the Green’s 
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function of the operator L — A1; in other words, we want to form the inverse 


of the operator L — Al. 
We denote by y, = y,(x,A), v = 1,2,...,7,,a system of solutions of the 
equation /(y) = Ay which satisfy the initial conditions 


Oforj # », 


Ye MOM) = Ni for j =v 


We apply the method of variation of the parameters to the equation 


lo) — =f, 
and obtain 
vo) =¥ cro) + ['4O(Y PO Jae (30 a) 
and also 
vo) =¥ creo — f s0(S F)ae. (30 b) 


Here Cj,...,C,,Cj,...,C, are constants and W 1s the Wronskian determinant of 
the functions y,,V2,...,Vn: 


yo- 1) ye- 1) ee yo 1) 
yo 2) yo 2) _ ya- 2) 
W = 
| V1 y2 w+ Vp 


and W,,W,,...,W, are the cofactors of the elements in the first row of W. 
Adding (30 a) and (30 b), and dividing by 2, we get 


y(x) = Cy) + f° elxseflede, 31) 
where C\,...,C,, are constants, and the function 
y,(x) yx) s+ Ya(X) 
(n— 2) (n— 2) (n— 2) 
_ po | wr) yer) Yn *(E) 
g(x,6) = + We a (32) 
| Yi) yl€) se WalS) | 


the positive sign being taken if x > &, and the negative signif x < €. 
We next demand that the function y tn (31) shall satisfy the boundary con- 
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ditions U,(y) = 0, v = 1,2,...,n;1.e. that it shall be a solution of the boundary- 
value problem 


Ky) = ay +f, U,(y) = 0, y = 1,2,...,n. 
This requires that 


LUO + [AOU ae =0 


By solving these equations for C,; and substitution in (31), we get 


w(x) = | GE (Ode (33) 
with 
G(x,6,) = "—" ex, £,a) (34) 
9335 A(A) 93 ? 
where 
Oi(y1)  Uy(y2) «.. ~Uy(y,) 
U, 1 Us 2 cee U, n 
AQ) = (V1) U2(y2) Or) (35) 
Uy 1) Uy 2) Uy ”) i 
and 
V(X) yx) vee Val) 2(x,€) 
Ui) U2) «.. UiQ,) U,(g) 
A(x,€,A) =| Us(y1) Us(ve) «.. Ualy,) U2(g) (36) 
Us») U0») Url») U gla) | 


If A is not an eigenvalue of the operator L, then A(A) # 0, and the formulae (33) 
and (34) are meaningful. Formula (33) shows that G(x,é,A) is the Green’s 
function of the operator L — Al. Hence: 

VII. The Green’s function G(x,€,A) of the operator L — X1 is determined by the 
formulae (32), (34), (35) and (36). 


8. The Analytic Nature of the Green’s Function of the Operator L — Al 
The functions A(A) and A(x,é,A), as may be seen immediately from (35) and 
(36), are integral, analytic functions of the parameter A. So (34) implies: 


VIII. Green’s function G(x,€,A) for the operator L — Al is a meromorphic 
function of the parameter X, and its poles can only be the eigenvalues of the 


operator. 
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Let A, be a simple zero of the function A(A). Then A, can also only be a simple 
pole of the function G(x,€,A), so that | 


Res 2 


G(x,€,A) = + Gs f,A), (37) 


where G,(x,€,A) is regular in a neighbourhood of the point A). 

If, then, A, is not a singular point of the function G(x,€,A), we have to put 
R(x,€) = 0. 

By a well-known formula in the theory of residues, we have 


h H(x,€,A) 

R(x,£) = (—1) Ay)” 

In the expansion of the determinant H(x,€,A) by elements of the first row, the 
coefficient of g(x,€) is A(A,) and = 0; hence the function R(x,€) and its first n 
derivatives with respect to x are continuous in the square a <x,€ <b. 
Finally, it follows from (36) that R(x,€) is a linear combination of the functions 
y(x),.-.,¥, (00); consequently, for a fixed &, R(x, €) satisfies the equation 


I(R) — A(R) = 0. 


We also derive from formula (36) that H(x,€,A), and therefore R(x,&) also, 
satisfy the boundary conditions. 

Therefore, for every fixed €, R(x,£) is an eigenfunction of the operator L for 
the eigenvalue Ap. Since, however, Ay 1s a simple zero of the function A(), there 
is only one eigenfunction yo(x) corresponding to it, up to a factor independent 
of x, for the operator L. We must therefore have 


R(x,£) = a(€)Vo(). (38) 


We now put G*(x,€,A) = G(é,x,A); then G*(x,€,A) is the Green’s function for 
the operator L* — Al, provided that the coefficients of the expression ((y) 
satisfy the differentiability requirements needed for /* to exist. 

From equation (37), on the other hand, we infer that 


G*(x,£,.) = REx ~ 


Tt G,(E,x,). 


So, for a fixed f R(é,x) isan eigenfunction of the operator L*, corresponding to 
the eigenvalue Ay. If we denote one of these functions by z,(x), we then have 


R(g,x) = b(€)z,(x). 


Hence 


R(x,€) = b(x).Z9(€). 
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Comparing this formula with (38), we find 


R(x,E) = cyo(x)z0(). (39) 


It remains to determine the constant c. 
From (37) and (39), we deduce that 


(A — Ag) 7 G(x,EA)yolEde = cve(x) |” vol QZK@ae 
+ (X= Ag) [? Gale, A)ya(E)de. (40) 


The integrand in the last term is a regular function of the parameter A in a 
neighbourhood about the point Ay; hence we find from (40) that 


lim (A — A) [” Gee Ayo(Cae = eval) [7 volO70ae. (A) 


A—>Ay 


On the other hand, 
(L — Al)yo = (Ao — A)Yo, 
whence 
l 
—1)-1y, = —__ ». 
(L 1) 1H yy — 120 
Since G(x,é,A) is the Green’s function for the operator L — Al, the last equation 
can be written in the form 


b 
[ Ge.trrlde = 5 00. 
a 0 
By substitution in (41) we obtain 
—yolx) = evalx) [7 volE)zl@ ae, 


so that c is given by 
l 


[, vol)z0l@ae 


and therefore 
R(x,£) = — role) 
I, YolE)Zo(S)aE 
It follows from (37) that: 
IX. For every simple zero Xy of the function A(A), 


Glx,é)) = — 20070) Gx, ,0), (42) 
(X — do) J vol)z0t8)a¥ 


where G,(x,é,A) is regular in a neighbourhood of the point do. 
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The formula becomes particularly simple if the functions y.(x) and Zo(x) are 
normalized so that 


[7 vo()z0@ae = 1. 
It then becomes 


— "ee 


G(x,E,A) = + Gy(x,€,A). (43) 


9. The Case of a Multiple Pole of the Green’s Function 


The above results may be generalized to the case of a multiple zero A, of the 
function A(A). In general, A, is then also a multiple pole of the Green’s function 
G(x,é,A). We shall again assume that G(x,€,A) is the Green’s function for the 
boundary-value problem 


ly) =f, Uy) = 90, v = 1,2,...,n, (44) 


where the coefficients of the differential expression /(y) and of the forms U,(y) 
are integral, analytic functions of the parameter A; the boundary-value problem 
(L — Al)y =f is obviously a special case of this boundary-value problem. We 
observe that all the results of sections 7 and 8 carry over completely to the 
general case, y,(X,A),Vo(X,A),---,),(%,A) now being understood to be a funda- 
mental system of solutions of the equation /(y) = 0, satisfying the same initial 
conditions as before, viz. 

l fork =», 


YEP) = lO tor k £ v. 


Let, then, Aj now be a pole of the »-th order for the Green’s function G(x,é,) 
for the boundary-value problem (44); then 


R(x, £) 4 Ry(x,£) 


R,- vals 2) 
OX — ry Oo —aAgesk , pf EE + Gy (x,€,4), (45) 


GUx,8,A) = 
G,(x,€,A) being a function regular in a neighbourhood of the point Ao. 
For a fixed é, the function G(x,€,A) satisfies the equation ((G) = Oin each of 
the intervals [a,€) and (€,5]. By substituting for the function G its expression as 
in (45), we obtain the relations, for A = Ag: 


IR) = 0, 
KR) +4 2 8 =0, 

(46) 
MRy a) + 5 <1(R,- yt. +t. 2 KR =0. 


(v — 7 or? 
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Moreover, it follows from the formulae (34) and (36) that, with a fixed é and 


for A = Ao, the functions R(x,é),R,(x,§),...,R,-1(x,6) satisfy the boundary 
conditions 


H 1 oe 
» 


my pl on UR, - p) = Q, Ro = k, p= 0,1,...,0v — 1), q= 1,2,...,77. (47) 


By making use of the relations (46) and (47) as well as the definition of the 
eigenfunctions and the associated functions ($2.3), we can satisfy ourselves of 
the truth of the following theorem, which we merely state without proof (see 
[45a]): 

X. The principal part of the Green’s function in a neighbourhood of the zero 
\9 of the function A(A) has the form 


> (ZA), ZO AG) + ZnO iO) 


Ala — Aye Q—Ayycr 
4 HOY impr) F Zi OY jmp 20) FF Zim TOV LCD (4g) 
X—A 
where 


V joV jryoe+oY jmy- ys J — 1,2,...,p 


is an arbitrary, canonical system of eigenfunctions and associated functions of 
the boundary-value problem concerned, while 


2 joZjrr-+ +92 jmp JI = 1,2,...,p 


is an appropriately normalized, canonical system of eigenfunctions and 
associated functions for the adjoint boundary-value problem. 


XN 
CHAPTER II 


ASYMPTOTIC BEHAVIOUR OF THE 
EIGENVALUES AND EIGENFUNCTIONS. 


EXPANSION OF A PRESCRIBED FUNCTION 
IN TERMS OF EIGENFUNCTIONS OF A 
DIFFERENTIAL OPERATOR 


$4. The Asymptotic Behaviour of the Eigenvalues and Eigenfunctions for 
Large Values of |A| 


1. Statement of the Problem 


For large values of |A|, approximation formulae and, indeed, asymptotic 
formulae can be given for the eigenvalues and eigenfunctions of a differential 
operator. Such formulae are not only of interest in themselves, but they also 
find application at a decisive point in the proof of certain theorems in the theory 
of differential operators, particularly in expansion theorems. 

It turns out that the behaviour of the eigenvalues A, and the eigenfunctions 
for an arbitrary differential operator as |A,| increases can, to a first approxima- 
tion, be characterized by the eigenvalues and eigenfunctions of the operator 
which is generated by the same boundary conditions but the simplest possible 


differential expression of the nth order, (y) = —. 


We shall first investigate the asymptotic behaviour of solutions of (y) = Ay 
for large |A|, where /(y) is an arbitrary differential expression. Then by sub- 
stituting the asymptotic approximation expressions so obtained into the 
equation A = 0 (cf. §2.1), we shall obtain the corresponding asymptotic 
estimates for the eigenvalues. 

First we put A = —p”; the equation /(y) = Ay then takes the form 


Ky) + p’y = 0, (1) 
or, in more detail, 
d"y d"—*y , 


[For simplicity, we have at first taken po(x) = 1; see §4.11(b) for the case 
42 


§4.1 THE PROBLEM 43 


Pox) # 1]. Without loss of generality, we may assume that p,(x) = 0. For if 
P;(x) # 0, then by using the substitution 


_ 
y= F exp(— | p(x), 
we find that equation (2) takes the form 


2 
2 + ple) ee yy . + DPi(x)y + p"y = 0, 


where p,(x),...,p,(x) are likewise continuous with respect to x in the interval 
[a,b], and p has not changed. 

Also without loss of generality we can replace the interval [a,b] by [0,1], for 
we can return immediately to the general case by making the substitution 


x =a+(b—ab)t. 


2. The Domains S and T 


We divide the complex p-plane into 2n sectors S,, k = 0,1,...,22n — 1), 
defined by 


(3) 


(Fig. 7). 

As we shall see later, the asymptotic formulae for the solution y of equa- 
tion (2) depend essentially on which sector S, the point p lies in. 

We denote by 


W1,Wo, ooe Wp 


the different nth roots of —1 arranged in an order in each case to suit later 
requirements. 
Later we shall need to use the following essential properties of the sectors S,. 


I. For each of the sectors S,, the numbers 01,wW9,...,0, can be ordered in such a 
way that, for all p € S,, the inequalities 
Rp) <R(pws) <... <A(px,) (4) 


hold, where A(z) means the real part of z. 
Proof. It is sufficient to prove the assertion for two adjacent sectors S,. For, 


multiplication by any particular nth root of unity, say 
e= e@intin 


transforms each sector S, into the sector S,, obtained from it by a rotation 
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S Wii 


; 
7 Z Vo Mi 


SY” 


“jy 


Fig. 7. Fig. 8. 


through the angle =, 0 <k <2n — 1,and we have k’ = k (mod2) ; and the 


numbers 


e—1w,,e7 1wo,...,€ tw, 


differ from the numbers w,,w»,...,w, Only in their order. 
The inequalities (4) can now be written in the form 


RB(pe.e~1w) <A(pe.e~ two) <... <A(pe.e7 4w,,), 


and by what has just been said, they provide the corresponding inequalities 
for the sector S,,. Since /is arbitrary, they give the corresponding inequalities 
for all S, with k’ — k = 0 (mod 2), provided only that they hold for the 
sector S,. Hence it suffices to prove the assertion for two adjacent sectors, say 
So,-1 and So. 

Suppose is odd: n = 24 — 1. Let, then, the order of the {w,} be chosen so 
that 


21 27 An 
arg w, = 7, arg w, = 7 + —, argu, =a ——, argzw, =r7+—, 
n n n 


dn 7 
arg We — Tt neo”? arg Ww, = n 


If we multiply any number p € S,,_, in turn by @4,w,...,w,, then the points 
Py>P20+++xP, SO Obtained will be situated as shown in Fig. 8. For the domain 
Son—- 1 therefore, the inequalities (4) are geometrically obvious. 

Similarly, we can show the validity of the assertion I for the domain S, and 
also for the case when nis even. 

Later we shall find it expedient to consider more general domains, namely 
those which are obtained from tlie sectors S, by a translation p > p — c, 
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where c is a fixed complex number. These new sectors with their vertices at the 
point p = —c will correspondingly be denoted by 7,, k = 0,1,2,...,2n — 1. 

Taking account of the way in which the T, are produced from the S, by 
translation, we see that, for p € T,, the inequalities 


A((p + c)uy) <A((p + c)we) <... <A((p + c)w,) (5) 


hold, for a suitable ordering of the numbers w,,w»,...,w,. 

In the sequel we shall let p vary in a fixed domain 7, and so we shall write 
simply S and T instead of S, and T,. The order of the numbers w,,ws,...,w, 
will be such that for p € T the inequalities (5) are valid. 


3. Reduction of the Equation /(y) + e”y = 0 to an Integro-differential Equation 


The derivation of the subsequent asymptotic approximation formulae depends 
on the fact that the equation 


Ky) + py = 0 (6) 


can be reduced to a certain equivalent integro-differential equation. To this 
end we put 


ny) = —pey""? — ... —Pr); (7) 
and we can then write equation (6) in the form 
y + p"y = m(y). (8) 


The homogeneous, linear, differential equation y” + p"y = Ohas, for p # 0, 


the fundamental system 
EPr® EPO ePOn™, 


If we regard (8) as an associated inhomogeneous equation with the inhomo- 
geneous term my), then the method of variation of the parameters gives the 
general solution y in the form 

x pa,(x—¢) POr(x— €) 
y = cyeP?* +. + 6,008 | nT aa my)dé, (9) 
where m,(y) denotes the value of m(y) for x = E, 

The equations (8) and (9) are equivalent in the sense that any solution of (9) 
with arbitrary constants ¢,,...,¢, is also a solution of (8), while conversely for 
each solution of (8) there exist constants c,,...,c, with which (9) is satisfied. In 
this, the constants c,,...,c, may depend on p. We keep k fixed at one of the 
values k = 1,2,...,n, and put 


for j = 1,2,...,k 
a ak ee , . (10 
c; =e; + y | oe mA(y)d€ for j = (k + 1),...5n (10) 


0 
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Equation (9) then takes the form 


_ | Ky o€,p)ine(y)deé 
0 . 


np” 


y = cye?®™ +. + ce? Or® + 


oe " Kylx,£,pyndydé, (11) 
np"—* Jy 
where 


k n 
K,(x,€, p)= » wy err 9), K,(x,€,p) — yy wer rae 9). (12) 
a=] k 


awxk+l 


Equation (11) is the integro-differential equation referred to. 


4. A Lemma concerning a System of Integral Equations 


LEMMA | 
If the system of integral equations 


V0) =f) + YP Auless Ode T= Aout 13) 


satisfies the following hypotheses: 
1. the functions f (x) are continuous in the internal [a,6]; 
2. for any fixed €, the functions A;(x,€,A) are continuous ina <x < & and 
E<x <b; 


3. for any fixed x and § (a <x,€ <b), the A;,(x,£,A) are regular analytic 
functions of the parameter 4, defined in a region of the complex d-plane 
extending to infinity; 

4. there exist constants R and C such that 

C 
|A ;(x,&,A)| < 1A 


for |A| > Randa <x,§ <b: 


then, for a sufficiently large Ro, there exists, for |A| > Ry, one and only one 
solution 


yilx) = y(x,A), i= I,...,7, 
such that the functions y(x,) are regular analytic functions of d for |A| > R, 
and such that 


] 
yi(x,A) = fx) + O (3) for A> o.} (14) 
1 The expression O(1/A") denotes, here and elsewhere, a function of the form f(X,r)/a*, 


where |f(x,A)| for a <x <b and for sufficiently large |A| always remains less than a 
constant. 
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Proof. If the system (13) has a solution y,(x),i = 1,2,...,r, then using the 
method of successive approximations we obtain the equations 


yx) =f) + YY” Alesha (Oe 


$Y FF Au Cb) ja (Erboy;, (Ede des 


Jisds =1 


A) + |” Aloo MG (OME + 


+ y { vee [ Ajj, (X,EpA)..-4; 5 (En-aF Sj, (Eder... 


Ji poeoos Jn=1 


+ y [. oe [ Ai (%EvA)---Aj, jn a (EnsEn tj, 4, (Eng dbs. db, 43. 


Jivserdngi=l 


We put B = max |y{x)|, a<x <b, i=1,2,...,r. 
Then, for |A| > R, the terms in the last row of this equation chain add up to 
an amount not greater than 


n+1 
jaye [((b — a)Cr}] B, 


and therefore tend to zero as n —> 00, provided |A| > Ro, where 
Ry = max {R,(b — a)Cr}. 


If this condition is satisfied, the solution is expressible in the form of an 
infinite series, viz. 


yi(x) = yx.) 
= fx) + YP? Aylok Nf (Cae 


+ y [- [. Ai (EAA; j,(Er$20F;,(F2)dbsdbo +... - 


Jvjg=1 


Conversely, it can easily be seen that this series converges uniformly for every 
case with |A| > Ry,a <x <4, and represents a solution of the system (13). 
Hence follow all the assertions of the lemma, and in particular the formula (14) 
and the regularity of the functions y,(x,A) for |A] > Ro. 
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5, Asymptotic Estimates for Solutions of the Equation /(y) + p”®y = 0 


First we need estimates for the kernels K, and K, in (11) and for their 
derivatives. 


LEMMA 2 


There exists a constant C such that for all p € T the following inequalities hold: 


 KiQ8.d) <Clpl"klerm-9| for 0 <é <x <1, usa) 


5 Kulséd) < Clon — Wler™*-9| for0 <x <e<1, (15) 


vy = 0,1,2,3.,.... 
Proof. We choose a constant C such that, for all j,k = 1,...,n and all x,é in the 


interval [a,b], 
| eCCos OW 9) <C: (16) 


this is possible, since the left-hand side of (16) is a continuous function of the 
variables x and &. It follows from the inequalities (5) that, for p € T and for 
a <k, 

K(pw,) < A(pw, + (p + cwy — ,)), 
which, forO <é <x <1, implies 


| ep Palx—8)| < |elerat + NOK Oa) I(e— 8) < Cer on) 


Therefore (see (12)) 


d” < - 
Keb] =| Pasremeno] < ck |operme“O 
a= 


xX 


which is the inequality (15a). 
Inequality (15 b) may be proved similarly. 


THEOREM | 


If the functions po,...,P, are continuous’ in the interval [0,1], then the equation 
y™ + pay?) +... + Pa-ry’ + Pay + p"y = 0 


has, for each region T of the complex p-plane, n linearly independent solutions 
Y1y---Yn» Which are regular for p€T and for sufficiently large |p|, and which, 


1 In fact, as will be seen from the proof given, the conclusion of the theorem remains true 
if po,...,Pn are arbitrary functions summable in the interval [a,6]. 
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with their derivatives, can be expressed in the form 


net +0(4)], 
P 


dy, 4X l 
dx . +o(°)], 


d"~*y, — .n=-1,p0,.x n—1 l 
aenet eet! + O(°)]. 


(17) 


Proof. Suppose that the equation (6) has a solution y, for which 


c, =Ofory #¢kandc, = 
Then 


O_x ] * ] 1 
y= erm a | KalsbedmeOddé — rs | KalesBiplmeCndae. (18 a) 


The results of differentiating this equation (m — 1) times, together with the 
original equation, give the system of equations 


v O 
Tot = prayers + —T5 [ SRE) my de 
dx Mp ° 


1 av 
od | TRIKE P) my (y, dé, v = 041,2,-.,(2 — 1); (18) 


np’~* ox” 


for, the result of differentiating with respect to the upper and lower limit x is 
equal to zero when the sum of the two expressions is formed; this can be seen 
by differentiating the integral in (9) with respect to the upper limit, and using 
the known result w}j +... +o, =Oforv = 1,2,.... 
If in (18) we put 
ayy 
dx” 


we obtain the system of equations for the functions Z,, = Z,,(x,p): 


= prerz. (19) 


LP pour ogev Ki 8 
zursp) = 0k — Lf” erponten Opry SEE |p 2, 5-60) 


Ox 


41 © PalS)Enn- Ep) fovea 1 p82, E0)| dé 


1 fi xe) ~y O K(x, §, 
+ eres Jo a Penn sl6) 
np x 


+1 * Pala Ep) on + 1 pézy a) de (20) 
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If we now put 
Qo” 
_i g-renten Op-v- 049 ORCS PA€) for € < x, 


K kvalX ,p) — 


an ae d ASE) 9.6 for E > x: 


k = 1,2,...,n3 v = 0,1,2,...,(n — 1); a = 2,3,...,n; 


then the system (20) can be reduced to the form 
1a 
zl) = OL +5 | Kinser Geddb. 21) 


For a fixed k, the n equations obtained from (21) by puttingvy = 0,1,...,@ — 1) 
form a system of integral equations for the functions Z,,(x.p), v = 0,1,...,7. It 
follows from Lemma 2 that all the functions K,,,(x,€,p) are bounded for 
0 <x,é <1,peT, and sufficiently large |p|. Hence all the assumptions in 
Lemma 1 of section 4 are satisfied. By this lemma, the system (21) has one and 
only one solution z,, = Z,,(x,p), Which is analytic in p and which satisfies 


Z4y(X,p) — Wy + o(2). (22) 

From this and from (19), the relations (17) follow immediately, and hence 
also the linear independence of the functions y,(x,A), k = 1,2,...,n. 

It remains to prove that a solution y,(x,A) of equation (6) exists, satisfying 
equation (18 a); to do this, it is only necessary to show that for arbitrarily 
chosen constants c, (i.e. for numbers independent of p) there is a solution of 
equation (6) which for these c, satisfies equation (11). 

The equations (10) represent a linear transformation of the c, into the c; (we 
remark that y and therefore also the expression m,(y) on the right-hand side of 
(10) depend linearly on ¢,,...,c,). Clearly it is sufficient to show that the deter- 
minant of the transformation (10) is different from zero for sufficiently large 
|p|, <7; we can then solve equations (10) with arbitrarily prescribed c; for c,. 
The solution y of equation (6), or, what is the same thing, of equation (9) for 
these values of c, will then be the required solution. 

If now the determinant of the transformation (10) vanishes for some p € T, 
then the equations (10) with cj; =c,; =... =c, =0 for this p have a non- 
trivial solution c;. The corresponding function y will consequently be a 
nontrivial solution of the equation 


l 
np" 


- 1 [. Ki(x,€,p) me ydé —_ 


np" 


y= =| Kiséipm()dé, (23) 


which is obtained from (11) on putting cy = cg =... =c/' = 0. 
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We shall show that, for sufficiently large |p|, this is impossible. 
By differentiating (23) (n — 1) times and using the substitutions 


Zy >= 


~ prerP™*)-1-y = 0,1,...,(n — 1), (24) 
we obtain for the functions z, the system of equations 
Z,(x,p) = (25) 


l x 
_ 7m I. p7 Ven Pants Ke As TRE Pe, n—2y (£)z, ~2(&,p) +... +p,(6)Z0(£,p)}dé + 


at [. p~ "em POKte— k) o° As oO Kx, 6,p) {p"-? 'Do(E)Z,—-o( Ep) +. + Dal 2)Zo(Esp)}dé. 


Let 
m(p) = max |z,(x,p)I, v = 0,1,...,(n — 1). 
OsxSl 


By using Lemma 2 to estimate the right-hand side of (25), we obtain the 
inequality 


le,(p) <| a ok fps +. + Pela 


+75 Cle —® [(lpel + + es s)de| me) 


Since the left-hand side attains its maximum m(p), it follows that 


m(p) < mo) [Ll bet ipta|at < < mo) F3 ; 


where C, is a certain constant. 

For large |p| this inequality is possible only if m(p) = 0; hence certainly we 
have z,(x,p) = 0. It follows from this immediately by (24), for v = 0, that 
y = 0. This proves the theorem completely. 

Remark. The equation 


U™ + p"U =0 
has, for p # 0, the fundamental system of solutions 
U, = eP Or k= 1,2,...,72 


From formula (17) we deduce that the solutions y, in Theorem | can be 
approximated by these functions U, for large values of |p]. 
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In other words, for large values of |p|, only the first and last terms in the 
equation 
y + poy"? +... tapay + p"y = 0 


play an essential part. 

In the case of a second-order equation of the form y” + p(x)y — p*y = 0, 
there are four of the domains S (fig. 9); by the theorem just proved there exist 
in each of them linearly independent solutions y, and y, which can be 
represented asymptotically in the form 


nef of) 
naefio()] 


For real and positive p it is often convenient to replace these solutions by the 
linear combinations of them 


Yi + Ve — cos (px) + (2), 
2 p 
v1 Je = sin (px) + 0(*). 
21 p 
The asymptotic behaviour of solutions of differential equations for large 


values of the parameter has been investigated in many works. Additional 
information and references can be found in the books [17], [97], [95b]. 


Fig. 9. 
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6. Refinement of the Asymptotic Formulae 


If the coefficients po,...,p, of the differential expression are not only themselves 
continuous but also have continuous derivatives up to a certain order m, then 
the asymptotic formulae already obtained can be made more precise. 

Suppose, for example, that the coefficient p,() has a continuous first 
derivative p,(é). 

We substitute in the right-hand side of (18 a), instead of the functions y, and 
their derivatives, the expressions from (17). Then using the estimates (15 a) 
and (15 b) for v = 0, we obtain: 


Vy = PO - | Kid pera *pylEd 


$+ | "Kylee E,petoat"plQdé + 0(5] ene’ (28) 
np J x Pp 


But by (12) 


k-1 x 
[" Kicegner>plOdg = cron ie eR 


+ 0 [, ps8), 


and by integrating by parts and using the inequalities (5) we can easily convince 
ourselves that the first addend in the brackets is an_expression of the form 


o(*). Hence 
p 
x x 1 
[Kent wer'plende = coon [ ralende + o(*)]. 
0 0 p 
and similarly 


1 ] 
J K,(x,€,per pol E)dé = e°"O (;) 
Substituting these expressions in (26), we get 


y,= oo 4+ “aol + (5) |. (27) 
where 


YuolX) = — = [ pues (28) 
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Similarly, by substituting the expressions from (17) into the right-hand side of 
(18 b), we obtain: 


\ 
AY, _ proter| 4 Yio) + (5). y=0,1,....4n —1); = (29) 
p 


in this case it can readily be checked that the first two terms in the brackets in 
formula (29) result from the fact that we differentiate the expression 


ePOKx | l + Jy of) 
p 


in (27) v times and, after taking out the factor p’w,e°°**, discard all terms of 
] 
order o(-:)- 
If p. and p; have continuous derivatives up to the second and first order 


respectively, then by substituting the expressions (29) in the right-hand sides 
of (18 a) and (18 b), we derive in a similar way to the above 


Repetition of the argument leads to the following 
Corollary. If the functions P»,P3,...,P, have continuous derivatives up to the 
orders m,(m — 1),(m — 2),... respectively in the interval {0,1}, then the follow- 
ing asymptotic formulae hold for the solutions y,yo,...,¥, constructed in 
Theorem 1: 


y= oon + aos + 2s +... + Honk T (=m) | (31 a) 


oY 
Vk = prayer 42a +e) +... + teal + o( )| (31 b) 


ax” p7t? 


k = 1,2,...,n; vy = 1,2,...,n — 1, 


where Vyy1(X),----Vivm(X) are continuous functions in the interval [0,1]. 


Remark 1. The formula (31 b) follow from (31 a) by termwise differentiation, 
taking the factor p’aye?°™ outside the brackets, and collecting all those terms in 
the brackets which contain factors \{p* with k > m into an expression of the 


l 
form o( =e] . 
We can easily satisfy ourselves that this is so, by considering that the term 


m+4+1 


or) in (31 a) is an expression of the form 
p 
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l 


x 1 
| [. e POX KI (x,€,p)p(E,p)dé — { 


eH PoM- OK (spol Ep |, 


where (€,p) denotes a certain bounded function; from the Lemma in section 5 
it then follows that the vth derivative of this expression with respect to x 


(v = 1,2,...,(7 — 1)) is of the form (st): 
p 


The functions y,, (x) can therefore be expressed in terms of the functions 
Vron(X)- 


Remark 2. We can determine the functions y,o,(x), up to additive constants, 
even more simply by substituting the expressions (31 a) and (31 b) into the 
equation /(y) + p"y = 0, and, after removing a factor p"e°°**, comparing 
terms with the same powers of p, from 1/pto 1/p™ inclusive. In this way we obtain 
a system of recursive differential equations of the first order, from which we can 
successively determine the functions y,o,(x), up to additive constants. 

In fact, after the substitutions mentioned, and the removal of the factor, we 
obtain the result 


A A+. +4 + 0(in) =0, (32) 
pp p p 


where A,/p” denotes the sum of all terms of the equation which contain 1/p” but 
not 1/p’*’. 
Obviously, equation (32) can hold only if 


A, = 0, Ag = 0,..., A,, = 0. (33) 


If we actually calculate the expressions for A,Ao,...,4,, We obtain from the 
equations (33) the following expressions for the functions yyo,: 


I x 
Yuo = 4% —_{ P(E)dé, (34 a) 
NW, 0 
yv-l v-a-1l x 
Yor —2 SY OL Pyare oO ¥ = 23,4, 84) 
a=0 p=0 0 


The constants a, (v = 1,2,3,...) could already be found from the equations 
(18). Thus we saw earlier that «, = 0, for example. In the same way the 
remaining constants «, can be determined. 


7. Normalization of the Boundary Conditions 


We wish to investigate the different systems U,(y),v = 1,2,...,n, of linear forms 
which define a given differential operator. If yf? or y{” appear explicitly in the 
form U(y) but y and y” do not, for any v > k, then we say that the form 
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U(y) has order k.1 We consider the forms U,(y) of order (nm — 1), if there are 
any. By replacing them, if necessary, by equivalent linear combinations, we can 
arrange that the maximum number of forms of order (n — 1) is <2. The 
remaining forms have orders <(n — 2); we apply the same process to the 
forms of order (n — 2) and reduce their number to a minimum; and so on. 

The operations described are referred to as the normalization of the boundary 
conditions, and the boundary conditions finally resulting are said to be 
normalized. From the way in which they are constructed it follows that the 
normalized boundary conditions must have the form 


Uy) = Uv o(y) + U'i(y) = Q, (35) 
where 

ky-1 

Uvo(y) = a) w) + »y a, jv ‘), (36) 
j=0 
ky -1 . 

Uyi(y) = By? + X By yy (37) 
j= 


n—-l>k, >k,2>... 2k, >O, Kao <k,, 


and for each value of the suffix v at least one of the numbers «,,8, is non-zero. 


8. Regular Boundary Conditions 


We consider a fixed domain S;,; as before, we number @,,w»,...,w, So that, for 
p E Sh, 
(pw) < A(pw.) <... < A(pa,). 


For our later purposes it is convenient to particularize a class of boundary 
conditions which we shall refer to as regular. This class is defined differently for 
odd and even values of n, as follows: 


(a) nis odd: n = 2p — 1. 
The normalized boundary conditions (35) are said to be regular if the numbers 
6, and 0, defined by the identity 


ole} Own y (x + SB,)wn Bios, By co, 
k k k k k 
Ogi)" py 1 (ot, + SPa)w,* Bye Bow," 
0) +05 = " so "| 38) 
k k 
O,W" ve XO 4 (a, + SB,,oon” nny B,con” | 


are both different from zero. 


1 yy, ¥, correspond to ya, ¥, of §1.3. 
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(b) nis even: n = 2p. 
‘The normalized boundary conditions (35) are said to be regular if the 
numbers 0_, and 9, defined by the identity 


Oy 
5 + & + Os = 


l 
k 
yw"... ows y (a, +58, )wi (a+ Ps) at Bywis » J. Byaog! 
| 
k 
AoW"... Ones? (ot. -+5B.)ar? (a2 Bo) t Bows,» Lee Bown? (39) 
] 
k n n n n 
Oy" ‘ee OL,” 4 (a, +58, us (o0-+ Pa) of Ban?» cee Bn 


are different from zero. 

This definition of regularity is independent of the choice of the region S for 
which the numbers w, were arranged in order. In order to show this, we first go 
back to section 2 and recall that there we were able to obtain the different 
ordered-arrangements of the w, belonging to individual domains S, by 
multiplying by an nth root of unity the values of the w, in the order in which 
they occurred for one or other of two basic, adjacent domains, say the domains 
S, and S,,,-,. If now all the w, in equation (38) or (39) are multiplied -by the 
same factor with modulus equal to 1, then the numbers @ will change by a 
common factor with modulus equal to |. In particular, the ratios @: 4, or 
6,:9,:8_,, and so also the roots of the equations 6, — @ =0 or 6,€? + 
Go€ + &_, = 0, remain invariant; we shall apply this result later. 

So it now remains to find out how the numbers @ change when we change 
from an appropriate domain S, with k even to a suitable chosen domain S,, 
with k’ odd. We shall, in fact, find that this change, too, can be achieved for 
each separate number @ by multiplication by a suitably chosen number with 
modulus equal to 1 but that now, in general, all the numbers @ will not be 
multiplied by the same factor. 

We discuss this first for the case when n is odd: 

If n =2u — 1, the numbers ¥, = et", k = 0,+1,42,...,4(u — 1) 
obviously represent all the distinct nth roots of —I. We can correlate any 
particular order of the w, to the ¥, simply by specifying the row of suffixes. For 
example, the row of suffixes in the proper order for the inequalities (4) to hold 
in the domain S,,_, runs: 


(0, +1,—1,+2,—2,...,(u ~~ 1), —(u ~ 1)), 
while that for the domain Sg Is 


(0,—1,+1,—2,+2,...,(—p ~ 1),(u _ 1)). 
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Now, we have 


k k k k 
XW,” eee yw! By," 1 eee Bw," 


k k k k 
a," 1... Ow," Brut ves nnn | 


| ad 


and the determinant for @, has the same form except that in the xth column the 
a, are replaced by f,. 

Let us suppose that 6) and 6, have been formed for the domain S4, _,, and let 
the corresponding numbers for the domain S, be denoted by 6, and 4). It is 
clear that 6,,6, are obtained from 6,,6, by interchanging the numbers w, in the 
second and third columns, in the fourth and fifth columns, etc. For ” odd, this 
reduces to an actual interchange of columns in the determinant for just one of 
the numbers @, or @,. This will become clear by considering the following 
arrays: 


(a,) w even, n = 4v — J. 


6,| a a a... oO a a fp... 8B B 
(0, 1,—1,....¢ — I),—-@ — 1),v,—,...,(2v — 1),—(2v — 1)) 
6, aa a... « ao Pp Bp... 8B B 
(a,) wodd,up=2v+1n=4 +1. 
6,| a a a..a a 8B Bp ..p 8B 
(0, 1,—1,...,»,—v,(v + 1),—(v + 1),...,2», —20)]. 
6. aa a..a B Bp p ..B 8B 


These arrays merely show, for each column of the determinants for 6, and 
6,, the #-suffix for the w, occurring in that column, and the letters « or B above 
and below this index shows whether one of the a,...,«, or one of the f,,...,8, 
goes with it. For the transition from S,,-, to So, the index row (0,1,—1,...) is 
simply replaced by the row (0,—1,1,...) and the (a,f)-distribution remains 
unaltered. 

This can be accomplished in case (a,) for 6,, and in case (a,) for 9, by a 
simple interchange of columns, and it follows that 


6, = —6, for n=4y —1, 
§6.= % for n=4v+1. 


To calculate the value of the other @ in each of these cases, we proceed thus: 
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A multiplication by e =e" transforms the system «,...,w, intoitself; only 
the sequential order is permuted—and this simply by adding / to each number 
in the index row. 

From the suffix row: 


a ao a Bp 8B 
0,1,—1,. a 1),—(@ — 1),v,—»,(v + 1),...,-(2v — 1)) 
the row: 
te oe, a o B OB.. 
(—1,0,—2,1,—3,.. ov — 9), —v,v — 1),—(@ + 1),y,. (20 — 2),2v — 1) 
can be obtained (with / = —1). 
We see that this last row can be brought into the form 
Oo a B B 


(0, —1,+1,—2,+2,...,—( — 1),(v — 1),-»,», -—@ +:12),0 + 1),..., 


—(2v — 1),Q» — 1)) 


by interchanging the members of its first (2» — 1) pairs. 
Hence it follows that 


6, __ — erri(kat weet k,)/n 8 for n=4y— 1. 


Similarly it can be shown that 
6, — ew 2riCkit «2.4 kndi/n 0, for n= 4p + l. 


These results imply that, for odd n, the definition of regularity is independent 
of the-choice of the zone S,. Moreover, it has been shown that: 

The root of the equation 6,€ — 6, = 0 remains invariant for the transition 
from S, to S, if k = k’(mod 2), but for the transition from S,,4, to Sy, it 1s 
multiplied by a factor e=?™("+---+*=)/"| where the plus sign is to be taken if 
n = 4y + 1, and the minus sign if m = 4v — 1. If €, € are the roots for S, 
with even and odd k, then £2 = e#274+..-+ king), 

For even n, n = 2u, we can attain our goal rather more conveniently. The 
numbers #, = e@t*™/", k = +1,43,+5,...,4(n — 1), now denote all the nth 
roots of —1, and so we can again characterize any special order of the numbers 
W,,...,@,, and consequently any particular domain S,, by giving the corres- 
ponding row of suffixes for its #,. For the domains S,,_1, So, Sy respectively we 
obtain in this way the suffix rows: 


( 1,—1, 3,—3, ’ (n ~ 1),—( ~ 1)) 
(—1, 1,-3, 3, »—(n— 1), @—1)) 
(—1,—3, 1,—5,3,—7,5,....-—(7—1), (@—3), (—1)). 
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We again distinguish two cases: 
(b,) pe = 2v, n = Ay, 


The row of suffixes for S, can be obtained from that for S, by interchanging 
the second and third columns, the fourth and fifth, ..., the (n — 2)th and 
(n — 1)th. 

If the corresponding columns are interchanged in the determinant (39) 
formed for the domain S,, andif sis replaced by 1/s, then clearly the determinant 
(39) for the domain S, is obtained. 

Since the number of column interchanges carried out is odd, we have 
D,(s) = —D,(1/s), where D,(s) denotes the determinant (39) for the domain S,. 


(b,) p=avt+in=4 42. 


In this case we obtain in exactly the same way 


Do, -31(5) = —Dy (1/s). 


From the equations obtained in cases (b,) and (b,) it follows, first of all, that 
the definition of regularity is independent of the special choice of the domain 
S;. 

We further notice that: 

The roots é’,” of the quadratic equation D,(s) = 0 for the regular boundary 
conditions remain invariant when k is changed to k’ with (k’ — k) even; but 
they go over into their reciprocals 1/€’ , 1/€” if (k’ — k) is odd. 

Finally we remark that: 

If the numbers «,...,«,,64,..-,8, are all real, then the numbers @ go over into 
their complex conjugates for the transition from S,,_, to So. In this case, then, 
the independence of the definition of regularity from choice of S, follows 
immediately. 

We shall now discuss a few examples of regular boundary conditions. 


(a) Conditions of Sturm Type for even n(n = 2p). 
This is the name given to boundary conditions of the form: 


U joy) = yh? + 5 ayo? = 
J = 1,2,...,p, (40) 
U,V) = VY + Byy? = 


withn — | Pky >k,>..>k, 20 n-loek>ky>...>k, S0. 
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Here half the boundary conditions contain only the function values of y and 
its derivatives at the point x = 0, and the other half of them only the values at 
the point x = 1. 

In this case, 


k k k k 
ws} Ory of wis 0 0 
k k k k 
O_, 4 Wy" wy w aw" 0 0 
— +9 + Os = + 
s 0 0 sak taki ki 
Oy, “Arr. Wyse Wy, 
k’ ] k? k? k? 
0 O sw out ee re 


(The + sign occurs because the rows are not in the same order as in formula 
(39)). Consequently 


0, = 0 
wy! wat 1 wa +1 wit was 2 wr! 
6,=+ . . , . . ees (41) 
wy" wit 1 Oty war wnt, 2 war 
wy? wy! wi, 1 we 
G6,=+ . Lee , . ves . (42) 
ws" wie whe lose win 


Hence conditions of Sturm Type are regular provided that all the determinants 
in the formulae (41) and (42) do not vanish. 
If, for example, 


k, =k, —1, ks =k, —2,...,k, =ky —e +1, 
ki =k, —1, kg =k, —2,...,k, = ki —v +1, 
these determinants are different from zero and so the conditions are regular. 
(b) Conditions of Periodic Type. 
This is the name given to conditions of the form: 


Uy =P —yP =0, =» =0,1,2....,(7 — 1). 
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Conditions of periodic type are regular. 


For, if n is even (n = 2), then 


I ae (1 —s) (1-4) l 1 


. S 
l 
6 ; Wy, eee Wi-l (1—s)w, l—- Wutl Wi+2 eee Ww, 
ee 
n-1 n-1 | n-1 l I n-1 n-1 n-1l 
Wy cee Wut ( —S)wy, _ Wi+2 Wite2 eee Wy 


=+Cc(1 — (1-4), 
5 
where C is the Vandermonde determinant for the numbers w),...,w,, and 
obviously does not vanish. 
Hence 
6, = +2C; 0, = 0, = +C £0; 


i.e. the conditions are regular, 
If nis odd (n = 2 — 1), then 


ad (1—s) a 
W1 eee WwW -1 ¢ —s)w W +] coe WwW, 
6, +6,5= + ° me = +C(1—s). 
n-1 n-1 n-1 n-1 n-l 
lw oe Oy (L—s)wy, ON EL see Wh 


We deduce that in this case, too, the conditions are regular. 
(c) Boundary Conditions for n = 2. 


The most general boundary conditions for nm = 2 have the form 


AYo + by + AVo + doy = 0, 
Co + AY + Co¥o + hy. = ‘ “) 
We distinguish the following cases: 
1. a,d, — b,c, # 9. 
By solving (43) for yg and y;, we can bring the conditions to the form 


Yo + Wo + a1 = 0, 
Yi + Heo + Moy, = 0. 
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Then 
Wy, We 
O_ 
—4+6,+ 65 = | —!_.,. 
5 SW, _— Wo 5 
S 


1.e. the conditions are regular. 
2. ad, — bic, = 0, |a,| + |b,| > 0. 
In this case we can transform the conditions (43) so: 
Yo + AY + ayo + boy, = 9, 
CoYo t Ay, = 0, 
(a, + shy)a, —(a + - b, wy 


6 
vs FF OS = 
Cy + sd, Co + do 


] 
= w1(b1Co + a,dp) (s + “| + 2(44,Co + byd))w; 


the conditions are therefore regular if b,c, + a,d) # 0. 
3, Qa) = 6b, =c, =d, = 0. 


Since the formulae (43) must be independent, in this case aydy — boc, # 0; 
hence the conditions (43) are equivalent to yy = 0, y, = 0. Hence 


1 1 
OL 
— +6 +645 = l =1 ~,, 
AY Ss = A) 
A) 


i.e. the conditions are regular. 
Hence the conditions (43) for n = 2 are regular in just these cases: 


l. ad, —_ bic, x 0; 
2. a,d, — b,c, = 9, |a,| + [by] > 0, byco + aydy # 0; 
3, Qy —_— b, = Ci — d, — QO, Ad — bolo 7 Q. 


9. Asymptotic Behaviour of the Eigenvalues 


The results of section 5 enable us to prove the existence of infinitely many 
eigenvalues for differential operators with regular boundary conditions and 
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to give asymptotic approximation formulae for the behaviour of the sequence 
of eigenvalues at infinity. It turns out that the principal terms in these approxi- 
mations are already completely determined by the numbers 65,8, or &,1,0_4, 
which are calculated from the prescribed boundary conditions by means of 
formulae (38) or (39). The further nature of the differential expression and 
boundary conditions does not affect them at all. 

We assume below (and for brevity we do not mention this again in the 
hypothesis for the following theorem) that the coefficients of the differential 
expression considered are continuous in the interval [0,1]. It should, however, 
be expressly mentioned that all the results of this section would still hold if we 
required only that these coefficients be summable in the interval [0,1]. 


THEOREM 2. 


A differential operator of the nth order which is generated by an expression of 
the reduced form described in section | and by regular boundary conditions has 
precisely denumerably many eigenvalues, whose behaviour at infinity is 
specified by the following formulae: 

For odd n such that n = 4v — 1, the following two sequences are obtained:? 


sn Inf? l 

A, = (—2kzi) f — ore + (7) (44 a) 
, sn In,é62 ] 

Ay = (2kni) 1 +- nee + o( 7) (44 b) 


k=N,N+1,...; 


and for odd n such that n = 4v + 1, the two sequences: 


; Ing l 
rN — n nN 0 a ? 
, = (2k7in) 1 + Ski + (75) (44 a’) 
n An n In,€'2 l 
MM = (—2krin) 1 — aE + o( 7) (44 b’) 


kK=N,N+1,...5 


here €% and & denote the two solutions defined earlier of the equation 
0,6 + 0) =O for S, with k odd and k even respectively; and N is a chosen, 
sufficiently large, positive integer. 


1In view of the ambiguous usage of the word ‘‘denumerable”’ in English mathematical 
texts, we have followed the author’s terminology. 
2 Here Ino~é is any fixed branch of the natural logarithm. 
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For even n, n = 2p, and 0§ — 40,0_, #0, the following sequences are 
obtained: 


Ai = (—])! a)2# Has (c:)| 
= (—1)"(2kn) {1 + (,4)} (45 a) 
N= (—1)"(2ken)*"{1 + H Tot + o()}. (45 b) 


where &' and €" are the roots of the equation 
O15" + O& + 6 = (46) 


for the domain So, i.e. are roots of the equation D,(é) =0, and where the 
upper or lower sign is to be taken according as n = 4v orn = 4v 4+ 2. 

For even n, n = 2p, and 0; — 40,0_, =0, the following sequences are 
obtained: 


= (—1)"(2ken) "1 + Hoe + o( an) (47 a) 
, ; ; | 1 
Ng = (<I kn) {1 = TBO (a) (47 b) 


k=N,N +1,..., 


where &€ is the double root, occurring in this case, of the equation (46) (relative 
to the @ for the domain S,). The signs are to be chosen in the same way as for 
the equations (45). 

In the first three cases, all eigenvalues of sufficiently large modulus are simple, 
but in the last case, all eigenvalues of sufficiently large modulus can be either 
simple or double. 

Proof. First let n be odd (n = 24 — 1). We consider a fixed domain T. 

Let the numbers w, be numbered in the order such that, for p € T, 


B(p + c)wy) <A((p + c)we) <... <A((p + c)a,). (48) 

We put 
~, =(p + c)e,, k = 1, 2,...,n. (49) 
The points f,,fo,...,0, lie on a circle of radius |p + c| and divide the circum- 
ference into n equal parts, each corresponding to an angle of an yr 1 


In the closed, right half-plane there lie certainly at most » points, since more 
than » points (Fig. 10) would lead to the contradiction 


pL > 7. 


T= 


7 
2u — | 


Hence, and by the inequalities (48), it follows that at least the first (« — 1) 
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Pky+t 


Phu 


Pk 


Prs 
Fig. 10. 


points /,,f>,...,),-1 Must lie in the open, left half-plane. Similarly, the last 
(u — 1) points 6,,43,Py+2)---5Pn Must lie in the open, right half-plane. In other 
words, we have: 


Bp) < 0, Rp.) < 0,...,A(B,-1) < 0, (50 a) 
AP +1) > Q, Bp, +2) > 0,...,2(p,,) > 0. (50 b) 


If p — © and remains in the domain T, then the left-hand sides of (50 a) and 
(50 b) tend to — and +-© respectively. 

For, A(6, -1) will only then not tend to — oo if the angular distance between 
P,-1 and the positive or negative imaginary axis tends to zero. But this would 
imply that, for sufficiently large |p|, the (4 + 1) points f,-1,/,,...,), lie on the 
arc 


T €, 


implying that, for arbitrarily small e > 0, 


27 


2 
7 + a va 


which is impossible. 
From the assertion just proved it follows that, if p > 00, p € T, then 


(a) e°! tends exponentially to zero if j <p 
(b) e? tends exponentially to infinity if j > p. 
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By Theorem | in section 5 there exists in the domain J a system of n in- 
dependent solutions ),,y.,...,y, of the equation /(y) + p"y = 0 such that 


Ve = eon + o(*)|. 
p 


Ex | 
yr = per jo + o(*) |. (51) 
ye-) _ prrerons| af + 0()]. 


With these solutions we form the determinant A(A) = det[U,(y,)]; as was 
proved in §2.1, the eigenvalues are the zeros of A(A). 
So we shall now concern ourselves with this determinant A(A). 
For brevity, we introduce the notation [a] = a + O ! . 
By substituting the expressions (51) into the normalized forms U,(y), we 
obtain the equations 
l 
voy ;) = (poo, + o(*)| = (pw ;)*[o,], 
I © 
Uno) = (ow jrere{8, + 0(2)| = GupherByh 
and hence 


U(y;) — Uvo(¥;) + Uv(y;) — (pw ;)**{[oy] + eP [By }}. 


If 7 <p, the function 
eP?1 == e~ CsePS 


decreases exponentially as p —> 00, p ET; hence 


U(y;) = (pw ;)**[a,] forj <p. (52 a) 
Similarly we find 

U,(y;) = (pw)? [B,] for j > pb. (52 b) 
Finally, we have 

U(y,) = (pm,){[o,] + e?°*[B (52 c) 


We substitute all these expressions in the equation 
A = det [U,(y)] =9 (53) 


and divide out the common factors p™,p",...,0°" of the rows and also the 
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common factors e?®#*1,e?®#+3,... ef" of the last » columns of the determinant 
D(A). The equation can then be written in the form 


A, = 0, * 
where 
[o]oy vee CA ({o] + [B,Je??*)eoit [Bion tee [Bion 
[ogjot? ... [oglok, ([oo] + [Bey [Bolo .. [Belon 


A, = _ (54) 


[doy [oxen a (onl + (Ble? on” [Bloaes + [Balen | 


From the definition of the numbers @, and @, (see (38) in section 8), we therefore 


have 
Ay = [8] + e?°+[8,]. 


For a root p of the equation A, = 0 with a sufficiently large absolute value, we 
have 


je. eH = 2M ol) = at + o(;)] 7 ef t o(7) ©) 


since, by the regularity of the boundary conditions, 6, # 0 and 6, # 0. Hence 


! 
p 


1 


Wn 


fingé + 2kri + o( )}. k = 0,+1,+2.,.... (56) 
We now show that zeros of the function AQ) which are given by formula 
(56) do actually exist. 
We put 


] , 
py, = — (2kai + In,£); (57) 
p 


w, 
then equation (56) takes the form 
1 
p =r + 0(*), k = 0,+1,+2,.... (58) 


The numbers p, obviously lie on a parallel to the bisector of the domain T, 
so that for sufficiently large & there is only one choice possible for the sign of the 
integer k if the point p,; 1s to lie in 7. Closer investigation shows that: 


(a,)n = 4 — 1. 
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For a T with even subscript, k must be chosen to be positive; and fora T 
with odd subscript, k must be negative, 


(a) n= 4 +1. 

For a T with even subscript, k must be chosen to be negative; and for a T 
with odd subscript, k must be positive. 

About each point p, we now describe a circle I’, with the same radius r in 
each case. By what we have just said, for a sufficiently large |k| these circles will 
lie entirely within the domain T. The equation (53) is equivalent to equation 
(55), and since £ = e’*®# the latter can be written in the form 


e@ulO~ Pk) 0(-) —(. (59) 


Outside the circles T, the modulus of the function 
f — eu? — Px) __ | — e@ule — Po) __ l 


is always greater than a positive constant. 
For, if we introduce a new variable ¢ by putting 


w,(p _ Po) — ¢ 


f=é&—-1, 


and the circles [, transform into circles Ty, with the same radius r about the 
points € = 2kzi. Since f = f(%) is a periodic function with period 27i, we have 
only to show that its modulus is bounded from below by a positive constant in 
the domain D bounded by the lines %f = +7 and the circle [9 (Fig. 11). In 
this domain the function f(g) never vanishes, however, and for sufficiently 
large |2C|, say |#l| > N, |f(Q| remains greater than a positive constant, 
because 


then 


lim [f(Q| = co 
{> +0 
lim = [f@®| =1. 
[> -@ 


7/7 
Co 


Fig. 11. 


70 §4 ASYMPTOTIC BEHAVIOUR 


This proves our assertion. We conclude from it that, for sufficiently large p 
the function A has no zeros outside the circles I’,. 

Let m be the minimum of the function |e?*?— “Pm i] on I’,. Since 

— p, =re® on IT,, m does not depend on k. For sufficiently large p, 
10 (Ip) <monl,. 

Hence, by Rouché’s theorem (see, e.g. [91]), it follows that, within I’,, 
the equation (59) has the same number of zeros as the equation 
eeu Px) _ | —= 0) i.e. just one zero, which we denote by p,. 

By (56), 


pi = —{2kmi + Ingé + o(=)}: 
Wy p 


(4) =o) 


but by the same formula, 


We have therefore 


n= pte +ing +0(!)) 
Ww 
or 
, Kari Inoé 
rian ot + oe. (7). (60) 


If we apply the same arguments to a domain 7, first with even subscript, then 
with odd subscript, and raise the result to the nth power, then, taking account 
of the proper sign of the number k, we obtain for each case the desired sequence 
of eigenvalues (44). 

That these eigenvalues are simple for sufficiently large k rests on the fact 
shown above that they are simple zeros of the determinant A(A). 


Now let n = 2p, i.e. is even. We again consider a fixed domain 7 for which 
the inequalities (48) hold. By the same considerations as were used for n odd, 
we conclude that 

Kp) < 0, Kp) < 0,...,A(py—1) < 0, (61) 
A py +2) > 0, A(py+s) > 0,...,A(p,) > 0, (62) 
and the left-hand sides of (61) and (62) tend exponentially to —coo and +00 
respectively, as p —-> 00 while remaining in the prescribed domain T. 
Hence, we conclude, exactly as in the case of n odd, that in the domain T 
Oy ,) = (pw ;)**[a,] forj7 <p —1, (63) 
U,(y;) = (pw ;)""e?°[B,] for j > w + 2, 


and moreover 


Uy) = (po,)"{[o,] + e?[B,]}, (64 a) 
U(Vu+1) = (pea 41){[a,] + 67% *[B,]}. (64 b) 
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When n = 2p, if w, is a root of the equation w” + | = 0, then so is —w,. 
Hence, and by the inequalities (48), it follows that 


RG) < 0, A641) Z Q, 
and 


Wy = Way Wo = Wye yy sey My = Wy 4d 
Hence (64 b) can be re-written as 
Uu+1) = (py 42) {l%) + e7?“[8,}. (64 b’) 


By substituting the expressions (63), (64) for U,(y,) into the equation A = 0 
and removing the common factors, we obtain an equation of the form 


A, = 0 (65) 
where A, = 
[ows]... [xo 3) wr{[oy] + e? +B} wi i{[o,] + 27 ?°*[B,]} [Bw a]... [Bion' 
[ww]... [age uJ wi{ [oo] + e°?+[Bo]} wit i{[%2] + e~ Pu Bo)} [Bow n'y 2] .-» [Bown'] 


[oh]... [ook a) wir{[on] + e?*[B,)} war allo) + e7?°*T BD} [Baeu's ol --- [Ben] | 
(66) 


By the definition of the numbers 6 ,0,,8_, (see (39), section 8), 
Ay = [4] + [8Je°* + [8-yJe“?°* 
and therefore 
Pu, = [8,Je2™ + [O]e?* + [01] = Ae°P* + Oye? 4+ O_, + o(*), (67) 
because, by virtue of the relation 2(p,) <0, we clearly have 
jerrn| = fem] Je] < Je], 
i.e. the function e°”* is bounded in the domain 7. 
We consider the quadratic equation 
0,€% + Oé + 0_, = 0; (68) 
let €’, €” be its roots, so that 
O,£7 + OF + O_, = ACE — €') (E — &"). 
Then (67) can be rewritten as 


ePuNy = Oy(eP — £") (oP — £") + o(7) 7) 
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The equations 
eP@u _ £' — 0 and ef? — &" = 0 


have respectively the roots ‘ 


an, €’ + 2kni), a (Inp-€” + 2kvi), kK =0,+1,+2,..., 
Wy Wy, 


but we are concerned only with those which lie within the domain 7. 


Fig. 12. 


Now the roots of both these types obviously lie on a parallel to one or other 
of the boundary lines of the domain S. Consider the domain S), for example. 
Then we have 


oO, = —je-**/" for n = 4p, 


w, =! forn = 4 +4 2, 


i 


and the above roots are given by 


jeinin Ino é" _ 2kn, jeitin Ing é" —Ikn forn= 4y, 
and by 
—i Ing €&’ + 2kzx, —Iln,é° + 2kn forn = 4-422. 


Hence it follows that, for n = 4v all roots with k < 0, and for n = 4v + 2 
all roots with k > 0, are situated inside 7, and at a positive distance from its 
boundaries, if Ty is suitably chosen, i.e. if c is suitably chosen. 

In the first of these cases, we replace k by —k and So arrive at the two 
sequences 


? J ? e “ l 4 . 
On Wu 


which, for k = 1,2,... lie inside T>, provided that the upper sign is chosen for 
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n = 4v and the lower sign if n = 4v + 2. Using p’, p” and (69), the equation 
Ay = 0 can be put into the form 


[eeu - Pi) _ 1] [eeule - Pid) —1}+ (5) —(). (71) 


About each of the points p,, p;, A = 1,2,..., we describe a circle T; or 7 with 
the same radius r. If we choose r to be sufficiently small, all these circles will lie 
entirely within 7). If we apply Rouché’s theorem again, as we did for n odd, 
then it follows that for sufficiently large |p|, p € T>, the equation Ay = 0 can 
have zeros only within the circles I’, and Ty, and that it has the same number of 
zeros, in fact, as the equation 


fe@ule - Pi) — ]] [eule ~ Pic) —1]=0 (72) 
has. 

Now suppose that 65 — 46,@_, #0; then ¢’ # é”. Consequently the 
numbers p,; and p;, are pairwise distinct. Hence, by choice of a sufficiently small 
common radius r it can be ensured that the circles [; and [’; do not touch or 
intersect each other. 

In each of these circles the equation (72), and therefore also the equation 

o = 0 for sufficiently large k, have exactly one zero, which we denote by 


Px OF Py. 
In the circle I’; the modulus of the factor e®#”~-*? — 1 remains greater than 
a certain positive constant; hence by (71) it follows that in the circle I, the 


equation A = 0 is equivalent to the equation 


Hence 


and similarly 


~) dkri Iny €' ] 
ee o(;:)} 
and similarly for px. 

On raising this to the nth power, we obtain the formulae (45). 

It is important to notice here that consideration of a domain 7, with k odd 
would yield no new eigenvalues. For if we consider, for n = 4y, instead of the 
domain T, a suitable domain, for example, 7}, then E’ and &” go over into 1/€’ 
and 1/é”, and so the logarithms Ing’ and Ingg" become —Inoé" and —Iné” 
respectively, apart from a possible addend 2k’zi. On the other hand, w, is also 
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to be replaced by —w,, as may easily be seen from the considerations of 
section 8. Hence, in order that the numbers — — l an, E' + 2kmi) shall lie in T,, 


it is now necessary that k be positive, and this implies that in formula (70) for 
n = 4y the plus sign now stands. 

Hence the p; and also the p; remain unchanged. Since the equation (71) also 
remains unchanged (up to multiplication by a constant factor), the same p; and 
p, result. The same thing holds for n = 4v + 2 1f we go over from the domain 
T, to T,,-1. Hence, for large |A|, the eigenvalues (45) are the only eigenvalues 
of the operator L. 

Now suppose that 65 — 40,0_, = 0. Then we have €’ = €”, and so py = px, 
and the circles I’; and [, coincide for each k. The equation A, = 0 therefore 
has, for sufficiently large |p|, precisely two zeros in each of these circles 
lr, =, = ly, though these may possibly coincide to form a double root. 

Let 6, be either of these zeros. The equation (71) then takes the form 


[ecu - P%) __ 1]° + o(*) — 0. 
p 
It follows that 
] 
ePvle— Pe) | 4 o(—=) = 0, 
le} 


and hence 


| 
The last part of the argument as in the case 65 — 46,0, # 0 enables us to 
derive the formulae (47). 
Remark. If the coefficients p,(x),...,p,(x) have continuous derivatives up to a 
definite order, then more precise asymptotic formulae can be derived, involving 
higher powers of 1/p (see, for example, [54] and [110]). 


10. Asymptotic Behaviour of the Eigenfunctions 


We shall now use the results of sections 5 and 9 to elucidate the asymptotic 
behaviour of the eigenfunctions for eigenvalues with very large moduli. 

We again assume that L is a differential operator in the interval [0,1] and is 
generated by a differential expression with continuous coefficients and by 
regular boundary conditions. 


Let )1,Vo..+-,Yn be linearly independent solutions of the equation 
Ky) + p"y =0 which satisfy the relations (17) in a certain domain T (ef. 
Theorem 1, section 5). An eigenfunction which belongs to a prescribed eigen- 
value A = —p” withp € T must be expressible as a linear combination of the 
functions )y,Va5--+5 Vn: 

Y= OVy + Cao Fee + Cans 
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where the coefficients c; are non-trivial solutions of the system of homogeneous 
equations 


Cope + U,('2)C2 + eee + CV Vn)Cn — Q, v= 1,2,...,n. 


For simplicity, we consider only a simple eigenvalue A, for which the rank of 
the determinant A = det [U,(),)], vA = 1,...,n, is equal to (n — 1). 
Then 


V1 }’2 yn 
U91) U2(¥2) U()’,) 
y= (73) 
| U0) Uy 2) Uy n) r 


is an eigenfunction belonging to the eigenvalue A if we assume that not all the 
minors of the elements of the first row of the determinant A vanish. (If the 
latter were the case, we should have to arrange the ,,y.,...,y, In a row, of 
which not all the minors vanish). 

We examine the cases for n even and n odd separately. 


(a) Letn be odd (n = 2p — 1) 

In (73) we substitute for y, and U,(y,) their representations (17) and (52) and 
divide through by the inessential factors p*,p",...,p°",e°@#*1,...,e°°". Then an 
eigenfunction y, is obtained, perhaps not normalized, for which the following 
formula holds: 


eP TT] 16. eh @H-P 7] ePu*T1] eP Prt - DY], erent Di q] 
[wp ]wy? vee [w.Juit 4 (CA + epee [Bo]wit) [Bolen see [Bolen 
Yo = 
I fa,jot ... [oJokny (fo) +e°*[B los) [Biloir  -  [Balen” 
(74) 


Now for sufficiently large |A| the number p must correspond to one of the 
p, or p, lying in T; we therefore obtain for a domain T with an odd subscript 


] 


Lt 


and for a domain J with even subscript 


p = ph = | Deni + Ing £ + (7) 
Oy 


for a suitable positive k; and in these formulae the upper or lower sign is to be 
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taken according asm = 4v — l orn = 4v + 1. We are therefore led to define 


\ 
] , 
pp = - (F2kmi + Ing &), py? = 7” (4:2kmi + Ing €), 
Lu rm 
in each case; and we then have 
' I ” l 
Pe = Pe + o(7), Ph = PR + (7) 
Then obviously 


‘ (1) 
ePOux — ePkOux — e@uPy *T1], 


el) = 1 + o() + (o(z)) +..=14 (7) = [1]. 


In particular, 


since 


ePOn — etrknt + Inger 7) _ €[1]. 


By substituting these expressions in (74) we obtain the following representa- 
tion for each of the eigenfunctions y{?, y belonging to the eigenvalues 
Axor, Of formula (44) 


Yo = 
PTT] 1. ePH PRL] HPA TT] OHH PR DED] Le omPk DET] 
[owp]oy? ee [oegheogt ey [w. + EBolw,z? [Bolo 4 ve [Bolo 

| [o, Jor” see [a Jon [~,, + EB,Jwx" [Balog a see [Balon” 


(75) 
k=N,(N+1),..., o =12. 


If we restrict ourselves to a closed sub-interval x’ <x <x” inside the 
interval [0,1], formula (75) can be simplified considerably. For, if we take 


e#"””* outside the determinant, then the th elementof the first row is [1], and 
the other elements of this row become: 


ev 2077) for v < wand e& Onde DTT] forv > p, 


and these expressions tend exponentially to zero when k —> ©. 
Hence 


Olyewy or owt 1 Bowhs, Love Bow, , 
_ (c) ] 
yO =(— 1) Berne ae + (5). 79 
p 
k,, . n wn 
OX, Wy ore 06,0 1 Bw, +1 ss: B ws 


It must be noted here that, in general, the numerical order of the numbers 
W4,.++,W, Will be different for the sequences y and y®, 
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(b) Let n be even (n = 2p). 


We repeat the above argument and obtain two sequences of eigenfunctions 


corresponding to the eigenvalues Aj and A, To the eigenvalue Xj, corresponds 
the eigenfunction 


Vi, = 
eM PTT] 2. ee PTL] eH TT) oe eee? TT]. em OLT) 
CAC oe. [x,]wr , [o, + E’BoJw i Ee + = = Bee Out] [B.Jwi?, » .. [Bales 
/ ] 
[,]o4" —_ [o,Jo,” 4 lo, + g B Jon” |, + é P| Ou" [Blo 2 .[B, Jeon” | 
(77) 
k = 1,2,3,.... 
For any closed sub-interval of the interval [0,1] this formula simplifies to 
Vi = 
| I 
Opus? ... Oger? 4 (x + z p wry Bows, » J. Bow, 
, ] 
— 1) \# —1,@uA;,% Of — 
(—1)4~*e | + ( i) 
| 01,03" vee yOu” («, + é! p) wiry Bon", » ree Ben” 
Ogi? ... Onwwi? y (a, + E’By)wi? Bows, 2 J. Bowon? 
] 
+(—1)te 7 OHPK* . . . : |- (7) (78) 
06,004" see 04,0" 4 (a, + gE" B,, Jeon wane vee nn” 


The formula for the eigenfunction y,, corresponding to the eigenvalue 4; is 


obtained by replacing p, and €’ by p; and &", | 
For n = 2, the formulae are particularly simple. If, for example, the domain 


T is the first quadrant of the p-plane, then we have w, = iandw,,, = —i, and 
formula (77) then has the form 


ph, I , — ip;.x s ] 
Vn, = (—iytelt + 7 By + (7) — (i)e k Ee +€ Bo + o(7)|. 


The results presented in this section are derived essentially from Birkhoff’s 
work [5]; the method of proving Theorem | stems from Stone [109a]. 
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Theorem | has been proved by Birkhoff under more general conditions; viz. 
the coefficients p,, k = 2,3,...,n can be analytic functions of the parameter p 
which are representable in the neighbourhood of infinity by 


. 


PAX,p) = » Ayy(x)p” 


y=0 


(and this representation may even be taken to hold in the sense of asymptotic 
convergence). The domains S and 7 also become more general, and instead of 
formula (17) the following formula is obtained: 


VW. = je(o) + o(*)| eXp f [ euoad, 


where w,(x),...,w,(x) are the roots of the equation 
w" + a,-1 (xo * +... + doo(x) = 0, 


and u,(x) are certain functions of x which are defined by formal substitution 
for y, in the differential equation. 


11. Various Generalizations of the Asymptotic Estimates 


(a) Arbitrary Interval. 


An operator L defined in an arbitrary interval [a,b],a < b, can always be 
brought, by the transformation x = a + #(b — a), into a form for which the 
above asymptotic formulae are applicable, provided only that the differential 
expression for L satisfies the conditions required in section 1. Correspondingly, 
in the formulae (17) and (75) to (78) x is to be replaced by (x — a)/(b — a), 
and in the formulae (44) to (47) for the eigenvalues the number 7 must be 
replaced by z/(b — a). 


(b). The Boundary-value Problem L(y) = Apy. 


Let L be an arbitrary differential operator for the interval [a,b] and in whose 
differential expression /(y) the coefficient po(x) is =1. Suppose that the function 
p(x) 1s real and continuous and has the same sign throughout [a,b]; without loss 
of generality, we may take p(x) > 0, for if this is not so initially, we have only 
to replace p and A by —p and —A. We introduce a new independent variable 1 
by the substitution 


1 P® woe 
c=] Vol at, (79) 


where 
h = [ V e(€) dé. (80) 


| With the variable ¢, the boundary-value problem L(y) = Apy can be written 
in the form Lyy = Ay, where L, is the operator which is obtained from J. by the 
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change of variable to ¢ and division by p(x). For its differential expression the 
coefficient of the highest derivative is again equal to 1. 


Then by a change of the dependent variable to y, where 


y = yexp (—+ f ecoar), 


we obtain finally from L, an operator L, for which all the assumptions of 
Theorem 2 are satisfied and accordingly the previous asymptotic estimates 
(17), (44)-(47), (75)-(78) may legitimately be applied, provided that the con- 
dition for regularity of the boundary-conditions, and the incidental assumption 
that the eigenvalues are simple, are satisfied. If we now transform these 
formulae back to the old variables, we obtain the corresponding formulae for 
the boundary-value problem L(y) = Apy. In particular, the formulae for the 
eigenvalues are obtained from (44) to (47) by replacing 7 by 7/{[2./p(x) dx}. 

Obviously we have also arrived at asymptotic approximations for any 
eigenvalue problem L(y) = Ay of the type considered in Chapter 1, for which 
the condition p)(x) = | is unnecessary. 

The important assumption, of course, here is that p,(x) is real in [a,b], does 
not vanish, and does not change sign. 

Considerable complications arise in the case where the function p(x) 
vanishes at some points of the interval [a,b], or, what comes to the same thing, 
if the coefficient p,(x) in an eigenvalue problem L(y) = Ay vanishes at certain 
points of (a,b). This case and some of its various generalizations were the 
subject of investigations by R. Langer and other authors (see the paper by 
R. Langer [57e], which deals with the case n = 2). 

The most general results in this direction have been obtained by V. S. 
Pugachev [92]. 

(c). J. D. Tamarkin [110] obtains asymptotic expressions for the eigenvalues 
of the generalized boundary-value problem 


y” + p,(x,A)y"? + eee + p,(x,A)y — 0; Uy) — 0, v= 1,2,...,n, 


where p,(x,A) is a polynomial of the kth degree in A, and the coefficients of the 


forms U,(y) are polynomials of the nth degree in A. 
Even more general results (for differential operators with partial derivatives) 


have been obtained by M. V. Keldysh [47]. 


§5. Expansion in Terms of Eigenfunctions 


1. The Basis of the Fourier Method 


The solution of partial differential equations by Fourier's method leads toa 
most important problem: the expansion of a given function in terms of the 


eigenfunctions of differential operators. 
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Suppose, for example, that a solution of the equation 
aN 


Cu =u 
ae ax + pe) a+, . + p,(x)u, axx <b, (1) 
is to be found which shall satisfy the initial conditions 
Ou 
ose; |B] = 900s 2) 
t=0 


and the boundary conditions 


n-1 Ou n-1 O’u 
2X an (2) _, +d Balas), _ 0, J = 1,2,...,72. (3) 


We shall seek a solution of equation (1) which satisfies the boundary con- 
ditions (3) and has the form 


u = y(x) (A cos pt + Bsin pt). (4) 


By substitution in (1) and (3), we find that the function y(x) must satisfy the 
differential mane 


Ky) = + ro) et 7+. + p,dy = —p*y (5) 


and the boundary conditions 


Uy) = by ox jy + y, B, yy = — (6) 


Ifnow y # O, then y is an eigenfunction of the boundary-value problem (5), 
(6) for the eigenvalue —p?. Let 


—Pi;—P2s—P3>--- 
be all the eigenvalues of this problem and 


y(x), Yo(x), Y3(X), ee 


the corresponding eigenfunctions, where each eigenvalue is repeated as many 
times as it has linearly independent eigenfunctions belonging to it. Then the 
infinite series 


(oe) 


= > »y,(x) (A, cos p,t + B, sin p,t) 


satisfies, at least formally, the equation (1) and the boundary conditions (3). 
It still has to satisfy the initial conditions. Substitution into the first of the 
initial conditions yields 


fla) = ¥ Agval2 (7 
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This equation represents a series expansion of the prescribed function f(x) in 
terms of eigenfunctions of the boundary-value problem. 

The question as to the precise basis of the Fourier method thus leads directly 
to the following problem: under what conditions may a prescribed function f(x) 
be expanded as a Series of eigenfunctions of a proposed boundary-value problem? 

The solution of this problem is easiest in the case of a self-adjoint boundary- 
value problem, i.e. if the expression /(y) and the boundary conditions 
U,(y) = 0, j = 1,2,...,n, generate a self-adjoint differential operator. 


2. The Case of a Self-Adjoint Operator 


Let L be a self-adjoint operator which is developed from a differential ex- 
pression (y) and the boundary conditions U,(vy) = 0, j = 1,2,...,2. Without 
loss of generality we can assume that Ly = 0 only for y = 0, i.e. that the 
boundary-value problem 


(y)=90, Uy) =0, fF =1,2,....0, (8) 


has only the trivial solution y = 0. For, otherwise, we have only to replace 
l(y) by the expression /(y) — cy, where cis a real number different from all the 
eigenvalues of the operator L. And we know that such a number exists because 
a Self-adjoint operator can have only denumerably many eigenvalues. 

If the boundary-value problem (8) has only the trivial solution, then the 
operator L has a Green’s function (§3, sections 1-6), and this 1s a Hermitian 
kernel. 

We consider an arbitrary function f(x) in the domain of definition of the 
operator L: this means that the function f(x) has continuous derivatives up to 
the nth order and that it satisfies the boundary conditions (6). If we put 


Lf =h, 
then 
fx) = f° G(x,Enla, 


i.e. the function f(x) is “‘sourcewise” representable’ by the continuous kernel 
G(x,8), | 

By reason of the Hilbert-Schmidt expansion theorem in the theory of 
integral equations (see [86a], for example), the function /(x) can be expanded 
in a uniformly convergent series of eigenfunctions of the kernel G(x,¢). But 
we have seen that the kernel G(x,é) and the operator L have one and the same 
set of eigenfunctions. This proves the following theorem. 


1 This terminology is commonly used by Russian authors where a function is expressed 
using a Green’s function as kernel. 
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THEOREM | 

Any function in the domain of definition of a self-adjoint differential operator 
can be expanded in a uniformly convergent, generalized Fourier series in terms 
of eigenfunctions of this operator. 

We remind ourselves that the domain of definition of a differential operator 
of the nth order consists of all functions which have continuous derivatives up 
to the nth order inclusive and which satisfy the boundary conditions belonging 
to the operator. The above assertion about expansion in eigenfunctions 
therefore applies to all functions having these properties. 

Since the domain of definition of the operator L is dense in L*(a,b), we 
deduce from Theorem | the immediate 


COROLLARY 


The eigenfunctions of a self-adjoint differential operator form a complete 
system in L*(a,b). 

If the eigenfunctions y,(x), yo(x),... are chosen so that they form an ortho- 
normal system, then the coefficients A, in (7) are determined by the formulae 


A, =(6.9n) = [, fOdrC)dx, 2 = 1,2,3,... 


These coefficients A, satisfy Parseval’s relation 


[2 yee = ¥ Lay 


Theorem | remains true in the case of the general boundary-value problem 
Ky) —Aox)y =0, UY) =0,  v = 1,2,...,n, 


where /(y) is self-adjoint and the boundary conditions are self-adjoint, pro- 
vided that the function p(x) is continuous and positive in the interval [a,b]. 
For, in this case, the eigenvalue problem L(y) = Apy is equivalent to 


(x) =A |” Gx pOyOae, 


and the integral operator occurring here has the kernel G,(x,¢) = G(x,£)p(&) 
and is a Hermitian operator in the space & of the functions f (x) with the scalar 
product 


(fof) = | ACfeedax, (9) 


le. G,(x,€) is a so-called symmetrizable kernel. 
So the Fourier coefficients A, are now given by the formulae 


A, = (f9n) = [fT o(adx, 
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where the y,(x),yo(x),... form a complete system of eigenfunctions of the 


1] , 
operator 5 L, orthonormal in the sense of the scalar product (9). 


Kamke [45] weakened the condition p(x) > 0, replacing it by p(x) > O with 
another restriction. He also generalized Theorem | to the case of a self-adjoint 
boundary-value problem 


L(y) + Al,(y) = 0 for any real A, U,(y) = 0, » = 1,2,...,n, 


under certain additional assumptions relative to the differential expression 
l,(y) and the boundary conditions. 


3. Expansion in terms of Functions of a Differential Operator with Regular 
Boundary Conditions 


The considerations in §5.2 do not attain the objective if the operator L 
is not self-adjoint. So we shall now adopt another method which depends on 
the analytic properties of the operator L — Al and on the asymptotic formulae 
obtained in §4. We shall assume that the differential expression from which L 
is generated has the special form described in §4.1 and that the boundary con- 
ditions are regular in the sense of the definition laid down in §4; we shall 
further assume’, as in §5.2, that A = O is not an eigenvalue of the operator, 
and therefore L has a Green’s function G(x, 6). 

In the complex A-plane we consider a sequence of circles I’, (kK = 1,2,...) 
with the origin as their common centre, and having the following properties: 


1°. The radius R, of the circle I’, increases without limit as k — oo. 


2°. There is a positive number 6 such that all the eigenvalues of the operator lie 
at a distance >5 from each of the circles I°,. 


By virtue of the asymptotic properties of the eigenvalues proved in §4.9 such 
circles I’, do exist. Let G(x,s,A) be the Green’s function of the operator L — Al, 
and in particular, let G(x,s,0) = G(x,s) be the Green’s function for the 


operator L. 
We apply the theorem on residues to the integral 
1 G(x,s5,A)dr 
L,==— ————— 
2ri lx A 

and obtain 

Mi H , 

Ty = Glxs) + 9, (10) 
v=1 y 


1 It follows from the argument at the beginning of §5.2, and from the asymptotic 
formulae for the eigenvalues that this assumption does not restrict the generality. 
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where H,(x,s) is the residue of the function G(x,s,A) at the pole A, and m, 1s the 


number of such poles in the circle I’,. \ 
We shall prove that 
lim J, =0 (11) 
k—o 


uniformly with respect to x and s in the interval [a,b]. The equation (10) will 
then lead to the series expansion 


with uniform convergence for all x and s in [a,)]. 

In fact, it follows from the asymptotic formulae for the eigenfunctions that 
the circles I’, can be chosen so that m,4; = 442; onthe other hand, it is easily 
seen that, in the squarea <x,s <), 


H,(x,5) 
A, 


tends uniformly to zero as k — o. 


The proof of the relation (11) is based on 


LEMMA | 
On the circumference YT, the function G(x,s,A) satisfies the inequality 


M 
|G(x,5,A)| < Ini i~ in (12) 


where M is a certain constant. 
Proof. We put A = —p”. Then, with a suitable choice of arg p the circumference 
I’, transforms into the arc y,: 0 <arg p < 2zx/n, |p| = ./|A|; y,, varies entirely 
within the two neighbouring domains S) and S, of the complex p-plane. 

In order to prove the inequality (12), we make use of formulae (32)-(34) of 
§3.7; 1t 1s convenient to deal with the cases n even and n odd separately. 
(a) Let n be odd; n = 2p — 1. 


Let the numbers w,,w9,...,w, be numbered in the order such that, for p € So, 
B( pa) < B(pw.) <... <A(pw,). 


Then, for p € Sy, we have 


0 
0. (13) 
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We denote the two parts of the arc y, which lie in S, and for which 
#(pw,) <0 or A(pw,) >0 
by Vp and vy, respectively (see Fig. 13). 


We evaluate the function G(x,€,A) on the arc y,, making use of the formulae 
(32)-(36) of §3.7. 


Fig. 13. 


We denote by W, the cofactor of the element y("~”). in the determinant 
yr PE) WP PME) PE) 
yer P(E) vB) we PE) 


—_ 
— 


! yy(é) yol€) YnlF) 
and put 


W(8) 
We U4) 


Then formula (32) of §3.7 can be written as 


z,(é) = 


o(x,é) = +4 y y,(x)2\(2). (15) 


By Theorem 1 of §4.5, we have, for p € So, 
yO(E = Peep "oi, j = 1,2)... v = 0,1,2,...% — 1). 


We substitute these expressions into (14) and cancel the factors 
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p,p2,...9p” 7, eP™, ePxS  ePnS in the numerator and denominator, so 
obtaining . 


B 
z,(£) = e7 Per —_ (16) 
p""* B 
where 
wit wy? wit 
n-2 n~2 n-2 
Wy Wo Wy, 
p= ; 
l ] J 


and f, is the cofactor of the element w?~* in this determinant. Hence 
n 0 for j = 0,1,2,...,(n — 2), 
Y od = (17) 
v1 fs 1 forj =(n — 1). 


The system (17) has a unique solution, for otherwise it would be satisfied by 


By _ WwW, 
p n° 
since w, = —1. Consequently (16) takes the form 
] 
2() = errs [-wh v= 12,1 — I), (18) 


From formula (15) and formula (32) of §3 it follows that 


Ue) = 2¥ Ualo)2(O + 4¥ Uno de8). (19) 


We consider the function G(x,€,A) for x > €; then the sign in (15) must be +. 
We multiply the first » columns of the determinant H(x,é,d) in the numerator of 
formula (34) of §3.7 respectively by 42z,(&), $2,(§),...,42,(6), and also the 
(u + 1)th to the nth columns by —4z, 4 ,(€), —42z, 40(),..., —4z,(6 respectively, 
and then add all to the last column. The elements of the last column then have 
the form 


py 
py ¥(x)z8) in the first row, 


and 


Y Unvde(6) — 


J 


J 


y Uvo(¥ )z/(8), v= 1,2,...,7, 
=U 


+1 


in the 2nd to (m + 1)th rows. By formulae (17) and (52) of $4 and formula (18) 
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above, these elements can be written in the following form, with the introduc- 
tion of P,,v = 0,1,...,7, 


1 1 
np" Po = np") y ePovis Yow], 
pp I y (1-8 k ky +1 . 
— = — eP®) ~ p “(—B wit + e 7 Pass ky OL wkytt ° 
np)” np NAY vor] Oe Pr tows] 


Moreover, we have the asymptotic formulae proved in §4.5 and 4.9: 


Vy — ePr"T1], y= 1,2,...,, 


p [aw 5”) for j —= 1,2,..., fb — 1), 
Uy) = ( p™([aywk] + e?[Byw}) for j = p, 
pve? Bow] for j =(% + 1), (u + 2),..., 
it n 
4A) =[]e™ Tl eb) + eth). 
v= S=ut 


We substitute these expressions into formula (34) of §3 in which the last column 
of the determinant H(x,€,,) has been changed in the above way, and distribute 
the factors of the denominator A()) as follows. We divide the (v + 1)th row by 
p*’ for v = l,...,n, the wth column by [6] + e°°+[6,], and the jth column by 
e°?i for j = (u + 1), (uw + 2),...,2. The formula then takes the following form 


pwo,x eP*T 1] pox <1 
ov Kt w'((aa) + e?°1Bi)) me Pp 
G(x.bd) = SO men (ol emma eel 
wi"((on} +e? “LB pl) ah 
| [~,0} ] ~ [65] + e?*[8,} [B,, n ] P,, 


By the conditions (13), all the exponents in this determinant have a real part 
less than or equal to zero; further, by the results of §4.9 (p. 69), the denominator 
(@.] + e?*+[6,] is bounded, in modulus, from below by one and the same 
number on all the arcs y;. Hence all the elements of the determinant, and 
therefore the determinant itself, are bounded on these arcs. 

Hence on the arcs y;, the inequality 


M 


= (20) 
|p|” 


|G(x,,A)| < 


holds, M being a certain constant. 
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We now show that the same inequality holds on the arcs y;. To do this we 
need only to multiply the elements of the first (uw —1) columns of the 
determinant A(x,€,A) by 42z,(€),4z.(€),...,42,-1(€), and those of the pth, 
(u + 1)th,....2th columns by —4z,(€),— Ee O.. ,—4z,(§), and add all to 
the last column. Then by repeating the previous argument, we find that the 
inequality (20) holds also on the arcs y,. 

Hence the inequality (20) has been shown to hold on the part », + y; of the 
arc y, lying in So. But the preceding arguments are clearly applicable to any 
domain S,. Hence it follows that (20) holds on the part of the arc y, lying in 
S, and indeed on the whole arc y,. Changing from p to A, we obtain the 
inequality (12), and the lemma is proved for n odd. 

(b). Let n be even; n = 2p. 

This case differs from the preceding one only in that, in two columns, the 

uth and (w + 1)th, the numerators 

eP“tl] — [é"] and e?™{1] — [€"] 
appear. These are bounded from below on the arcs »,, and therefore the 
inequality (12) holds on the arcs y,. Similarly, it may be shown that it holds on 
the arcs yx. 

By applying the lemma just proved, we obtain for the integral J, the estimate 


1 M M 
7 RR, ORR = RT» 


from which follows smmediavty the uniform convergence J, > 0 as k — oo. 
The considerations on p. 84 therefore lead to the following theorem. 


IZ;,| S57 


THEOREM 2 


The Green's function for a differential operator L generated by regular boundary- 
conditions can be expanded in a uniformly convergent series 


G(x,5) — __ y H%,8) 
yv=1 Xr, 
H,(x,€) being the residue of the function G(x,é,A) at the pole d,. 
If all eigenvalues of the operator L are simple zeros of the function A(A), and 
if the coefficients of the corresponding differential expression / satisfy the 
differentiability conditions necessary for the existence of the adjoint expression, 
then (see §3.8) 


(21) 


—fT WX, é) — y(x)z,(8), 


where y,(x),Z,(x) are eigenfunctions of the operators L and L* corresponding 
to the eigenvalues A, and 4, and are normalized so that? 


J, wCZAO)dx = 1. 


* We shall assume these normalization conditions to be satisfied in Theorems 3 and 4, 
which follow. 
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If we apply Theorem 2 to this case, we obtain 


THEOREM 3 


If all eigenvalues of an operator L of the form considered in §4.1, which is 
generated by regular boundary conditions and for which the adjoint differential 
operator exists, are simple zeros of the function A, then its Green’s function can 
be expanded in a uniformly convergent series 


(ne) = 5 *eve(6) (22) 


y 


From this theorem we easily obtain one on the representation of a given 
function. 


THEOREM 4 


Let L be an operator generated by regular boundary conditions with a differential 
expression of the form considered in §4.1 and for which the adjoint differential 
operator exists. Let all eigenvalues of L be simple zeros of the function A. 
Then any function f(x) in the domain of definition of the operator L can be 
expanded in a uniformly convergent series of the eigenfunctions: 


fo.6) 


fx) = dL »y0); 


v=1 


4 


= |) £2, 


where y,(x), z,(x) are the eigenfunctions corresponding to the eigenvalues 
A,, A, of the operators L, L*. 
Proof. We put Lf = h; then 


f(x) = [} GE n@)dé. 


In this we substitute for G(x,€) its ex pansion (22). Since the series is uniformly 
convergent, we may integrate termwise. Hence 


fix) = Y a0), 


y 


where 


toe _ i] _ lips 
a“, — A, I. h(é)z,()dé — A, (h,z,) ~~ A, (Lf,Z,) ~~ A, (f,L*z,) a A, (f,4,2,) 


— (f,2,) ° 
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The asymptotic formulae obtained in §4 also enable the convergence to be 
investigated of an expansion in terms of eigenfunctions of a given function 
which does not belong to the range of definition of the operator L. The results 
found are similar to those for ordinary Fourier series; see the papers by 
Birkhoff [5], Tamarkin [110], and Stone [109a]. 


Remark. It follows from the Corollary to Theorem 1 in §5.2 that every func- 
tion f belonging to L7(a,b), and not merely functions f belonging to the 
domain of the operator, can be expanded as a series in the eigenfunctions of a 
self-adjoint operator L and this series is convergent in the quadratic mean. 
Since the eigenfunctions of a self-adjoint operator are orthogonal to one 
another, the expansion also converges unconditionally. 

The question arises whether a similar proposition holds for a non-self- 
adjoint operator L. G. M. Keselman [49] and V. P. Mikhailov [75] have given 
an affirmative answer to this question on the assumption that the operator L 
is generated by regular boundary conditions, and that, for an operator of 
even order, 02 # 0_,0,. Their results contain sufficient conditions to ensure 
that the eigenfunctions of the operator L should form a Riesz basis in the 
space L?(a,b) — (see N. B. Bari [3]). We say that a complete system 


(Pn(X)n=1 C L7(a,b) 
is a Riesz basis if for all fe L*(a,b) 
Y | Preorseodx f < 00, 
and if, for any sequence of numbers ¢,, ¢c, ... such that 
Ylel< 00, 
there exists a function fe L*(a,b) such that 


b 
[, Fedex = ey, 2 = 1,2,3,... 


If {p,(x)} 1s a Riesz basis, then the following lower and upper estimates 
hold, instead of Parseval’s relation: 


m [7 LM) Pale <¥ | [? Fodescaddx <a [f(a Pa 


where the numbers m and M,M >m> 0, are independent of f. 
We remark that operators whose eigenfunctions and associated functions 
form a Riesz basis are spectral operators: see N. Dunford’s survey [21]. 
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4. Expansion in terms of Eigenfunctions when the Boundary Conditions are 
Separable 


The boundary conditions are said to be separable if each separate U,(y) con- 
tains the values of the function and its derivatives at only one of the two 


boundary points of the interval. Separable boundary conditions therefore have 
the form 


Uy) = » oY” = Oforj <m, 


1-1 (23) 
U,(y) = », B yy? = 0 for j => IM, 


For m = 0 or m1 = n, they are equivalent to the Cauchy initial conditions 
y® = 0 or yi) =0, k = 0,1,...,(1 — 1). 


It may readily be seen that, if0 < m <n, the separable boundary conditions 
are not regular; the only exception is the case where vis even and m = 4n. In 
this case the boundary conditions are of the Sturm type and are therefore, 
under certain conditions, regular (see §4.8). 

In the case of separable boundary conditions, the concept of the weight of 
a matrix is useful.! Let A be a matrix of rank m: 


Mo ses SH n-1 


nO eee Omn-1_Jj 


then at least one of the determinants 


(24) 


Om di mn ym 


does not vanish. We call the sum j, + ... +/,, the weight of the determinant 
(24), and the weight of the matrix A is the greatest weight of its non-vanishing 
determinants (24). 


LEMMA 2 
Among the determinants (24) of the matrix A there is just one non-vanishing 
determinant whose weight is equal to the weight of the matrix A. 


} This concept and also the proof given later in this section are due to M. V. Keldysh [47]. 


92 §5 EXPANSION IN EIGENFUNCTIONS 


Proof. We prove this assertion by induction. For m= 1, the assertion 1s 
clearly true. . 

Suppose it holds when the number of rows in the matrix A is less than m; we 
shall show that it then holds for a matrix of m rows. 

Let the number p be chosen so that the (p + 1)th column contains at least 
one non-zero element, while all the later columns contain only zeros, 1.e. 


Clearly, p + 1 >m. 
We consider the determinant 


Oy p—-mt1 1p 


Xm, p—-m+1 vee mp 

if it is not equal to zero, the proposition obviously holds. 

Suppose, however, that it vanishes, and let g be the rank of its matrix, so that 
q < m; then (m — q) of its rows are linear combinations of its remaining rows. 
Without loss of generality we may assume that these are its last (#2 — q) rows, 
since an interchange of rows in the matrix A merely changes the signs of the 
determinants (24). We add to the last (m7 — qg) rows of matrix A certain linear 
combinations of its other rows (this does not alter the determinants (24)) to 
bring A into the form 


4 4 , 4 
19 ceo O1 p—m 1 p-mt1 eee Oi 0 eee ) 
, 4 / 
A! “40 q,p—m Og,p—-mt1 gp 0 0 
/ 
%7+1,0 %o+1,p—m 0 0 0 0 
, ¢ 
%m0 m,p—m 0 0 0 v 
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is obviously equal to g, and the rank of the matrix 


/ / 
Wo+1,0 +e. Co+1,p—m 


/ 
moO Oy, »p—m 


1s equal to (#7 — q). By the hypothesis of the induction, the matrix A; has just 
one non-zero determinant of the maximum weight; let this be 


/ t 
1 im-gti “ee 7 ina 


/ / 
OG, im-qtr “ee Caim 


Similarly, the matrix A, has just one non-zero determinant of maximum 
weight, say 


t / 
Wgtli, +t WG4tljimag 
D’ = 
q 
oa ou! 
m,Ji m,Jm=q 
Then the determinant 
t , t a! 
Ljy ljm-q lJIm-qt+1 I,Jm 
/ , a’ a’ . 
D Cod Xo,jm=q QiJm-qt1 4,Jm 
7 , r 
Cots Co+lim-aq 0 0 
, , 
Xm, Js Om im-4q 0 0 


is obviously non-zero, and all other non-vanishing determinants of the matrix 
A’ have a smaller weight. For, any such determinant A has at least one sub- 
determinant A, formed out of the first g rows of the matrix A and which, 
together with its complement A,, is non-zero. The weight of the determinant 
A is equal to the sum of the weights of A, and Aj. 

But by definition of D, and D,, we have 


weight of A, < weight of D,, (25) 
weight of Aj < weight of Dy, (26) 


and the equality signs are possible only if A, coincides with D, or Aj with Dy. 


Hence 
weight of A < weight of D, 
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and the equality sign appears only if it occurs in each of (25) and (26). This is 
clearly only possible if A coincides with D. 

Hence Lemma 2. 

We now consider the matrices 


X10 1 n-1 Bn +1,0 Bmtin-1 


A=|. . and B= 
mn-1 Bro see Ban-1 


formed from the coefficients of the boundary conditions (23). 

Since these boundary conditions are linearly independent, the rank of each 
matrix A and B is equal to the number of its rows; hence we may speak of the 
weights of these matrices. In future, we shall denote the weight of the matrix A 
by A and that of the matrix B by /. Further, we shall put a = 0 and b = 1; as 
before, this causes no loss of generality. 


LEMMA 3 (M. V. Keldysh) 


An operator L which is generated by the differential expression Ky) = y™ 
and the separable boundary conditions 


n 
», % no? =O forj <m, 
k=0 


n-1 
> Buy? = for j > m™, 
kK=0 
which are not equivalent to the Cauchy initial conditions has infinitely many 


eigenvalues 4, 
For these eigenvalues the following asymptotic formulae hold: 


A, = Py (27) 
where 
ye n h+l 
_. 2 n 10 ] 
Py Jum y > (28 a) 
sin —— 
n 


ifn — im = 2« and if v is sufficiently great, and 


yh _ hd 
py = menitin, " +O ! 28 b 
sin (2x + 1)z y ( ) 
n 


ifn —m = 2« + 1 and if v is sufficiently great. 
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Proof. We put A = p"; the equation y™ = p"y has the general integral 


y= > cero, 


jel 


where w,, j = 1,2,...,1, denote the different nth roots of unity. By substituting 
this expression in the boundary conditions (23) we obtain a system of homo- 
geneous equations 


-1 
» ¢q > Bjwgp'e?? =0, j=(m+1), (m+ 2),...,n 


for the unknowns ¢,,C»,...,c,. Hence the eigenvalues for the operator L in 
question are determined by the equation 


n-1 ke n—-1 kk 
>» 04. 4p vee >» 1 ,VnP 
k=0 k=0 
n-1 n-1 kk 
mk p™ cee > mk?OnP 
=0 = 
Doe) =|. 7 =0 29) 
y Bin+1,nip eo? vee » Bra+1,k@nP eo” 
k =0 k =0 
n-1 n-1 
»y Brxwoip eo? ves » Brrconprerr? 
| K=0 k—0 


Clearly, we need discuss only the case 1 <m < 4n, since the case 4n << m 
<n — 1 can be reduced to the first case by the substitution x, = 1 — x. 


We put 


A;; = .y Op EPs i = 1,2,...,m; j = 1,2,...,n 
- (30) 
B= 2, Biwip’, i=(m+1),...,n; j= 1,2,...,n. 
By Laplace’s expansion theorem, 
Axj;, «+ Atjm | | Bmttimer ot Bmttin 
A(p) =D] - J dale . Pim ts Fi? (31) 


Am), cee A inj Bn, jm i Bin 
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the summation being over all possible permutations J,, jo,---sJmsJm+12Jm +29°*+9Jn 
of the numbers 1,2,...,7 satisfying j; <j, <<... <j, and j,j4, <<... <J,. In 
other words, D(p) has the form ‘ 


N 
D(p) = i P,(p)e**, (32) 


where P,(p) is a polynomial in p of degree < n(n — 1), and all the exponents 
p, are different. 

We are concerned with those terms of this sum in which the real part of the 
expression w,; , +... +, takes the greatest value. We shall shortly see 
that for each such term the polynomial P,(p) associated with the term consists 
of just a single product of two determinants from (31). Hence it suffices to 
show in each case that neither of the determinants 


| 
Aj, “ee Aj;,, Buti, jm+ Butlin 
; . . (32 a) 
A mj, Anim Brim +1 Bin 


in this product vanishes identically. We shall prove this assertion for the first 
of the determinants; the proof for the second of the determinants is similar. 
We have: 


hy h 
Aj; eee Ay,, Oth, eee Cn, Wj eee w 
_ y 4 oo. pir tee thm (33) 
hy shes.eestm 
h h 
| Anj, vee A nim nh, » Cnhy, || Oi, a Din 


But if Aw, , +... +,,) takes the greatest value, then the points 


JIm+i 
Oj 4p jme yee Mz, On the unit circle are situated as far to the right as they can 
be in any of their possible arrangements, and consequently, when numbered 
appropriately, they lie within an angular distance of 27/n from one another 
(Fig. 14). Ifthe numbers w,,w»,...,@, are numbered cyclically in a suitable way, 
then 


— fk-1 — 
oO; = C Wj > k= 1,2,...,7, 


where £ = e?"", Introducing the notation £, = ¢*, we have 


hy hm | 
Wi ooo Wi, 1 l 
Ch, cee Chim 
— — hy+..4+h 
a=|. bee = wptetim #0, (34) 
h, h m-1 m-1 
Wi, eee w i” | ; Ch, ooo Ch. 


since the numbers ¢,,, ¢,.,...6,,. are all different. 
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WwW. 
JIm+2 


JIm+t1 


Fig. 14. 


Hence we have only now to show that among the various non-zero deter- 
minants 


On, bn 


(35) 


Xmh, Omhin 


there is just one for which the exponent /, +... +h,, of p takes the greatest 
value. But this exponent is equal to the weight of the determinant (35); hence 
our assertion coincides with the enunciation of Lemma 2. Hence the first, and 
similarly the second, determinant does not vanish. On the other hand, the 
values of w, . +... + ,, for the separate combinations with greatest real 
part are obviously different. Hence in the corresponding polynomials P, only 
one single product of determinants appears, and it has been proved that 
the terms which have the greatest value of Rw, ,. + ... + w;,) cannot vanish 
identically. For (n — m) even, there are just two such combinations (cf. 
Fig. 15); for (n — m) odd, just one. In the latter case, we notice that the two 
combinations having a not smaller real part of w; , +... + ;, also have 
non-vanishing coefficients in formula (32), since these too consist of only a 
single product of determinants. 
We examine the cases (n — nz) even and (n — m) odd separately. 


1. (n — m) even, n — m = 2k. 
By numbering the numbers w, suitably we obtain two sums 


W — + W_ K+] ove + We-1> W_K+] + We Ke+2 + vee + Wx 
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We 
Wy 
Wy 
Wo Wo 
WH] Woy 
W—2 
Fig. 15. 


with the greatest real part (see Fig. 15, drawn for « = 2). We put arg p = 
and consider the function D(p) in the sector S: — <op< - In this sector 
at least one of the exponents 

(wig + Wig ty Hee F Wy —a)P, (M41 + Onde +. + Ox)p 


has areal part greater than the real part of all other exponents in the expression 
for D(p). Let h be the weight of the matrix 


X10 O1n—-1 
A= , 
Omn0 Om,n-1 
and / the weight of the matrix 
Bin +1,0 eee Bm+i.n-1 
B= 
Bo see Buin-1 
Then by the above investigation, in the sector S we have 


D(p) _ goons tocer tate’ gh HY 4- o(*)] 
p 


+ goers sa-nistonsedogi ght 4 o(*)]. 
p 
Now «’ and §’ clearly differ only by a factor of (—1)"~ te?" + Dis, 


B’ — y’. (—1)"~4e22#@ + Din 


because we obtain the coefficient 8’ from the coefficient «’ by replacing the 
order j,,/o,--.,J, In the determinants (32 a) by the order js, js,...,/,,/,;. Hence the 
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sign of the product of determinants in formula (31) changes by (—1)"7}, 
while each of the numbers w j.18 multiplied by = e®™", so that, for example, 
a in formula (34) is multiplied by e?*/", Hence we have 


D(p) = apt tel 8 Omer ton tone | + (5) 
p 


+ (Ay Fetath Delon 00 + (7) |} 
p 


We put w, =o + i7; then w_, =, =o — ir, where 


7 = sin —. (36) 
Hence we have 


D(p) _ a ph *tglo-mt oo-e + o(*) 
p 


+ (—1) —1,2itp + 2ni(h+1)/n 1 + o(7) |} (37) 
so that the equation D(p) = 0 is equivalent to the equation 
(—1)"-? + o(*) + e2'lte +7 (h+D] — (. 
p 
Hence 


e2ilte +7 (h +1] _ (—1)" + o(*) 
p 
and so for the required eigenvalues we have 


2i| tp + (h +)) = (2 + mri + 0(4), vy = 0,-+1,+2,.... 
Therefore 
‘ +5 
pe EE oft 
7 p 
and so 
v5 a 
p= a——— o(*). 
T V 


On the other hand, by applying Rouché’s theorem, it may easily be shown 
that (cf. the detailed investigation in §4.8 for the case of regular boundary 
conditions), for sufficiently large |», there is in the neighbourhood of the point 
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just one zero p, of the determinant D(p), and 


y 42 _ Atl 
2 n ~ 1\, 
ar I, of! 
i.e. by (36), 
ype hs! 
— 7 2 he o(*) 
Py . Dear Vv 
sin —— 


2. (n — m) odd;n —m =2« +1. 


2 , 
In the sector S: — — < go < 0 the expressions 
n 


(we fois fo + oy)p and (p41 + Octo $+ + x42) 
have a real part greater than any other expression 
(0 jpeg Hove #0 y,)P 
By considerations similar to those of case 1 we find 


D(p) — elon tenner tteay | 4 o(*)] 
Pp. 


+ (AI teen tecnesta teen dng ET + o(*) |. 
_ \p 


Hence it follows that D(p) has infinitely many zeros p, in the sector S, where 


in = 2 ne 1 
Py — 7 n “gp et] + o(*). 


This concludes the proof of Lemma 3. 


We now consider briefly the question of expansion in terms of eigen- 
functions and associated functions. 

An expansion in terms of the eigenfunctions of a differential operator does 
not, in the general case, have the usual orthogonal properties. However, as 
we proved in §5.3, this expansion does have properties similar to those of an 
ordinary Fourier series expansion if the operator is defined by regular bound- 
ary conditions. If the boundary conditions are separable, the situation is 
more complicated. 
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W. Ward, [117], investigated the eigenfunction expansion for the boundary- 
value problem 


y? +rAvy = 0, pO) =... = y""0) = y(1) = 0. (38) 


He proved that an expansion in a uniformly convergent series in the eigen- 
functions is possible only for functions which satisfy certain analyticity condi- 
tions, and consequently such an expansion is far from being always possible 
for all functions in the domain of definition of the corresponding operator L. 

On these questions of eigenfunction expansions there is quite an extensive 
literature (see [17], [95], [122], [S0a, b], [43], [44], [102], [116]). 

Here we shall only go Into some results due to A. P. Khromov [50c] and 
W. Eberhard [24], who make use of Lemma 3 (M. V. Keldysh) in their 
investigations. 

A. P. Khromov considered the operator L corresponding to the differential 
expression 


Ly] = y + polxyy? +... + p(y, py e C0, 1) (39) 


and general separable boundary conditions. It is not difficult to see that in 
this case the assertion of Lemma 3 remains true. In what follows we need the 
concept, introduced by M. K. Fage [26], of an /-analytical function. We mean 
by this a function f to which the differential operator / can be applied arbi- 
trarily often and which satisfies inequalities of the form 


(fy tte |< cmon +9)? 


p=0,1,2,.., ¢=O0,1,...,(n — 1), 


(with the exclusion, possibly, of p = 0 and q = 0) in an arbitrary interval 
[x,8] < (0,1); here C is a constant which depends only on /, f, « and f. It 
turns out that the properties of an /-analytical function are in many respects 
similar to those of an ordinary analytical function. 

Let {y,(x)}f_-, be the sequence of all the eigenfunctions and associated 
functions of the operator L under consideration, numbered in order of increas- 
ing absolute values of the eigenvalues. A. P. Khromov has established the 


following proposition. 


Theorem 5. 


ce 
If the series ) ayp,{x) converges uniformly on some interval [Xo,x,] | (0,1), 
then: e=t 
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(a) the several series 


~ d 
2M (5) Pea p=0,1,2...,  g=0,1,.-.(2—1) 


converge absolutely and uniformly on every interval [0,a] < (0,x;). 
(b) the sum f(x) of the series 

2 PAO) 
is an l-analytical function; 


(c) the functions 


fx), WfO)L PLO)... 


satisfy the boundary conditions at x = 0. 

Conversely, if a function f(x) which is summable on the interval (0,1) is an 
l-analytical function in the interval [0,a), 0 <a< 1, and if the functions 
T(x), fo], PL/Od),... satisfy the boundary conditions at x = 0, then f(x) 
can be expanded in the system {p,(x)};=1 as a series which converges uniformly 
in every interval [0,b] < [0,a)  [0,R), where the number R is defined by 


l im | d ee 1/k 
aa~imla l(a) ee 
[x] being the integral part of the number x. 

W. Eberhard’s results ({28] —[31]) have much in common with those of 
Khromov. Eberhard considered an operator L corresponding to the differ- 
ential expression (39) with coefficients p,(x) which are holomorphic in the 
circle |x| < 1 in the complex x-plane, and to general separable boundary- 


oO 
conditions. He proved that, if the series )° a,@,(x) converges uniformly on 
k=1 


the interval [x,,x,], then this series converges uniformly in every closed sub- 
domain in an open polygon (and consequently that the sum of the series is 
regular in this polygon) whose vertices are at the points 


COS 2na/n p(y linn 
° 9 


x}. 
cos m/n 


in the case where n — m = 2, 
and at the points 
cos (2x + 1)a/n “avian 


vy = ],...,11 
cos 7/n oe 


1° 


inthe case where n — m = 2x + 1; 


here, as usual, mm is the number of boundary conditions at x = 0. Eberhard 
has also obtained certain necessary conditions, in terms of ordinary analy- 
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ticity, for the function f(x) to admit a uniformly convergent expansion as a 
series in the system {p,(x)}°_,; and here certain additional restrictions, as well] 
as holomorphism, are imposed on the coefficients ofthe differential expression 
(39). 


5. The Case of Multiple Poles of the Green’s Function; m-fold Completeness; 
Keldysh’s Theorem 


We again consider, as we did earlier, the generalized eigenvalue problem 


Ky) =f, (60) 
U,(y) = 0, v = 1,2,...,n; (61) 


let Aj,A2,Ag,... be all the eigenvalues of this problem. They may in general be 
multiple zeros of the characteristic determinant A(A) and consequently 
multiple poles of the Green’s function. Let A be any one of these eigenvalues, 
and let the eigenfunction and its associated functions which occur in the 
canonical system belonging to A be 


VoVise++9Vq-1: 


We construct m systems of functions 


d’ t t 
(v,k) A k — —_ 
y = Peraer 1 Vk-17y Toe +94) | , v =0,1,2,...,0m — 1). 


Each of these systems 
yO yo yor -1) 


a ] 


will be called a derived system from the original system y,y4,...,¥q—1- 

The system of all eigenfunctions and associated functions of the boundary- 
value problem (60), (61) is said to be m-fold complete if every system of m 
functions 


(fi) fol). Fm OO}, 


each of which has continuous derivatives up to order 7 inclusive and satisfies 
the boundary conditions (61), can be represented as the limits of m uniformly 
convergent sequences of finite linear combinations 

fx) = lim Y Cy 


N—-o i,k 


whose coefficients are independent of v. 

The considerations in §5.3, §5.4, together with certain refined results in 
the theory of analytic functions, and a formula for the principal part of the 
Green’s function (see §3.9), enable the following general theorem to be derived, 
which we merely state here without proof. 
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THEOREM 6. (M. V. Keldysh, [47]) 


Let 
iy) = y™ + pax, rA)y""? + >.. + Lp(a) + A” y, 


where p,(x,A) is a polynomial in A of degree < Km. and suppose the boundary 
n 


conditions may take the form 


» a, yt"? =0, j= 1,2,....); O<p <n, 
k=1 


Y Bvt =0,  §= (P+ VE + 2pm, 


or 
y? = yh, k= 0,1,2,...,(7 ~ I). 


Then the system of eigenfunctions and associated functions for this boundary- 
value problem is m-fold complete. 


In the article by V. N. Vizite1 and A. S. Markus [114] it 1s proved that 
Theorem 6 remains valid also in the case of regular boundary conditions, 
where, when 7 is even, it is required that 06 # 0_,0,. Further, it remains true 
in the case of general self-adjoint boundary-conditions; the conditions of 
periodic type, which feature in Theorem 6, are a particular case of these last. 

An expanded exposition of M. V. Keldysh’s results (in particular, those 
on which Theorem 6 is based) is given in the monograph by I. Ts. Gokhberg 
and M. G. Krein [33]. In this monograph have been collected the most 
recent results in the spectral theory of completely continuous non-selfadjoint 
operators. (We remark that the resolvent, which we have discussed above, of 
the operator L is a completely continuous operator.) 

In addition to the problem of eigenfunction expansions, a problem of great 
interest is that of determining a differential operator from the specification 
of its spectral characteristics, e.g., from the eigenvalues and norms of the 
eigenfunctions (with known initial values). Particular aspects of this problem 
are considered in V. M. Levitan’s book [63f] and in the articles by Z. L. 
Leibenson [59a, b]. 

Surveys of the general state of the theory of linear non-selfadjoint opera- 
tors, both differential and those specified abstractly, are given in the papers 
[21], [81a], [48], [73]. 

The articles [52], [72a], [100], [23b], [59a] may also be consulted on the 
subjects discussed in this chapter. 


CHAPTER III 


DIFFERENTIAL OPERATORS IN A SPACE 
OF VECTOR-FUNCTIONS 


86. Basic Concepts 


1. Linear Differential Expressions in a Space of Vector-Functions 


Let R,, denote an 7-dimensional complex vector space; i.e. R,, consists of all 
vectors 


VY = Ovo Yn)» 
where each y, is a complex number. 
Functions y = y(x), of the real independent variable x, whose values are not 
numbers but vectors in R,, are called vector-functions. A vector-function is 
therefore simply a system of #7 complex-valued functions 


YX) = 04), ¥2(%),- Pm) 5 
each of the scalar-functions y,(x) 1s called a component of the vector-function 
y(x). 
The function y(x) is said to be continuous at the point x, if all its components 
are continuous at Xp. 
Similarly, a function }(x) is said to be differentiable if each of its com- 
ponents is differentiable, and by definition 


y'(x) — (1(X),¥2(%),- ; Vm(X))- 
Derivatives of higher order are defined in a similar way. 
It may easily be seen that: 
(Vtz)y=y tz, Ay HV ty, 
(y,2)' = 0.2) + 2) 
where 


(9.2) = vase). 


In addition to vector-functions we shall also be concerned with operator- 
functions. The values of operator-functions are linear’ operators in R,,. Such 


1 The reader is reminded that the author’s definition of linearity, on p. 2, includes 
homogeneity. 
105 
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operators can be represented by means of square matrices A(x) = [a;,(x)] of 
order m, whose elements are scalar functions. We shall not, from now on, make 
any distinction between operators and the matrices which represent them. 

Essentially, operator-functions are also vector-functions, since the aggregate 
of all linear operators is a vector-space of dimension m?, Consequently an 
operator-function A(x) will be said ‘to be continuous at the point xz, if all its 
functions a,,(x) are continuous at Xo, and to be differentiable at x, if all the 
functions a,,(x) are differentiable at x9. So, by definition, A’(x) is the matrix 
whose elements are a;,(x). We see without more ado that the following rules 
hold: 


(A+B) =A'+B’, (AA) =A HAA, 
(AB)’ = A'B+ AB’, (Ay)’ = A’'y + Ay’. 


We denote by C™ the aggregate of all vector-functions y(x) which have 
continuous derivatives up to the mth order inclusive in a fixed interval [a,5]. 
Let Po(x),P(x),...,P,(x) be operator-functions which are continuous in [a,b], 
and suppose det P,(x) # 01in [a,b]. An expression of the form! 


I(y) = Po(x)y + PydyO"P +. 4 PCy (1) 


is called a linear differential expression in the space of vector-functions. 

It is defined for an arbitrary vector-function y(x) belonging to C”, and 
itself represents a vector-function, continuous in the interval [a,b]. The theory 
developed below for such differential expressions and the differential operators 
corresponding to them is formally similar to the theory already presented for 
scalar functions. 

We remark that, essentially, /(y) is a system of m differential expressions of 
the nth order which depend on m scalar functions y,(x),¥2(x),...,)',(x). For 
n = 1, we consider, for example, the expression 


Ky) = y’ + P(x), 


and resolve the vector-function z = /(y) into its components by putting 


I(y) — (1,(y),lo(),. , otal): 


Then we have 


m 


ly) — yj +) Pi QOYns J — 1,2,...,277, 
where P,(x) are the elements of the matrix P(x). 


1 We remind the reader that, e.g., the expression A(x)y'"~»)) denotes the result of applying 
the operator A(x) to the vector y'"~ (x). 
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2. Boundary Conditions 


We denote by 
VarVarVa 3s VoVosruYo > 


the value of the vector-function and its first (1 — 1) derivatives at the points a 
"-1) are vectors in the space R,,. We put 


and 6 respectively, so that y,,...,y} 
U(y) = Apa + Av +... +A, yO" 
+ Boy, + Buryp +... +B, ”, (2) 


where Ag,...,A,-1,Bo,..--B,-1 are fixed linear operators in the space R,,. If 
several such forms are given, U,(y),U,(y),...,U,(), then equations of the form 


U(y) = 0, U2(y) = 0,...,U,(y) = 0 (3) 
are called boundary conditions. 

We denote by Z the aggregate of all functions ye C™ which satisfy the 
conditions (3). Let L be an operator, with J as its domain of definition, which 
is defined by the equation 

Ly = Ky); 
then L is called the differential operator which is generated by the differential 
expression I(y) and the boundary conditions (3). 

We put 
Uy) — Ay oa + cee + Ay ne + By ob + see + By n-ws”, 

y = 1,2,...,q. 
We will assume in the definition of Z that the forms U,(y) are linearly in- 
dependent; this implies that the rank of the matrix formed from all the elements 
of the matrices [A,,], [B,,], viz. 


Ajo). . 541 ,-1)Brio;- , By n-1 
A49;. , Ao n-1Bo95- , »»Bon-1 


| A joy--+sAgn-1»Bgo-+»Ban—1 


is equal to mq; for each form U,(y) has m components. 
From now on we shall mainly be concerned with the case g =n. 


3. The Homogeneous Operator Equation 


We want now to examine the following homogeneous equation in an unknown 
operator function Y(x): 


WY) = P(x) Y™ + Py(x) YC D+... + P(x) Y =0. (4) 
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As before, it is assumed that P,(x),...,P,,(%) are continuous and that det P, # 0 
in the fixed interval [a,5]. 

We mention some of the most elementary properties of solutions of this 
equation; they resemble the properties of sdlutions of ordinary, homogeneous, 
nth-order equations. 

Solutions Y,,Yo,...,Y, of (4) are said to be linearly independent if the 
equation 

YC, + ¥,C, +... + Y,C, = 09, 


where C},C,,...,C, are constant operators, holds only for C, = C, =... 
= C,, = 0; otherwise they are linearly dependent. 

A system of n solutions Y,,Y,..., Y, is linearly independent if and only if the 
determinant of the n® X n?® matrix 


Y, Y, Y,, 
Y; Y, Y, 
L ye-» ye-» . ye-v | 


formed from the matrices Y is non-zero. 
If the solutions Y,, Y,,...,Y, are linearly independent, then any other solution 
Y has the form 


Y — YC, + YC, + eee +- Y,C,: (5) 


Consequently, formula (5) defines the general integral of equation (4). 

If Y is a solution of equation (4), and if c is a fixed, constant vector, then we 
get a solution of the equation y) = 0 on putting )(x) = Y(x)c. To see this, we 
need only apply both sides of equation (4) to the vector c. 

Hence, and from formula (5), we can at once infer: 

Any solution of the equation I(y) = 0 has the form 


l= YiC) + Y Co + eee + Y Crs (6) 


Where Cy,Co,...,C, are arbitrary constant vectors in R,,. 


4. The Homogeneous Boundary-Value Problem 


The problem of determining a vector-function y which shall satisfy the con- 
ditions 

Ky) =0, (7) 

Uy) =90, v= 1,2,...,n, (8) 


is called the homogeneous boundary-value problem. 
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We consider the n? x n?2 matrix 


U(¥,) UY.) ... U,(Y¥,) 
UY) U¥2) ... Ua Y;) 

U= (9) 
U,,( Y)) U,( Y.) U,,( Y,) 


The results in §6.3 imply: 


I. A homogeneous boundary-value problem (7), (8) has a non-trivial solution if 
and only if the determinant of the matrix U vanishes. 


5. Lagrange’s Formula; the Adjoint Differential Expression 


In addition to the assumptions made in §6.1, we now further require that, 
for k = 0,1,2,...,n, the coefficient matrices P,(x) shall each be continuously 
differentiable (n — k) times. We denote the scalar product of the vectors 
y, ZE R,, by (y,z). Integrating by parts, we obtain 


[. ((y),z)dx = P(,0) + [. (y,1*(z))dx, (10) 
where P(7,{¢) is a bilinear form in 


N= VarVar-- Ve? Yu Voss Ve Y) 
and 
C = (2,,24,...,26" -,2,,25,...,20 9), 
and where! 
I*(z) = (—1)(P#2) + (— 1)" UP #2) -Y oe Pz. (11) 


The differential expression /*(z) is said to be adjoint to [(z). Formula (10) will 
be referred to as Lagrange’s formula. A differential expression /(y) is said to be 
self-adjoint if [*(y) = I(y). 

Repeating the argument in §1.5, the following general representation for a 
self-adjoint differential expression is obtained: 

II. Any self-adjoint differential expression is a sum of Self-adjoint differential 
expressions of the form 


lay) = (Py), lay 1) = FLGP YO PY + PYM YOY, 


where P = P(x) denotes an operator function all of whose values are Hermitian 
matrices in R,,. 


11f A = [a,,] is a square matrix of the nth order, we denote by A* the adjoint matrix, i.e. 
the matrix B = [b,,] for which b,, = a,,. A matrix is called Hermitian if A* = A. 
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6. Adjoint Boundary Conditions; Adjoint Operators 


We supplement any given set of linearly independent forms U,,...,U,, to form a 
complete system of linearly independent forms U,,U,,...,U2,. Repeating the 
argument of §1.6, we can transform Lagrange’s formula to 


f, (U(y),z)dx — (U,, Von) + (U2, V2,-1) Toe (U2,,V 1) + I. (y,P*(z))dx, (12) 


where V,,V,...,V2, are linearly independent forms in the variables 


, (n—- 1) t (n-1 
ZarsZigevey Ze 7 zh... Ze DY. 


The boundary conditions 
V, = 0, y = 1,2,...,n, (13) 


(or any conditions equivalent to them) are said to be adjoint to the boundary 


conditions 
U, = 0, v = 1,2,...,n. (14) 


The operator generated by the differential expression /*(y) and the boundary 
conditions (13) is said to be adjoint to the operator L generated by the 
differential expression /(y) and the boundary conditions (14). It will be denoted 


by L*, 
It follows from formula (12) that, for the operators L and L*, the equation 
b b 
[, (Ly,z)dx = [, (y,L*z)dx (15) 
holds. 


Conversely, if a differential operator L, satisfies the condition 


b 
[ Laz) = | (Liz)dx 
for all functions y,z in the domains of definition of L and L, respectively, then 


L, = L*, 
For vector-functions (x) we define the scalar product <y,z> by 


b 
(v2) = [ (v@),2())dv. 
Then formula (15) can also be written as 
(Ly,z> = (),L*z). (16) 


An operator L is self-adjoint if L* = L. In other words, an operator L is 
self-adjoint if it is generated by a self-adjoint differential expression and self- 
adjoint boundary conditions. For a self-adjoint operator L, equation (16) 
takes the form 

<Ly,z> = <y,Lz). 
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The homogeneous boundary-value problem L*z = Ois said to be adjoint to 
the homogeneous boundary-value Ly = 0. By an argument essentially the 
same as that in §1.7 we obtain: 


The rank of a homogeneous boundary-value problem is the same as that of the 
adjoint homogeneous boundary-value problem. 

In particular, 
A homogeneous boundary-value problem has a non-trivial solution if and only if 
the adjoint homogeneous boundary-value problem has a non-trivial solution. 


7. Eigenvalues and Eigenfunctions of a Differential Operator 


A number J is called an eigenvalue of an operator Lif there is in the domain of 
definition of La vector-function y # 0 such that 
Ly = dy. 


All the basic theorems in §2.1, 2.4, 2.5 remain true for operators in the space 
of vector-functions. In particular, the eigenvalues are the zeros of the 
characteristic determinant 


U(Y,) ... U¥n) 
A(a) = , . , (17) 
U,,( Y;) U,( Y,,) 
where Yj, Yo,..., Y, are linearly independent solutions of the operator equation 


KY) —aAY =0. 

These solutions can be chosen to be integral, analytic, operator functions of 
the parameter A; hence A(A) also is an integral, analytic function, and Theorem 
I of §2.1 carries over to the present case. 

The argument in §2.2 also applies here, so far as the generalized eigenvalue 
problem is concerned. 

We may observe that the ordinary eigenvalue problem in the space of vector- 
functions is equivalent to a certain generalized eigenvalue problem for scalar 
functions. To illustrate this, we will examine the eigenvalue problem 


d 
os + A(x)y = dy, O<x <l, (18) 


y(0) = yQ), (19) 


where y is a vector-function in the 2-dimensional space R, = yx) = {1(%), 
yo(x)}, and A(x) is the operator defined in this space by the matrix 


AC) 0 1 
OT aoe 0 | 
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Consequently 
A(x)y = {yolx), Px) yi(%)} 


and written in terms of coordinates the problem (18), (19) takes the form 


Go + ya) = AW), 
(20) 


as + p(x)y,(x) = Ay,(x), 


yi) = y,(1); y2(0) = »,(1). (21) 
By differentiating the first equation in (20) to obtain 


d*y,  dyr dy, 
dx? + dx A dx’ 


and then using this equation with (20) to eliminate y, and a we obtain for y, 


the differential equation 


d? d 
=o — 2A — + [A? — p(x)]y, = 0. (22) 


By the first equation in (20) again, the boundary conditions may be written in 
the form: 


y,(0) = y, (1), 


dy, _[a%_, (23) 
dx * x=0 ax M4 x=1 


The original problem therefore reduces to the following generalized eigen- 
value problem: 


d? d 
3 2A — + [A? — p(x)]y, = 0, 


dy, 
— ‘i ~ anf 


y,(0) = yi); be — Ay, |. 


This remark often enables a generalized eigenvalue problem in the space of 
scalar or vector-functions to be reduced to a certain ordinary eigenvalue 
problem in the space of vector-functions. 

In connection with the generalized eigenvalue problem in the space of 
vector-functions, see also [95a,b, 122). 
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8. Differential Operators of the First Order 


The eigenvalues and eigenfunctions are determined most easily in the case of 
a differential operator L of the first order. The operator is generated by the 
differential expression 


Ky) =y' + PO) (24) 
and the boundary conditions 


Ay, + By, = (Q, (25) 


We want to find the eigenvalues and eigenfunctions of the operator L. To do 
this, we first determine the solution of the equation 


vi + P(x)y = Ay. (26) 


. . . . \ 
Putting y = e**z in (26), we see that z has to satisfy the equation 


z’ + P(x)z = 0. (27) 


To each vector c € R,, we assign that solution z(x) of (27) which satisfies the 
initial condition [z],_, = c. Since the solution is unique, the correspondence 


Cc > 2(x) 


is one-to-one. Moreover, it is obviously linear: c, > z,(x) and c, > z,(x) 
imply Ac, + uc, > Az, (x) + wz,(x) for arbitrary numbers A,u. Hence for each 
fixed x, this correspondence defines a linear operator Q(x): 


z(x) = Ox)ec; 
and, by the definition of Q(x), O(a)c =. 
Hence the general solution of equation (26) is of the form 
y = e*Q(x)e, (28) 


where c is an arbitrary vector in R,,. 
Substituting this expression in the boundary conditions (25) and removing 
a factor e“7, we obtain 


[A + &?-2BO(b)c = 0. (29) 
Next we consider the generalized eigenvalue problem 
[A + pBQ(d)|c = 0 (30) 


for the operators A and BQ(d) in the finite-dimensional space R,,. Let 
Lz sf4os+++9éy be all the eigenvalues, and cj,Cg,...,c, the corresponding linearly 
independent eigenvectors for this problem. Equation (29) shows that e*?~ is 
an eigenvalue and c the corresponding eigenvector of problem (30). 
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Hence all the eigenvalues of the operator L are of the form 


l 
A = bh—a In p; + 2kai, J = 1,2,...,7:K = 0,+1,+2,+3,...; 


and the corresponding eigenfunctions are 
l 
jexo(, 4 Ing; + 2k) | O(x)cj. 


Many authors (see, e.g. [6]) have investigated the more general eigenvalue- 
problem 


d 
ay, + Plxly = O(a), (31) 


Ay, + By, = 9, (32) 


With operator functions P(x), Q(x) continuous in [a,5]. 

The corresponding results are formulated most simply for the case where, 
for each x in the interval [a,5], all the roots of the characteristic equation of the 
operator Q(x) are different, say, 1,t49,...5/4,- In this case, they can be numbered 
so that for each x the functions p,(x),p.(X),...,4,(%) are continuous. Hence 
there exists a non-singular matrix S (i.e. one for which det S(x) # 0) such that 


S~ (x) QO(x)S(x) = M(x), 


where M(x) is a diagonal matrix having p,(x),19(X),...,4,(x) as its diagonal 
elements. 

We assume that Q(x) is a differentiable function in the interval [a,b]. Then 
S(x) and M(x) are also differentiable. Using the substitution y = S(x)), we 
can reduce equation (31) to the form 

& + Poop = AMOS, 


where 
P(x) = S7*(x)S’(x) + S7(x)P(x)S(x). 


Henceforward, therefore, Q(x) may be taken to be a diagonal matrix. It 
turns out that for large |A| the solution of (31) behaves approximately like that 
of 

dy 
de MODY: 


More precisely: there is a fundamental matrix Y(x) of solutions of equation 
(31) for which 


roy =[ 109+ of] (yf 


where Y,(x) is a definite matrix independent of A. 
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This result enables asymptotic formulae for the eigenvalues and eigen- 
functions to be obtained, and also a theorem to be derived on expansions in 
terms of eigenfunctions of a boundary-value problem, exactly similar to the 
expansions in §$4 and 5. For more details, see [6]. 


§7. Green’s Function for a Differential Operator 


1. Inversion of a Differential Operator 


Let L be a differential operator for which the boundary-value problem Ly = 0 
has merely the trivial solution. Then L has an inverse operator L~!. We shall 
show that L~! 1s an integral operator, the kernel of which is represented by an 
operator function in two variables. This function is called the Green’s function 
for the operator. 

The precise definition of Green’s function for the present case is similar to 
that given in §3.3 and now runs: 

An operator function G(x,é) in the variables x,€ is called the Green’s function 
for the operator L if it satisfies the conditions: 


1°. Each value of the function G(x,§) is a linear operator in the space R,,. 

2°. The function G(x,€) is continuous and has continuous derivatives with respect 
to x up to the (n — 2)th order inclusive in the squarea <x,€ <b. 

3°. For each fixed value of € in the interval (a,b), the function G(x,€) has con- 
tinuous derivatives with respect to x of the (n — 1)th and nth orders respectively 
in each of the intervals [a,é) and (€,b]. The derivative of the (n — 1)th order 
makes a jump in value, [P,(§)]~1, at x = €: Le. 


on-1 n-l 


0 _ 
sani GE + 0,8) — gaa GE — 0,6) = [Po(O)I-?, 


4°. For a fixed € in the interval (a,b), the function G(x,€), regarded as a function 
of x, satisfies the boundary conditions U,(G) = 0,v = 1,2,...,n, and, in each of 
the intervals [a,é) and (€,b], the operator equation (G) = 0, Iy) and Uy) = 0 
being the differential expression and boundary conditions which generate the 
operator L. 


Repeating the arguments of §3.3-3.5, we obtain the following results: 


THEOREM | 
If the boundary-value problem Ly = 0 has only the trivial solution, then the 
‘operator L has just one Green’s function. 


THEOREM 2 
If the boundary-value problem Ly = 0 has only the trivial solution, then the 
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equation Ly =f, where f is any continuous vector function f(x), has just one 
solution, which is given by the formula’ 


v(x) = | OO fOae. (1) 


Theorem 2 means that L~? is an integral operator with the kernel G(x, €). 

The kernel G*(x,é) = [G(€,x)]* is called the adjoint of the kernel G(x,€). A 
kernel is said to be Hermitian if it is the same as its adjoint: G(x,é) = [G(é,x)].* 
Theorems III and IV of §3.5 and all the theorems in §3.6 remain valid. 


2. Green’s Function for the Operator L — A1.? 


We seek an expression for the Green’s function for the operator L — Al. Let 
Y, = Y,(x,A) be a fundamental system of solutions of the equation ((Y) = AY. 
We choose this system so that the Y,(x,A) are integral, analytic, matrix- 
functions of the parameter A (it is always possible to do this). We write |W| for 
the determinant of the matrix 


y“ -1) ya -1) 
yy") 1. yu® 

VW = ; (2) 
Y, Y,, 


and denote by V,, v = 1,2,..., n, the mth-order matrices consisting of the 
cofactors of the elements Y, in the determinant |W]. Further, we let W,, 
v = 1,2,...,n denote the transposed matrices of V,. 


We put 
] 
Z, = Ww] W,, (3) 
LY! ¥,(0)Z,(6) for € <x, 
g(x) ={ (4) 
—4) Y\x)Z,(6 for é > x, 
U,( Y,) vee U,( Y,) Wi. Wr 
U = ,U = 
U(Y,) ... U,(Y,) Wa... W,, 


Here the W,, are m X m matrices. 


1 For a fixed x and &, P(x,£)f() means the result of applying the operator P(x,é) to the 
vector f(&). 


2 In this section 7.2, the author assumes that A is not an eigenvalue. 
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I The Green’s function G(x,§,A) for the operator L — Al is given by the formula 


GOxEA) = a(x) — YYW ,U(g). (5) 


J,;v= 


The proof of formula (5) is similar to the proof of formula (34) in §3.7. If 
we apply the method of variation of the parameters to the equation 
Ky) — Ay =f, we obtain for its solution y the formula 


y= Vode, + f? ecogayOae, (6 


the analogue of formula (31) in §3.7; the c, here, of course, are constant 
vectors. By substituting this expression for y into the boundary conditions 
U,(y) = 0,» = 1,2,...,n, we find these vectors c,; then by substituting the 
expressions for c, into (6), we arrive at formula (5). 

If we write A for the determinant of the matrix U, then 


] 
W iv —_ A V vs (7) 


where V,, is the transpose of that mth-order matrix consisting of the cofactor 
of the element U,( Y,) in the determinant A. Formula (5) can then be rewritten 
as 


GlnEd) = BOE) — 5D VOM WV). (8) 


3. The Analytic Nature of the Green’s Function for the Operator L — Al 


The functions Y,(x) and consequently also g(x,€,A), V,,,, as well as U,(g) and 
A are integral, analytic functions of the parameter A. Hence it follows from (8) 


that: 


II. The Green’s function for the operator L is a meromorphic matrix-function of 
the parameter X, and only eigenvalues of the operator L can be poles of this 
function. 

We want to look at the case when A, is a simple zero of the function A in more 
detail. 

We assume that the adjoint operator exists as a differential operator of the 
type under consideration here, and let z(x) be one of its eigenfunctions corres- 
ponding to the eigenvalue 29; we denote by y(x)z*(&) the matrix composed of 
the elements y,(x)z,(€). Using the result in §7.2, and repeating the argument 
in §3.8, we can deduce: 
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III. At a simple zero Xo of the function A(A), 
y(x)z*(£) 

(X — Ao) [° (,2z)dx 


where G,(x,€,d) is regular in the neighbourhood of the point Xo. If, in particular, 


[,2)dx = 1, 


G(x,€,A) = + G,(x,€,A), (9) 


then 


Gx) = — FE 2 + Glee). (10) 


For differential operators of the first order in the space of vector functions, 
the Green’s function was the subject of investigation in papers by Bunitski 
[13], Bocher [11], Birkhoff and Langer [6], and Bliss [8]. 


$8. Asymptotic Behaviour of the Eigenvalues of a Differential Operator 


1. Statement of the Problem 


In §4 we examined the asymptotic behaviour of eigenvalues and eigenfunctions 
of a differential operator in the space of scalar functions. The results obtained 
enabled us to derive in §5 some general theorems on expansions in terms of 
eigenfunctions. In this section we shall be concerned with the analogous 
problem, the asymptotic behaviour of the eigenvalues and eigenfunctions of 
a differential operator, in the space of vector functions. 

To solve this problem, we first investigate the asymptotic behaviour of 
solutions of the matrix equation 


KY) + p"Y =0 (1) 


for large values of |p|. The decomposition of the complex p-plane into the 
domains S, or T, (see §4.2) plays an important part here. 

For simplicity of exposition we discuss the case where P(x) = 1, and there- 
fore take as our differential expression 


KY) = YOLP\Y" V+... +P,Y. (2) 


(The possible generalizations to the case P)(x) # | still remain open for con- 
sideration). Further, without loss of generality, we may assume that P,(x) = 0. 
For, otherwise, we can apply to equation (1) the substitution Y = UY, where 
U is a solution of the matrix-equation 


J 
U'+ PiU =0 


satisfying the condition 
det U # 0. 
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We then obtain for Yan equation of the form 
KY) + p'Y=0 


in which the coefficient of Y“~» vanishes. 


Finally, we can assume a = 0 and b = 1, and so we are considering the 
interval [0,1]. (See §4.1). 


2. Asymptotic Behaviour of the Solution of the Matrix Equation /(Y) + e"Y = 0 
for large |o| 


Let then 
KY)= YO P,Y"-9 +... 4+ PLY. (3) 
We consider the matrix-equation 
KY) + p"Y =0. (4) 
By applying to this equation essentially the same considerations as in §4.5 we 
obtain the following theorem. 
THEOREM | 


If the matrix-functions P.,,...,P,, are continuous in the interval [0,1],1 then the 
matrix-equation 


YOu P,Y"-I4 ..4P,¥+p"Y =0 


has in each domain T of the complex p-plane n linearly independent solutions 
Y,,Yo,.-.,¥, which are analytic in p€T and which for sufficiently large |p|, 
péT, satisfy the relations* 


Y,= coon + (5). 
pP 


dY, sans ; 
ae ~ Pe |oat + 0(,) 


d"-*Y t ¢ eed OX n- l 
pena = Pte E at o(;)|: 


(5) 


A(x, 
Here (5) is a matrix of the form ae p) with A(x,p) uniformly bounded: 
pP 


A,jcp)| <M for |p| > R, 0 <x <1, 


where M and R are independent of x. 


1 See the footnote, p. 48. | 
2 Here 1 stands for the mth-order unit matrix. 
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If the functions P,(x) are not only continuous but also sufficiently often 
differentiable, then the asymptotic formulae (5) can be made more precise, and 
formulae of the form \ 


_ pa,x Y Via) 4 fa) Y a) J 
Y (x) =e 1B a re + O(- | 


can be obtained, with similar formulae for the derivatives Y,v = 1,2,...,.2—1. 
(See §4.6 for the details). 


3. Normalization of the Boundary Conditions 


A given differential operator is characterised by the boundary conditions 
Uy) = ~~ 0, y= = I, 2,.. 
The number k is salted the order of a form U(y) 1f U(y) contains at least one 
of the vectors y and y™ but does not contain the vectors y$” or y™ forv > k. 
We consider forms U(y) of order (n — 1), if there are any; they have the 
form 


Uy) — Ay nie" 0 + By, yr" ) +... 


The rectangular matrix [A,,-1,B8,,-1] has m rows and 2m columns. The 
maximum number of linearly independent rows of 2m elements is, however, 
2m; if, then, we replace the rows in the forms of order (n — 1) by linear com- 
binations of these rows (if this process is necessary), we can arrange that not 
more than two forms of order (nm — 1) occur. 

Continuing in the same way with the remaining forms, we can, after a finite 
number of such steps, reduce the boundary conditions to the form 


U,(y) — Uo(¥) + Uy(y) — 0, (6) 


where 


ky -1 
U,0(y) —= Aye + d, A, V9, 

ont (7) 
Un) — By + d By jy, 

j= 


n—]1 2k Sk >... > k,, 70, kyae > ky, 


and where, foreach vy, v = 1,2,...,n, at least one of the matrices A,,B, is different 
from the zero-matrix. 

The operations just described are referred to as the normalization of the 
boundary conditions, and the finally resulting boundary conditions of the 
form (7) are called normalized boundary conditions. 
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4. Regular Boundary Conditions 
We consider a fixed domain S,, and number a ,w.,...,w, so that, for pe S, 
B(pa,) < Apa) <... < A(pw,). 


The asymptotic formulae which follow are going to be derived for a 
particular class of boundary conditions which we shall call regular. The 
definition of regular boundary conditions depends on whether 7 is even or odd. 
(a). Suppose n is odd; n = 2n — 1. 

The normalized boundary conditions (6) are said to be regular if both the 
numbers 9, and @,, defined by the equation 
+ 45+... 46,5" = 


k k k k k 

A,wy! Aywty (A, -t. SBy)w,? Byory eee Bw, 

A,wi? A,wtt, (A, + 5B,)w* Bow" Bow” 

QWy 204 2 TT SDo)W,, ages 2Wp, 
(8 a) 

kn kn kn kn Kp, 


do not vanish. 


(b). Suppose n is even; n = 2. 
The normalized boundary conditions (6) are said to be regular if both the 
numbers 6_,, and 6,, defined by the equation 


6_ 5" + Os FP st +O" = 


1 
Aywy ... Aywrt, (Ay + SB)wi (4, +2 Blo, Bywtg Brow 


] 
A,w,? cee A,wit_, (A, -- sB)wi fe -+- ~By)of Bow'*, 9 eee Bowe 


| e eee ° ry ° o ooe e 


l 
A,win _ Aon y (A, + sB,)wy" (4, +- ; B, jo Bon _ B wir 


do not vanish. 


We can see at once that regularity so defined does not depend on the particu- 
lar domain S, selected. For, the numbers 4 and @, (or @_,, and 6,,) are 
determinants, which are obtained from those used in §4 to define the numbers 
6, merely by replacing a, and 8, by the matrices A, and B,. Hence precisely the 
same laws hold for their transformation due to a change in domain S;; in 
particular, for odd n, they are multiplied by a factor of modulus 1. For n even, 
the whole of the calculation carried out in §4.8 can be analogously repeated. 
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It follows, therefore, that here too the roots of the equation 
O_s7™ +A yo} toe + In S™ = 0 


transform into their reciprocals when the domain is changed from Sp to S, or 
S,,-1. For odd n, there are two classes of equations 


05 + Os +... + 6,5" = 0, 


with the same roots in each case in each class. For S, with & odd, the roots are 
é)  €®- and for S, with k even, the roots are &”,...,¢0. 


ni 9 


5. Asymptotic Behaviour of the Eigenvalues 


In questions concerning the asymptotic behaviour of the eigenvalues, the 
equation 


6, + Os +... + 9,5" = 0 (9 a) 
for odd n, and the equation 
6_ Ss + Oma TE tw. +O” = 0 (9 b) 


for even n, both play an important part; the numbers @, here are the numbers 
defined in the previous subsection 4. 

If the boundary conditions are regular, these equations are respectively of 
order mand 2m, and all their roots are different from zero. 

In the following theorem we assume (though for brevity we do not mention 
this again in the hypothesis) that the coefficients of the differential expression 
considered are continuous matrix-functions in the interval [0,1]. It should, 
however, be particularly mentioned that all the results of this subsection 
would still hold if we required only that these coefficients be arbitrary, sum- 
mable, matrix-functions in the interval [0,1]. 


THEOREM 2 


Let L be a differential operator of the nth order, defined in the interval [0,1], 
whose differential expression contains no derivative of the (n — 1)th order, and 
whose boundary conditions are regular. 


(1). If nis odd, then to each simple root £5 corresponds a sequence X, and to 
each simple root &? corresponds a sequence X° ), of eigenvalues of the operator 


L, and 


1) 
19 = cFakwip| ve MA? 4 of 1] 


Nn | (2) 
AQ) = (4:2ka!) [ 4" nt + (2) 


k = N,N + 1),..., (10a) 


where the upper or lower signs hold according asn = 4v + 1 orn = 4v — 1. 
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To an r-fold root & or €®, on the other hand, correspond in general r 
Sequences of eigenvalues x) or XY) which for large values of k satisfy the 
following formulae: 


Ing 5? l 
Ne — (F 2kniy | + ao es S) + (ja sn) 


2kai J =jos(Jo + 1),....CJo +r— 1), 
» caf, , MI €® 1 \J/k=NMN+4)),.... 9 (la) 
Nej — (+:2k7i) E + a’ O jit 
Here it may be taken that &) =). =... =, 4. 


The same rule for upper or lower signs holds as for formulae (10a). These 
eigenvalues MY =Jo(Jo + 1),.-.,4o + 7 — 1), may coincide for certain j and 
k to give multiple eigenvalues. 

(2). If is even, then to each siniple root € of the equation (9b) for the domain S, 
corresponds a sequence A, of eigenvalues of the operator L, and 


n Ino € 
2kri 


\n J 
A, = (2kni) E - + (23) | k = N,N + 1),..., (10 b) 
where the upper or lower signs hold according as n = 4v or n = 4v +2. 
To each multiple zero € of equation (9b), with multiplicity r, correspond r 
sequences of eigenvalues \,,; of the operator L, and 


n Ino €é 
2kri 


Mj = (keri) + + (zr) |p f= h2 sr k= MN Dyes (11b) 


and again the sign is — or + according as n = 4v orn = 4v +2. 


For sufficiently large ||, the operator L has no other eigenvalues and the 
multiplicity of any particular eigenvalue X is equal to the number of eigenvalues 
from the above sequences which happen to coincide there, i.e. the multiplicity 
is not greater than the multiplicity of the corresponding root & of equation (9). 


Proof. First let n be odd, sayn = 2u — 1. 
We consider a fixed domain 7; let the numbers w ,w»,...,w, be numbered in 


sequence so that, for p € T, 
R(p + c)oy) < A((p + c)w2) <... <A((p + c)w,). 


By applying formula (5) of subsection 2, and repeating the argument in the 
proof of Theorem 2 in §4.9, we obtain 


U,( Y;) = (pw ;)[A,] for] < by, 
U,(¥;) = (pw,)e?°[B,] for j > pw, (12) 
UY.) = (pa, tTAy] + 2? “(BT 


124 $9 EXPANSION IN EIGENFUNCTIONS 


We substitute all these expressions into the equation A = 0; it follows from 
the definition (8a) of the numbers @, that, after removing the factors 
pp")... pkn ere tt. e?, this equation takes the form 


[60] + [OuJe@? + [Bale*? +... + [B,]e"™? = 0. 


Hence, and by arguments similar to those used in proving Theorem 2 of §4.9, 
follow all the assertions of the Theorem, for odd. 

The case of n even is treated in similar fashion. 

By using Theorems 1 and 2 we can obtain asymptotic formulae for the 
eigenfunctions, but we shall not go into this question in detail. 


§9. Expansion in Terms of Eigenfunctions of a Differential Operator 


1. The Case of a Self-Adjoint Differential Operator 


Let L be a self-adjoint differential operator in the space of vector-functions. By 
repeating the argument of §5.2, we obtain the following theorem. 


THEOREM 1 


Any vector-function in the domain of definition of a self-adjoint operator L can 
be expanded in a uniformly convergent, generalized Fourier series in terms of 
the eigenfunctions of L. 

If y,(x),e(x),... Is an orthonormal system of eigenfunctions of the operator 
L, so that all eigenfunctions of LZ for arbitrary eigenvalues can be expressed as 
linear combinations of these vectors, then the expansion referred to in Theorem 
1 has the form 


o@) 


SO) = 2 nal) (1) 
with 
en = |, Learn(x))dx. (2) 


We recall once more that f(x) and y,(x) are vector-functions; if, then, we put 


f(x) —= {f1(x), fo(x),- ae Sn(X)}s 
Vrlx) — {Vn1(X)VnolX)>--- Yam} a= 1,2,3,. a) 


then the representation (1) is equivalent to the system of expansions in 
uniformly convergent series 


10 6) 


f (x) = dy, CrVnjX)s J — 1,2,...,71, 


a= 


the same coefficients c, occurring in each of the series. 
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It follows from Theorem 1 that the functions y,(x),yo(x),y3(x),... form a 
complete orthonormal system in the space L?(a,b) of all measurable vector 
functions f(x) which satisfy the condition 


[fax < + 0. 


[| f(x)|? means here the scalar product of the vector f(x) with itself.] 
It follows that 


fi \feolax = ¥ heal 3) 
with ; 


on =f, Ped.r(oae. 


2. Expansion in terms of the Eigenfunctions of a Differential Operator Generated 
by Regular Boundary Conditions 


The considerations in §9.1 are not applicable if Lis not a self-adjoint operator. 

So we now make use of a different method which depends on the analytic 
properties of the Green’s function for the operator ZL — Al and on the 
asymptotic formula obtained in §8. Here we can assume that the boundary 
conditions which generate L are regular; further, without loss of generality, 
we may assume that the boundary-value problem Ly = 0 has only the trivial 
solution. Then ZL has a Green’s function G(x,§). As in §5.3, we consider a 
sequence of circles [',, k = 1,2,..., with centre at the origin, and having the 
following properties: 


1°. The radius R, of the circle I’, increases without limit as k ~ oo. 


2°. There is a positive number 6 such that all eigenvalues of the operator L 
have a distance >6 from each of the circles [’,. 

By virtue of the asymptotic properties of these eigenvalues proved in §8.5, 
such circles do exist. 

Let G(x,é,A) be the Green’s function for the operator ZL — Al, and, in 
particular, let G(x,€,0) = G(x,é) be the Green’s function for the operator L. 
We consider the integral 


I G(x,€,A) 
be = Fi rm ( (4) 
LEMMA 
On the circles T,, the function G(x,€,A) satisfies the inequality 
M 
IG j(2,8,A)| < rt (5) 


where M is a fixed constant, uniformly for k = 1,2,3,...anda <x,€ <6. 


126 §9 EXPANSION IN EIGENFUNCTIONS 


Proof. If we put A = —p", then, for a suitable choice of arg p, the circle I’, is 
transformed into a circular arc y, subtending an angle 27/n at the centre, and 
vy, goes through two neighbouring domains §),5, of the complex p-plane. We 
want to estimate the value of the function G(x,£,A) on the arc ,. 


Let the numbers w,,w.,...,w, be numbered in sequence so that, for p € S}, 
(pay) <A(pa) <... <A(pu,). 


We use the formula for G(x,€,A) from §7.2, and substitute therein the asymptotic 
expressions for Y,(x,A) found in §8.2. We get 


Z(E) — e7 Pers a [—w,1], (6) 
np 
and so, by (4) in §7.2, 
_ | pay »| 
e(x,f) = F Inphat (> e? [1], (7) 


where the — sign is to be taken if £ < x, and the + signif é > x. 
For the rest of the proof itis best to deal with m odd and n even separately. 
First let n be odd; n = 24 — 1. Then, for p € S,, we have 


Apo) < 0,..., A( pw, -1) < 0, (8) 
A(pw,, +1) > 0,..., A(pw,) > 0. 


We denote by y, that part of the arcy, which lies in S, and on which KA(pw,) 
> 0, and by y, the part of y, which lies in S, and on which A(pw,) <0. We 
want to estimate the value of the function G(x,€,A) on the arc y,. By using 
formula (12) of §8.5, it can easily be shown that, on y,, 


YO), = O(,) 9) 


To obtain an estimate of G(x,,A) on y,, we first transform the expression for 
G(x,g,A) in formula (5) of §7.2. By virtue of formula (4) in §7.2, for g(x,€) we 
have the equations 


20.) =), ¥,0)Z.0), KL =4+D YNZ; 


consequently 


Ula) = 4 Y Walt) — UW VIZ. 
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We write this expression in the form 


Ug) = 4 Y[- UW.) + 2U (VIZ) +4 Y (WY) ~2U VIZ 


a=n+1 


-4 * U(YIZAE) + Y Un VOZ8 


bY YOOZO— YL UFZ. (10) 
By (6), we have, however, fora <p, 


ky 
Ui ( Y)Z,(E) = — orl wy” r 1 pp@alt 7 OB] 3 


hence, and by (8), it follows that 


U,(Y,)Z,(é) = (5 | 1) 


Similarly, for a > pu, we get 
I] 
Ul Ye)Z(8) = Of 57-1) 


So (10) can be written in the form 


U,(g) — —4 y U,( Y)Z,(€) + 5 ») ; U,( Y)ZAE) + («-i=): (11) 


If this representation of U,(g) is now substituted into formula (5), §7.2, for 
G(x,§,A), we obtain 


Lu n l 
GwEd) = abd) F4¥ HIZO — 4 VVC + O(Fa), 
j= J=ur 
since, by definition of the inverse matrix, 


Ld WU Y,) — Ojyl, 


] ] 


Combining (12) with (4) from §7.2, we arrive at the formulae 


and further, by (9), 


Dy YZ IE) + O( ;) for & <a, 
G(x,8,A) = 0) 
_ >» Y (x)Z,(£) — 043] for €> x. 
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By (8) and the asymptotic formulae for Y (x) and Z (8), it follows that on the 
arc y, for € < x, 


_ I pope $) — i 
ES VZAO = — pays Yoder? = O( sa), 


j= 


and for € > x, 


YYZ) =— AY [ole = 0( 5). 


joutl J=ut+1 


In conjunction with (13) this gives the desired relation 


occa) = of) 


on the arc y,. 

Similarly we can verify this relation on the arc y,, and on that part of the 
arc y, which lies in the domain S4. 

This proves the lemma for the case of n odd. 

For n even, the proof runs on exactly similar lines. 

Using these Lemmata, the following theorems may be deduced exactly as 
in §5.3: 


THEOREM 2 


The Green’s function G(x,€) of a differential operator which is generated by 
regular boundary conditions can be represented as a Series 


= ee 


G(x,£) = 


which, at least for a certain definite order of its terms, converges uniformly. 
Here H,(x,€) denotes the residue of the function G(x,€,A) at the pole X,. 


THEOREM 3 


[f all the eigenvalues of an operator L generated by regular boundary conditions 
are simple zeros of the function A(A), then its Green’s function G(x,€) can be 
represented as a Series 

5 be ye eeder(e) 


y=l Xr, 


G(x,¢) = 


which, at least for a certain definite order of its terms, converges uniformly. 
Here y,(x) and z,(x) are the eigenfunctions, of the operators L and L* 
respectively, associated with the eigenvalues , and i,; and those eigenfunctions 
may both be so normalized that 


[, 0@.2.00)dx = 1. 
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Let L be an operator which is generated by regular conditions and of which all 
the eigenvalues are simple zeros of the function A. Then any function f(x) in the 
domain of definition of the operator L can be expanded as a Series in eigen- 
functions of L: 


fl) =. er4(a), where c, = [* MO. 


The series converges uniformly, at least for one particular order of its ternis. 
Here y,(x) and z,(x) are eigenfunctions of L and L*, for the eigenvalues i, and 
A, respectively; these eigenfunctions may be normalized so that 


| (y,(x),2z,(x))dx = 1. 


For simplicity we have restricted ourselves to the ordinary eigenvalue- 
problem and have not discussed the generalized problem Ly = AMy at all, nor 
the problem Ly = AB(x)y in particular. For the latter problem, if certain 
restrictions are imposed on B(x), it can be assumed that B(x) ts a diagonal 
matrix which has p4(x),po(x),...,4,(x) as its diagonal elements, the »,(x) being 
continuous and all different (see $6.8). We can obtain, for the solution of the 
corresponding differential equation and for the eigenvalues, certain asymptotic 
formula in which 


exp f [, le (atl 


takes over the réle played by exp {pw (x — a)} previously. These formulae can 
be used to obtain expansions in terms of eigenfunctions. It may be mentioned 
that no one has investigated up to now either this case or the general case 
Ly = AMy in any detail (except for the special cases n = 1 or m = 1). 


2. 


3. 
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